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Many dissociation reactions are accompanied by a change in the volume 
of the gas mixture. For instance, during the dissociation of carbon dioxide 
gas, there is an increase in the number of moles, and consequently an 
increase in the volume of gas mixture, by 4 mole for each mole of fully dis- 
sociated carbon dioxide gas. 
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For dissociation reactions which take place with an increase in the num- 
ber of moles, the composition of the gas mixture will depend upon the 
pressure to which it is subjected. The increase in pressure leads to the sup- 
pression of dissociation reactions and increases the content of the products 
of complete combustion in the gas mixture. In other words, increased pres- 
sure reduces the degree of gas dissociation if it is accompanied by an in- 
crease in the number of moles. For dissociation reactions which take place 
without increase in the number of moles the composition of the products 
is independent of pressure. Since the majority of dissociation reactions of 
rocket fuel combustion products take place with an increase in the number 
of moles, the increase in the combustion pressure decreases the degree of 
dissociation somewhat (but not significantly). 


3. COMPOSITION AND TEMPERATURE OF COMBUSTION PRODUCTS IN 
Rocket Morors. COMPLETENESS OF CHEMICAL ENERGY RELEASE 


Ordinary rocket motor fuels are composed of four elements: carbon, 
hydrogen, oxygen, and nitrogen. 

If the combustion was not accompanied by dissociation, the products 
of combustion would consist of carbon dioxide, CO.; water vapor, H.0O; 
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Fig. 5.20. Dissociation of water vapor, H.O, and carbon dioxide, CO», with tem- 
perature. 


and molecular nitrogen, N2.* However, even at moderate temperatures— 
about 2800°K—considerable dissociation of carbon dioxide and water 
vapors may be observed (Fig. 5.20). In addition, carbon monoxide, CO; 


* For the stoichiometric relations of components. 
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hydroxy] radical, OH; molecular oxygen, O2; and hydrogen, H2; are formed. 
At a still higher temperature, the presence in combustion products of nitro- 
gen oxide, NO; and gases in an atomic state: hydrogen, H; oxygen, O; and 
nitrogen, N; also becomes significant. 

The composition of the combustion products in the chamber is deter- 
mined by the equilibrium constants corresponding to the dissociation 
reactions, taking into account the pressure in the chamber, pe. Aside from 
that, the composition of gases depends, obviously, on the relative content 
of various elements in the fuel. 

Since the composition of combustion products depends on temperatures, 
the solution of Eq. (5.20) or (5.21) requires a considerable amount of cal- 
culation. Ordinarily, it is necessary to assume a temperature, then to find 
the composition of the products of combustion, i.e., partial gas pressures, 
according to the composition of fuel and equilibrium constants, and, finally, 
to check the basic equation of combustion 


Eg = Ng. 


The actual temperature and its corresponding composition of combus- 
tion products is determined, therefore, by a method of selection. 

Calculation results show that because of the sharp increase in the degree 
of dissociation of the combustion products with increase in temperature 
above 3000°K, increase in the availability of chemical energy in the fuel 
does not yield a proportional increase of temperature in the chamber. 
For instance, combustion temperatures of two fuels may be compared: 
kerosene plus nitric acid, and kerosene plus oxygen. The calorific value of 
the second fuel is 2400 kcal/kg, which is approximately 70% higher than 
the calorific value of the first (1400 kcal/kg). However, the combustion 
temperature of kerosene plus oxygen fuel (8550°K) is only 15% higher 
than the combustion temperature of kerosene with nitric acid (8050°K). 
This is a direct result of the intensive dissociation of the combustion prod- 
ucts and reduction in the completeness of the release of chemical energy. 

Let us further examine the dependence of combustion temperature on 
pressure. 

This dependence is well illustrated by a graph in Fig. 5.21, drawn for 
the combustion products of the fuel: kerosene plus liquid oxygen with 
a = 0.7. From this dependence it may be seen that, initially, at low values, 
the increase in pressure leads to a rapid increase in temperature, and that 
subsequent increase in pressure leads to a less rapid increase in temperature. 

The temperature and the composition of combustion products, besides 
being affected by pressure, are also affected by the relationship of com- 
ponents within the fuel. 

Not taking dissociation into account, the maximum combustion tem- 
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perature and the maximum quantity of heat released during combustion 
would occur when the theoretical relationship of combustible and oxidizer 
1S vo 1.e., When a = 1. However, the dissociation phenomenon reduces the 
combustion temperature, and the unequal stability of combustion products 
with respect to dissociation results in the fact that the maximum tempera- 
ture and maximum heat release correspond to values of a < 1 in contem- 
porary liquid rocket fuels, i.e., fuels which have a deficiency of oxidizer and 
an excess of combustible. 

Figs. 5.22 and 5.23 show the dependence of the combustion temperature 
and the quantity of released heat, Q, on the oxidizer excess coefficient a 
and pressure for kerosene plus oxygen fuel. The shift of maximum tempera- 
ture and heat release in the direction of a < 1 is explained by the fact that 
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Fia, 5.21. Dependence of oxygen plus Fic. 5.22. Dependence of oxygen plus 
kerosene fuel combustion products tem- kerosene fuel combustion temperature on 
perature on pressure when a = 0.7. the oxidizer excess coefficient a and pres- 


sure p. 


in the presence of insufficient oxidizer the relative content of carbon mon- 
oxide, stable with respect to dissociation, is increased in the combustion 
products. As the curves of Fig. 5.23 show, the loss of heat due to dissocia- 
tion in the combustion chamber is quite high—at a pressure of 14.7 lb/in.? 
it constitutes more than 30% of the calorific value of fuel whose heat con- 
tent is 2400 kcal/kg. 

We have examined the effect of dissociation using as an example a fuel 
which has the maximum value of heat content. Inasmuch as the intensity 
of dissociation is reduced at lowered temperatures, the influence of disso- 
ciation is reduced for fuels with lower calorific value (oxygen plus alcohol 
or nitric acid plus kerosene), even though it still remains considerable. 
For instance, the loss of heat in the chamber due to dissociation for the 
oxygen plus alcohol fuel constitutes 18-25% of the calorific value, and for 
nitric oxide fuels it is 12-18%. 
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The composition of combustion products determines the value of their 
gas constant. The simultaneous influence of the composition and tempera- 
ture on the heat capacity of gases leads to the corresponding change in the 
adiabatic exponent k. The influence of temperature and pressure of the 
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Fig. 5.23. Dependence of the quantity Fic. 5.24. Dependence of combustion 
of released heat during the combustion of | products gas constant R and adiabatic 
oxygen and kerosene fuel on oxidizer excess exponent k on temperature T and pressure 
coefficient œ and pressure p. p. 


combustion products on the value of gas constant R and adiabatic exponent 
k is illustrated by the graph in Fig. 5.24. As is shown in this graph, the gas 
constant of the combustion products increases with increasing temperature. 
This happens because, as a result of dissociation, there is a reduction in the 
content of polyatomic gases in the products of combustion. For the same 
reason, the adiabatic exponent of the products of combustion increases 
with an increase in temperature and a decrease in pressure. 


VI. Flow of the Combustion Products through the Nozzle 
of a Rocket Motor 


A. Basic Laws of Gas Flow 


1. Gas Ftow PARAMETERS 


The products of combustion, formed in the chamber of the motor, flow 
into the nozzle, where they are expanded and accelerated. In the process of 
motion along the nozzle the flow parameters change. Parameters of gas 
flow in the nozzle, aside from the parameters of state p, v (or p), and T, 
mentioned earlier, also include velocity w, and the area of the stream cross 
section S. 

During the flow along the nozzle, reduction in gas temperature and 
pressure take place. Reduction in temperature leads to the reduction in the 
degree of dissociation of a gas which was strongly dissociated in the cham- 
ber—recombination of atoms and radicals into molecules. 

The reduction in pressure does not hinder this process significantly. It 
is clear that the recombination results in an additional (with respect to the 
combustion chamber) release of heat and facilitates more complete con- 
version of the fuel chemical energy into kinetic energy of the gas flow emerg- 
ing from the nozzle. 

For the time being, let us not take into account the recombination 
phenomenon, and let us examine the flow of a gas of constant composition 
along a nozzle. As a basic equation for determining the parameters of flow, 
let us use the equation of the gas state 


p/p = gkT (6.1) 
or, in another form, 
pv = RT. 


The second equation, which contains the parameters of the gas flow, is 
the equation of the thermodynamic process. We shall dwell on this equa- 
tion in some detail. 

The change in the gas state during thermodynamic processes can take 
place in various ways; for instance, at constant volume, at constant pres- 
sure, or at constant temperature. Depending on this the interrelation 
between gas parameters varies. 
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The more general form of the thermodynamic process equation is 


p/p” = constant 


or 


pu” = constant. (6.2) 


In the future we will agree to discuss only those processes whose exponent, 
n, remains constant for the entire process. 

Taking various values of the exponent n, it is possible to describe basic 
thermodynamic processes which take place in a gas. For instance, taking 
n = 0, we will get p = constant. Consequently, Eq. (6.2) will express the 
equation of an isobaric process. Such a process takes place, for instance, in 
the combustion chamber of a rocket motor. This process in pv coordinates 
is shown in Fig. 6.1 by a line I (isobar). 


2 





Fig. 6.1. Curves of the basic gas thermodynamic processes. 


When n = 1, we will get pv = constant or, taking into account the 
equation of state, T = constant, i.e., the equation of an isothermal process 
(line II in Fig. 6.1). 

When n —> œ, Eq. (6.2) represents the equation of an isochoric process, 
v = constant. 

Line III in Fig. 6.1 (isochor) divides two characteristic areas: the area 
of expansion processes (increase in the specific volume v), and the area of 
compression processes (reduction in specific volumes v). 

Of all the possible gas processes, we will be most interested in the 
process of the expansion of a gas of constant composition taking place 
without heat exchange between the gas and the surrounding medium. It 
would seem that such conditions for gas flow, especially in the nozzle of a 
rocket motor, cannot exist, since the gas is at a high temperature and comes 
in contact with rapidly cooled walls of the nozzle. However, in actuality 
this is not the case. The time during which the gas is in contact with nozzle 


BASIC LAWS OF GAS FLOW 163 


walls is very short. Due to the high flow velocity, this time, even in motors 
having considerable dimensions, represents a magnitude in the order of 
0.001 sec. Besides, a considerable portion of the gas passes at some distance 
from the walls and releases heat to them only to a small degree. 

For the above reasons the amount of heat given up by the gas to the 
walls is insignificant by comparison with its total supply, and the process 
of the gas expansion in the nozzle of a rocket motor may be considered as 
taking place without heat exchange to the walls. Such an expansion process 
is called adiabatic. 

An adiabatic process taking place in the gas is characterized by the 
manifestation of the law of the conservation of energy in its simplest form. 

Reduction in gas temperature and pressure takes place during adiabatic 
expansion. Concurrent with the reduction in temperature there is a reduc- 
tion in the internal energy of the gas. The reduction in pressure also causes 
a reduction in the potential energy due to gas pressure. According to the 
law of conservation of energy, the difference in energies which the gas pos- 
sessed at the beginning and at the end of the process is entirely converted 
into the work of expansion of the gas. 

The work of gas expansion is variously exploited in different machines. 
In the piston engine, for instance, it is converted into the work of displacing 
a piston; in reaction motors, into kinetic energy of a gas flow. The amount 
of work obtained during adiabatic gas expansion is easily calculated by 
determining the change in the internal and potential energies of the gas 
during the expansion process. This we shall do below during the derivation 
of the energy equation. 

The exponent n for an adiabatic process is equal to the adiabatic 
exponent 


n= kh = Cy/cy. 


As can be seen from Fig. 6.1, the expansion curve with an exponent, 
n =k > 1 (line IV), passes between the isotherm and isochor but with a 
steeper slope than the isotherm. 


2. STABLE AND UNSTABLE Gas FLOW 


It is not sufficient, in determining the relationship between the param- 
eters of a moving gas, to use only the thermodynamic relationships. 

Of all the possible flows, let us first emphasize the stable or steady 
state flow. Steady state flow is that flow which is characterized by unchang- 
ing gas parameters (velocity, pressure, temperature, density) at any point 
in the gas stream. If these parameters are not constant with respect to time, 
then the flow is considered unstable. 

In many technical problems involving motion of gases, the gas flow may 
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be considered as stable, even when it is not entirely stable, at least as an 
average. This considerably simplifies the solution of many practical 
problems. 

It is this type of flow which we will consider as taking place in rocket 
motor nozzles. At the same time, the period of starting and stopping the 
motor or the time of transition from one working regime to another, when 
the gas rate of discharge varies, will be excluded from this examination, 
since such changes vary the parameters of gas flow. 

An example of an obviously unstable flow is the gas flow in the pulsating 
ducted motor. In this case, gas parameters vary not only along the length 
of the motor but also with respect to time. 


3. VELOCITY DISTRIBUTION IN THE STREAM Cross SECTION. 
ONE-DIMENSIONAL FLOW 


Gases, like liquids, have viscosity. A characteristic of viscosity is the 
coefficient of viscosity u, determined by the relationship 


tT = p(Aw/Ay), 


where 7 is the sheer force related to a unit area, arising between two layers 
of gas or liquid moving parallel with respect to each other and located at a 
distance Ay from each other; and Aw is the difference in velocities between 
these layers. 

The gas viscosity is considerably less than the viscosity of liquids and 
becomes apparent only when there is a large difference in velocities; for 
instance, immediately at the surface of a body placed in the gas stream. 
The particles of gas come in intimate contact with the surface and would 
seem almost to attach themselves to it and remain immobile. 

As the distance from the surface of the body increases, the velocity of 
gas particles increases rapidly, and then remains constant and equal to the 
velocity of flow. 

The thin layer of gas in which the velocity change from zero to the 
velocity of the free stream takes place is called the boundary layer. The 
gradual increase in velocity within the boundary layer is explained by the 
action of internal friction forces within the gas. 

The thickness of the boundary layer at the surface of the body in the 
gas stream increases in the direction of flow, and with great linear dimen- 
sions of the body can reach large magnitudes. 

The motion of gases and liquids, as is known from the course in physics, 
may be laminar and turbulent. In laminar—stream line—flow, streams of 
gas do not mix, whereas during turbulent flow there is an intensive mixing 
of fluid volumes. 
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Laminar flow may be either stable or unstable. Turbulent flow is always 
unstable. However, since the mixing of fluid layers takes place within vol- 
umes considerably smaller than the over-all flow dimensions, the turbulent 
flow may be considered, on the average, as being stable. 

During gas flow along a wall in the boundary layer, the flow is laminar 
for some distance, and then becomes turbulent. The change from laminar 
to turbulent flow depends on the conditions of gas flow. 
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Fig. 6.2. Traverse of velocity distribution in a turbulent flow. 


The law of velocity change across the thickness of a boundary layer is 
different for laminar and turbulent flows. In the turbulent boundary layer, 
due to the intensive mixing of gas layers, the velocity increases considerably 
faster than in the laminar layer. 

The conditions of gas flow along the rocket motor nozzle are such that 
the boundary layer is always turbulent, with a rapid increase in velocity 
to the full velocity of the stream (Fig. 6.2). Therefore, in examining the 
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Fig. 6.3. Schematic of a one-dimensional flow in the nozzle of a rocket motor. 


flow of gas along the nozzle of a liquid rocket motor, the boundary layer is 
usually not taken into account, considering that, at any point of a given 
section, the velocity is the same and is equal to the velocity in the center of 
the flow where the influence of the walls is not felt. 

The nozzle of a rocket motor is a duct with changing cross section, 
in which radial gas flow must take place toward the axis of the duct in 
the convergent part and away from the axis in the divergent part. How- 
ever, the radial velocities of gas in the nozzle may be neglected as a first 
approximation. 
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Therefore, at each section of the motor, the velocity is considered to 
be constant and equal to the velocity at the axis of the duct (Fig. 6.3). 
Such fluid flow is called one-dimensional. 

If the radial velocity is not considered, then the calculated (theoretical) 
velocity of flow will be greater than actual. Consequently, we can consider 
that the radial flow results in some loss in exit velocity and specific impulse 
of the motor. This loss is included in the total losses which accompany the 
flow of gas along the nozzle. 


4. EQUATION OF DISCHARGE 


Let us examine two sections, perpendicular to the direction of velocity 
of one-dimensional flow (Fig. 6.4), and calculate the mass of gas passing 
through both sections during the time At. The mass discharge is determined, 
in this case, by the volumetric discharge SwAt, multiplied by the density p. 
Therefore, the gas mass passing through the first section of one-dimensional 
flow will be p:w,S,At, and through the second, p2w2S2At. 
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Fra. 6.4. Derivation of consumption equation. 


However, in a steady state flow the gas parameters cannot change at 
any point between the first and the second section during the time At. 
Therefore, there will be no concentration or diminution in the gas mass in 
the volume bounded by the sections. The passage of gas through section 
1-1 must therefore equal the discharge through section 2-2, from which it 
follows that for one-dimensional gas flow in the steady state 


pwS = constant. (6.3) 


The derived equation is called the equation of discharge. It represents 
the expression of the law of conservation of mass for the case of gas flow. 

For an incompressible fluid p = constant. The equation then takes the 
form 


wS =_constant. (6.4) 
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At low velocities, gas may be considered an incompressible fluid. From 
Eq. (6.4) it follows that, for this case, the velocity will be inversely propor- 
tional to the area of the cross section of the stream. For a compressible gas, 
due to change in density p, the picture changes not only quantitatively 
but also qualitatively. As we shall see later on, at supersonic velocities in a 
divergent duct the gas velocity does not diminish, but increases. 


5. EQUATION OF ENERGY 


Let us examine the energy relationships characteristic of gas flow. 
These relationships result from the law of conservation of energy. Let us 
consider that a heat exchange between the flow and the walls (surrounding 
medium) is absent, i.e., we shall examine an adiabatic gas flow. In this 
case, the energy, E, of some gas mass will remain constant for any section 
of the stream. 

Let us calculate the value of E for 1 kg of gas at some section of the 
stream. 

Gas energy is composed of kinetic energy, potential energy due to pres- 
sure and weight, and the internal (heat) energy of the gas. First of all, at 
velocity w, the gas has a kinetic energy equal to 


mu? /2 


or, for 1 kg of gas, w?/2g, and, in heat units, Aw?/2g. 

The potential energy due to pressure of 1 kg of gas in heat units is equal 
to Apv. The potential energy of weight, mgz, for 1 kg of substance is equal 
to the height, z, of the center of gravity of the gas mass with respect to 
some datum line. For low density gas flows, the potential energy due to 
weight is negligible compared to the potential energy due to pressure. 

Finally, the internal gas energy 


6 fT es ie 


Therefore, if the potential energy due to weight is neglected, the total 
energy of 1 kg of gas 


E = Aw®/2g + U + Ap. 


The sum of the internal energy, U, and the potential energy due to 
pressure, A pv, is called the heat content, H (see Chapter V). Therefore 


E = Aw*/2g + H. (6.5) 
For adiabatic flow, E = constant. Consequently 


H + Aw*/2g = constant. (6.6) 
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or, for two arbitrary sections, 0—0 and 1-1 of the gas stream, 
Ho + Aw,?/2g = H, + Aw,?/29. (6.7) 


The derived equation represents an expression of the Jaw of conserva- 
tion of energy, which states that energy does not disappear or appear, but 
is transformed from one form into another. Equation (6.7) will be referred 
to from now on simply as the energy equation. 

Equation (6.7) is often used in determining the velocity of gas flow. 


w, = V we + 2g(Ho — Hy)/A. (6.8) 


If the flow of gas from a vessel of large dimensions is considered, where 
velocity wo is low, then 


Wi = /2g(Ho — H,)/A. (6.9) 


Let us present the energy equation in a different form, which will be 
useful to us in the future. In order to do this, we will take advantage of 
relationships presented in Chapter V, and will express the heat content of 
the gas in terms of the parameters of state 


H = kRT/(k — 1) = kp/gp(k — 1). 
The energy equation can now be rewritten in the form 
kp/p(k — 1) + w?/2 = constant. (6.10) 


In deriving the energy equation for incompressible fluid the heat energy 
should not be taken into account, since it does not change during the motion 
of such a fluid. On the other hand, due to the large specific weight of liquids, 
it is imperative to consider the change of the potential energy due to weight. 
The equation of conservation of energy for incompressible fluid is written 
in the following form: 


p/y + w*/2g + z = constant. (6.11) 


The equation of energy for an incompressible liquid Eq. (6.11), is called 
Bernoulli’s equation. 


B. Gas Flow in a Supersonic Nozzle 


1. VELOCITY OF SOUND IN GASES 


A very important gas characteristic is the speed of sound propagation. 

By the term speed of sound, we mean the velocity of propagation of 
longitudinal vibrations in the medium. This includes not only vibrations 
which are audible to the human ear as sound, but also gas vibrations whose 
frequencies are beyond the threshold of hearing. 
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Let a motionless gas mass at pressure p, density p, and temperature T, 
be confined in a cylindrical tube (Fig. 6.5). Further, let an impulse—for 
instance, a short shock—be imparted to the gas with the aid of a movable 
piston at the left end of the tube. The gas in the vicinity of the piston will 
compress and then, in expanding, will set in motion particles of gas located 
to the right. A wave will start to travel along the tube from left to right. 

At a certain time t,, the wave will reach section 1-1. After At sec (at 
time é), it will be displaced to section 2-2. 
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Fic. 6.5. Derivation of the speed of sound. 


Beyond section 2-2, the gas pressure will be the same as before the 
imparting of the impulse to the gas. To the left of section 2-2, the pressure 
will be greater than in the undisturbed gas. Let us designate it by pı, (pı > p). 

The velocity of wave propagation obviously will be c = Az/At, where 
Ax is the distance between sections 1-1 and 2-2. 

The increase in the gas mass in the tube volume corresponding to the 
region Az is due to the inflow into this volume, from the left, of a certain 
gas mass with some velocity w. The magnitude of this mass is 


AM = prSwat ) 


where S is the area of the tube cross section, and pı is the gas density to the 
left of the wave front. 

On the other hand, the increase in the gas mass can be expressed in 
terms of the change in density within the volume SAz 


AM = (p — p)SAz. 
Equating the right members of the last two equations we find 


w = (pı — p)e/pr. (6.12) 
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In order to eliminate the unknown velocity w, we will take advantage 
of the momentum theorem. 

The gas mass in the volume SAz will be pSAz. This mass during the 
time At begins to move with the velocity w. The change in the momentum 
must be equal to the impulse 

pS(w ~ O)Ar = (pı — p)SAt, 
from which 
Pi — P = pwe. 


Substituting into this equation expression (6.12) for w we get 


c = Vpi(~m — p)/p(p — p). (6.13) 


For weak gas disturbances such as sound vibrations, p, does not vary 
appreciably from p, and pı from p 


pı = p + Ap; Pi = p + Ap (6.14) 


c= VAp/Ap. (6.15) 


The quantity Ap/Ap depends on the process of the gas compression. 
Newton, who first derived this expression in 1687, supposed that the gas 
temperature remained unchanged during the passage of a wave. In this case 


and 


p/p = constant 
from which we get 
Ap/Ap = p/p 
and 
c= V p/p. 

This formula gives a value for the speed of sound in air almost 15% less 
than that obtained experimentally. In his time, Newton explained this 
discrepancy by the presence in the atmosphere of suspended solid particles 
and water vapor. 

Considerably later, in 1810, Laplace showed that the process of gas 
compression during the passage of a wave should be regarded not as isother- 
mal but as adiabatic, inasmuch as during rapid compression and expansion 
of the gas the heat transfer within a gas does not have time to take place. 

In the case of an adiabatic process in a wave 


p = p* constant 
Ap = kp*Ap constant 
from which 
Ap/Ap = k(p/p) 
and 


c= Vk(p/p), (6.16) 
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or, in accordance with the equation of state, p = pgRT, 
c = VkgRT. (6.17) 


The last two formulas yield a higher value for the speed of sound 
and are in good agreement with experiments. 

Therefore, the speed of sound in gas depends not on absolute values 
of pressure and density but on their relationship, i.e., on the temperature. 

For air, k = 1.4, R = 53.3, and the expression for the velocity of sound 
takes the form 


c = 49.0 yT. 


The velocity of sound in air at 32°F equals 1080 ft/sec. For combustion 
products in the chamber of a rocket motor when T = 5410°R, k = 1.2, 
R = 62.6, and speed of sound c = 3610 ft/sec. 

The speed of sound has a definite physical significance, representing the 
velocity of propagation of weak disturbances in a gas. The velocity of 
sound also has a definite energetic significance. In order to clarify it, let 
us rewrite the expressions for the square of the sound velocity and for the 
heat content of a gas. 

c = kgRT 
H = kART/(k — 1). 


Eliminating from these two equations the term T, we will find 
e = HA/g(k — 1). 


The obtained expression means that the square of the speed of sound is 
a measure of the heat content of a gas. 

The concept of the speed of sound has a tremendous significance in 
aerodynamics and gas dynamics. The flow of gas past bodies, through tubes 
and fittings, and, in general, the behavior of any gas motion, depends upon 
the ratio of gas velocity to the velocity of sound in a gas. Depending upon 
the value of this ratio, it is customary to speak of subsonic and supersonic 
regimes of flow and flight velocities. The ratio of flow velocity to the 
velocity of sound is cutomarily represented in all aerodynamic and gas 
dynamic calculations by the letter M, and is called the Mach number 


M = w/c. (6.18) 


Initially (in 1868), the relationship of gas velocity (or a body moving 
through the gas) to the velocity of sound was introduced into scientific 
circles by a Russian scientist—ballistician N. V. Mayevsky. Later, it was 
also used by an Austrian physicist, Mach, and is widely known in tech- 
nology as the Mach number. 
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Let us return to expression (6.13). 

If the pressure pı and density pı differ considerably from the pressure p 
and density p in the undisturbed gas, then the disturbing wave is called a 
strong, or shock wave, as differentiated from the weak accoustical (sound) 
wave. The velocity of propagation of the shock wave is 


c = VIl + Ap)/p](Ap/Ap). (6.19) 


By comparing expressions (6.15) and (6.19) it follows that the velocity 
of propagation of the shock wave is always greater than the velocity of 
sound. 

If a strong disturbance is imparted to the gas, i.e., if a large difference 
in pressures, pı — P, 1s created, then the formed wave, during its propaga- 
tion, will partially dissipate energy imparted to the gas. The strength of 
the wave, measured by the difference of pressures pı — p, will be decreasing. 
Its velocity will decrease correspondingly and the shock wave, after some 
time, will turn into a weak wave propagating with the velocity of sound. 


2. TuE Maximum EXHAUST VELOCITY 


Let us imagine a vessel, for instance, the combustion chamber of a 
rocket motor, which contains an immobile gas (wo = 0), with unchanging 
parameters, po, po, To. Let the gas flow out of this vessel into a region where 
the parameters of this gas will be p, p, and T. 

It follows from the energy equation that the velocity of gas will be great- 
est in those sections of the stream in which the heat content is least. A maxi- 
mum velocity will be obtained if all of the heat content is converted into 
the kinetic energy of the stream of the exhaust gas. At the same time, 
absolute temperature of the gas must become equal to zero. From the 
equation of the adiabatic process which is assumed to take place during 
gas flow from the vessel, 


T = To(p/po)*-?" (6.20) 


it follows that in order to get gas temperature equal to zero it is necessary 

for the pressure, p, in a gas stream to equal zero also. Therefore, the maxi- 

mum velocity of a gas can be obtained by exhausting into a vacuum. 
According to expression (6.9) the value of the maximum velocity is 


Wmax = V (29/A)Ho, 


where Ho is the heat content of the immobile gas in the vessel. 
The term Wmax may be expressed in terms of parameters of state and 
velocity of sound in a motionless gas 


Wmax = V 2kpo/(k — 1)po = V 2kgRT/(k — 1) = a V2/(k — 1). (6.21) 
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For air at room temperature, Wmax = 2,460 ft/sec; for combustion 
products of rocket motor fuel (To = 5410°K, R = 62.6, k = 1.2), Wmax 
~ 11,500 ft/sec. 

The maximum exhaust velocity, as follows from formula (6.21), depends 
only on temperature To, and is independent of pressure. From the energetic 
point of view this is completely clear. The maximum velocity will result 
when the initial heat content of the gas is completely converted into kinetic 
energy, and the amount of heat content is determined only by the initial 
temperature of the gas. At maximum velocity there is a complete conversion 
of the thermal random motion of molecules into the directed motion of the 
flow. 

It would seem, at first glance, that with increased pressure the flow 
velocity Wmax must increase, inasmuch as, conventionally speaking, the 
force pushing the gas out of the vessel increases. However, with the increase 
in pressure, the density po increases in the same proportion (when To = con- 
stant) and, consequently, the mass contained in a unit volume also in- 
creases. It is understandable that the increased pressure imparts to the 
mass, which has been increased in the same proportion, the same velocity 


W max: 


3. DEPENDENCE OF GAS PARAMETERS ON LOCAL FLOW VELOCITY 


Let us examine the manner in which the parameters of a moving gas 
will change with respect to flow velocity. Let us write the energy equation 
for two states of flow 


kgRTo/(k — 1) = 4w? + kgRT/(k — 1). 


Here, w and T are related to some arbitrary section, and To represents the 
temperature of gas when w = 0. For temperature T we get the following 
expression: 


T=T p. 
B (1 — QkgRTo/(k — 5) 


or 





T = ro(1 Sage ) (6.22) 


2 
Wmax 


Consequently, the higher the flow velocity of a gas the lower its tempera- 
ture. 
In an adiabatic gas flow 


p/po = (T/T); p/po = (T/T), 
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Consequently, 
ap? ki (k—1) 
p = Po (1 — =) (6.23) 
2 \1/(k—1) 


Therefore, both the pressure and the density of the gas are reduced 
with the increase in flow velocity. In reaching maximum velocity, the 
pressure and density of the gas are reduced to zero. 


4. Tor RELATIONSHIP BETWEEN THE LOCAL SPEED OF SOUND 
AND FLow VELOCITY. CRITICAL VELOCITY 


As we already know, the velocity of sound in a gas is determined only 
by the temperature. However, the gas temperature may be different at 
various points within the stream. Consequently, the velocity of sound will 





Wimax 


Fia. 6.6, Dependence of the speed of sound in a gas on flow velocity. 


also be different. Therefore, in considering a moving stream, we must con- 
sider not only the local velocity, temperature, pressure, and density, but 
also the local velocity of sound. 

The energy equation, after introducing the concept of maximum veloc- 
ity, may be written in the following manner: 


tw? + kgRT/(k — 1) = Whar (6.25) 
From this, taking into consideration that kgRT = c?, we get 
c= V (Wax — w)(k — 1)/2. (6.26) 


The relationship between the speed of sound and gas velocity is shown 
in Fig. 6.6. When w = Wmax, the speed of sound is reduced to zero, inasmuch 
as, under these conditions, T = 0. 
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The flow velocity may be either greater or less then the local speed of 
sound. The flow velocity equal so the local speed of sound is called the 
critical velocity. It may be found by means of expression (6.26) by letting 
C= W = Wer. 


Wer" = Wraax (ke = 1)/(k + 1): 
But since 
Winax = ZkgRTo/(k — 1), 
then 
Wer = BZkgRTo/(k + 1) = 2c/(k + 1). (6.27) 


Therefore, the value of w., depends only on the temperature of gas in 
the vessel from which flow takes place. 
Let us establish the local gas parameters at critical velocity. According 


to expression (6.22) 
2 
fies = T (1 a = ) 
Wmax 





or, since Wert = Whax(k — 1)/(k + 1) 
Ter = 2To/(k + 1). (6.28) 
Considering expressions (6.23) and (6.24) we also get 
Per/Po = [2/(k + 1) 4%; (6.29) 
Por/ po = [2/(k + 1)]VE-”. (6.30) 


In general, for air (k = 1.4) 
Per = 0.528 Po; 


Te = 0.833 To; 
po = 0.634 po; 


for combustion products of a rocket motor (k = 1.2) 


Der = 0.565 po; 
Toe = 0.909 To; 
Per = 0.621 po- 


Therefore, in order to get supersonic airflow exhausting to the atmos- 
phere, it is necessary to have a chamber pressure approximately twice that 
of the atmospheric pressure. 


5. CONFIGURATION OF A SUPERSONIC NOZZLE 


Up to now we have discussed dependence of the gas fiow parameters 
on velocity, but we have not examined the problem of reaching this velocity. 
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It is well known that the increase in velocity of an incompressible fluid or 
gas may be obtained by the constriction of the duct. However, experience 
shows that in a constricted duct, velocities greater than critical may not 
be obtained. Supersonic flow velocities of a gas can be obtained by means 
of the Laval nozzle, which consists of a duct whose cross section first 
decreases and then increases (Fig. 6.7). 

It has been found that, if gas flows through a Laval nozzle in such a way 
that it reaches critical speed at the narrowest section, then beyond the 
throat the flow velocity, w, begins to increase. In this manner, it has been 


Fic. 6.7. Laval nozzle configuration. 


determined by experiment that for supersonic flow, the flow laws are 
directly opposite to the flow laws of a subsonic stream. In other words, a 
stream flowing with subsonic velocity is accelerated in a convergent duct 
and decelerated in a divergent section. Supersonic flow, on the contrary, 
is decelerated in a convergent duct and is accelerated in a divergent duct. 

Let us examine this question in more detail. Let us write the discharge 
equation, Eq. (6.3), in the form 


S = constant/pw. (6.31) 


For an incompressible fluid (p = constant), the velocity is inversely 
proportional to the area of the duct section. For a compressible fluid at 
normal velocities, the density, p, as follows from formula (6.25), decreases 
with increase in velocity, but decreases so unnoticeably that the character 
of the gas flow remains qualitatively the same as for the incompressible 
fluid, and the area of the flow section decreases with increase in velocity. 
However, at high, supersonic velocities the value of p decreases with in- 
crease in velocity faster than the increase of w. Consequently, the flow 
cross section S, in this case, must increase. 

Let us examine the dependence of pw on w. According to expression 
(6.24) 
poll — w?/wha) 0. 


> 
I 


Consequently 
pw = pow ~ w?/whrax)/O, 


The term pw may be called the mass specific discharge, which means 
the discharge of a gas mass in a unit of time through a unit area of cross 
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section. When w = 0, the specific discharge is reduced to zero. When 
W = Wmax, it also reduces to zero since, at this velocity, the density p 
becomes zero. In the interval between these two limiting values of velocity 
there obviously exists a maximum for the function. 

Fig. 6.8 shows the dependence of the specific discharge on velocity. 
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Fra. 6.8. Relationship between the quantity pw and the flow velocity w. 


Let us see at what value of velocity the maximum specific discharge 
takes place. Differentiating the expression for pw with respect to w, and 
equating the derivative to zero, we get 
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Consequently, either w = Wmax, or the expression within brackets in 
Eq. (6.32) equals zero. 


In the first case, pw = 0, and we obtain a minimum shown on the curve 
of Fig. 6.8. In the second case, when 


w = Winax(k — 1)/(k + 1) (6.33) 


the specific discharge pw has its maximum value. 
However, the velocity determined by the relationship (6.33) is none 
other than the critical velocity equal to the local velocity of sound. The 
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maximum point on the curve of Fig. 6.8 separates subsonic flow velocities 
from supersonic flow velocities. 

In examining the curve pw, it is not difficult to determine how the trans- 
verse section of the duct must vary in order to obtain supersonic velocities. 
The cross sectional area, according to Eq. (6.31), must first decrease and 
then increase. At the point where the term pw has its maximum value, the 
area of section S must be minimum. The local velocity of sound and the 
critical velocity are reached at this point. Therefore, the minimum section 
of the nozzle is called the critical section. In this manner we see that attain- 





Fig. 6.9. Variation in flow parameters along a supersonic nozzle for different values 
of various gas parameters at the entrance. 


ment of supersonic velocities is possible in the Laval nozzle. However, 
this is not always possible. If the difference in pressures at the entrance 
and exit of the nozzle is insufficient for creation of critical velocity at the 
narrow section, it will be impossible to obtain supersonic velocity. 

Fig. 6.9 shows possible cases in the operation of the nozzle. The solid 
curves correspond to the basic operational case. The relationship of pressure 
at the entrance to the exit pressure is less than a critical ratio 


p/po < [2/(k + 1)]He-», 


In the meantime, the gas pressure drops continuously and the velocity 
increases. 
In cases in which the pressure ratio 


p/po > [2/(k + 1), 


i.e., when the exit nozzle pressure becomes too great or too small at the 
entrance, the critical velocity at the narrow section of the nozzle will not 
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be attained, even though the velocity there will be maximum. In the 
diverging portion of the nozzle the velocity will be decreasing and the 
pressure increasing (dashed curves in Fig. 6.9). 


C. The Performance of a Rocket Motor Nozzle 


1. THE AREA OF THE CRITICAL NOZZLE SECTION 


The geometrical dimensions of a supersonic nozzle of a rocket motor 
are determined by the areas of the critical and exit sections. 

First, let us find the required area of the critical section. For this, let 
us express the mass rate of fuel (gas) discharge G, in terms of the gas 
parameters at the nozzle critical section, 


G = 9PerWerSer- 
But, according to Eqs. (6.30) and (6.27) 
per = pl2/(k+1)V/FY and wet = V2kgRTo/(k + 1). 


Therefore 
G = gpSec2/(k + DIYO- V 2kgkTo/(k + 1). (6.34) 
Since 
po = po/ghTy 
then 


G = (poSer/WVRT»)[2/(k + 1)]}/4- V/2kg/(k + 1). (6.35) 


Equation (6.35) relates the gas discharge, G, to the area of the critical 
section of the nozzle, Ša, and to the gas parameters po and To in the combus- 
tion chamber. The last, as we already know, are determined from the com- 
bustion calculations in the liquid propelled rocket. 

It may be seen, incidentally, from expression (6.35) that the increase 
in pressure, po, leads to the proportional increase in discharge G, while the 
flow velocity We remains unchanged. In this manner, the volumetric gas 
discharge Werer, through the critical section, and consequently through the 
entire nozzle, at constant chamber temperature remains constant. This 
means that at an increased pressure, the discharge G increases only because 
of the increased density of the gas p in the stream. At constant temperature, 
increase in density leads to the proportional increase in pressure at all 
points within the stream. It should be noted that the increase in tempera- 
ture To at constant gas discharge rate also leads to the increase in pressure 
within the chamber and the nozzle. 

It follows from expression (6.35) that at mcreased chamber pressure 
the value of the required discharge G may be obtained with smaller dimen- 
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sions of the critical section Se, resulting in the reduction of all other nozzle 
dimensions. 

The composition of the combustion products influences the required 
area of the nozzle section by the values of the adiabatic exponent k and 
the gas constant R. The greater the amount of substances with a small 
number of atoms in the molecule and with low molecular weight, which are 
characterized by large values of k, in the products of combustion, the lower 
the gas discharge through a nozzle will be with a given area of section Ser. 

In many types of rocket motors (specifically, aircraft), it is necessary 
to regulate the thrust. This is most simply accomplished by varying fuel 
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Fic. 6.10. External view of a four chamber motor. 


discharge G. However, such regulation of thrust, especially if it is performed 
within wide limits of discharge variation, leads to a considerable change in 
chamber pressure. Reduction in the chamber pressure affects the combus- 
tion process adversely and, specifically, can give rise to pulsating burning. 
Reduction in chamber pressure also leads to the impairment of nozzle per- 
formance (see page 187). Therefore, if it is necessary to regulate thrust 
within wide limits, the motor nozzle must be designed with a variable criti- 
cal section in order to maintain constant chamber pressure. 
Mechanization of the nozzle (variation of the critical section during 
operation of the motor) represents an extremely difficult problem. There- 
fore, in practice, this problem is solved by designing motors with several 
combustion chambers. When it becomes necessary to reduce the thrust, one 
or several chambers are cut off and the remaining chambers continue to 
operate with normal or near normal combustion chamber pressures. Motors 
with two, and even four, chambers are known to be in existence (Fig. 6.10). 


2. AREA OF THE NozzLE Exit SECTION 


It has been shown above that the attainment of the maximum flow 
velocity, Wmax, is related to the complete expansion of the gas to p = 0, 
and, consequently, requires an infinitely large area of the exit section of the 
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nozzle S,. This is a very important condition, since, as we now see, for the 
attainment of maximum velocity the mere presence of a vacuum in the 
environment is insufficient. This is understandable, since for the attainment 
of maximum velocity it 1s necessary to cool the gas to T = 0 within the 
stream, i.e., to obtain a vacuum, not in the environment, but at the exit 
plane of the motor nozzle. 

The effort to increase flow velocity may lead to the unreasonable 
increase in the dimensions and weight of the nozzle. For this reason, the 
maximum velocity remains unattainable in practice and the value of flow 
velocity is limited by the consideration of the reasonable dimensions of 
motor nozzles. 

Let us find the connection between the dimensions of the exit section 
of the nozzle, Se; flow velocity, w.; and the gas exit pressure, pe. In order to 
do this, we will express the flow velocity, we, making use of Eq. (6.9), in 
terms of heat content of gas in the combustion chamber and upon exit from 
the motor. 


we = V (29/A)(Ho — H.). 


Since 
HS Cyl 5 
Te = To(pe/po)*-?"*; 
and 
Cp = kKAR/(k — 1); 


then the exit velocity 


We = V[2kgRTo/(k — DJL — (pe/po) 4]. (6.36) 


Further, let us express the time rate of discharge of the fuel, G, in terms 
of gas parameters at the exit plane of the nozzle 


G = gpeWeSe. 
Taking into account that 
Pe = po(De/Po)*/* 


and also using expression (6.36), we find 


G = gpoSe(pe/po)"* W{2kgRTo/(k — VILL — (pe/po)*-Y*]. (6.37) 


Equating the right side of Eq. (6.37) with the right side of Eq. (6.34), 
we get 


Se VIG = 1)/2/12/& + Dee? 


(6.38) 





Sa V (De/ p) E — (pe/ po) EDT 


As is seen from formula (6.38), the pressure ratio po/p, for a given gas 
with a constant adiabatic exponent k depends only on the ratio S,/Ser. 
The last we shall call the nozzle expansion ratio. 
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Relationship (6.38) enables determination of gas parameters at any 
section of the nozzle S,. For this S, must be substituted for Se, and p+ for 
pe in Eq. (6.38). Knowing the values of S+ and Ser, Pz/po can be calculated 
for a given section and then, with the aid of formulas for the adiabatic 
process, all the other flow parameters may be found for a section with an 
area Sz. 
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Fig. 6.11. Variation in gas flow parameters along the rocket motor operating on nitric 
acid and kerosene (k = 1.2). 


Equation (6.38) may be used for the subsonic as well as for the super- 
sonic sections of the nozzle. In this manner, it 1s possible to determine the 
changes in parameters of gas flow along the entire length of the nozzle. 

A typical graph of parameter variations of gas flow along the gas duct 
of a rocket motor is shown in Fig. 6.11. The variation in gas parameters 
along the length of the combustion chamber is shown by dashed lines, 
inasmuch as these variables are not sufficiently clear. From this graph it 
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is seen that the more rapid variation in parameters (with the exception of 
the Mach number) takes place in the region of the critical section of the 
nozzle. In a divergent portion of the nozzle, parameters vary at an ever- 
decreasing rate. 


3. DIVERGENCE OF A NOZZLE. 
THERMAL EFFICIENCY OF A Rocket MOTOR 


The expressions obtained above for the gas flow parameters permit an 
evaluation of the degree of recovery of the gas heat content in a rocket 
motor. In moving along the nozzle, the velocity, and consequently the 
kinetic energy, of the stream increases. The increase in kinetic energy is 
accompanied by a corresponding decrease in the heat content of the gas; 
the higher the speed attained by the gas at nozzle exit, the greater is the 
part of its heat content that is converted into useful kinetic energy. 

The complete conversion of the entire heat energy into kinetic energy, 
i.e., the attainment of maximum velocity Wmax, at the exit plane of the 
nozzle is impossible, as we have already seen. This would require building 
a nozzle with infinitely large dimensious of the exit plane. Therefore, in 
actual motors, it is necessary to limit the conversion to only a part of the 
total heat energy. The greater the expansion ratio, i.e., the ratio S./Ser, 
the lower the temperature and pressure of gas at the nozzle exit are, and the 
smaller the portion of available energy that is lost in the gas stream leaving 
the nozzle. 

The degree of conversion of heat energy into kinetic energy is defined 
by the thermal efficiency of the motor, 7. It is natural to define thermal 
efficiency as the ratio of kinetic energy of the gas stream to heat content 
of the gas at nozzle entrance 


ne = A(w?/29)/Ho. 
Since, according to the energy equation, 


A(w?2/2g) = Ho — H. 
then 
n= (Hy — H)/Ho = 1 — H./Ho. 


Ordinarily, thermal efficiency is expressed as a ratio of pressures at the 
entrance and exit of the nozzle po/pe, or its reciprocal, p./po 


ne = 1 — (pe/po)*-?"*, (6.39) 


The relationship of thermal efficiency of the motor to pressure drop and 
adiabatic exponent is shown in Fig. 6.12. 

We must also note that thermal efficiency may be expressed as the 
ratio of the square of exhaust velocity to the square of maximum velocity. 
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The pressure ratio, Po/Pe, is directly connected with the term Se/Se in 
Eq. (6.38). The interrelationship of these terms is influenced by the value 
of the adiabatic exponent of the combustion products. Fig. 6.13 shows the 
variation of po/pe with S./Se for two values of k. 

Relationship po/p- is one of the basic parameters in designing the nozzle. 

The calculation of thermal efficiency by formula (6.39) is possible only 
when the pressure at the nozzle exit plane is equal to atmospheric pressure. 
Otherwise, other, more complex relationships must be applied, which take 
into account the effect on the operation of the motor of excessive static 
pressure or rarefication at the exit plane of the nozzle. 
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Frc. 6.12. Relationship of a rocket Frc. 6.13. Relationship of nozzle pres- 
motor thermal efficiency to the pressure sure ratio po/pe to nozzle expansion ratio 
ratio of the nozzle po/p,. for the two values for two values of k. 
of k. 


The curves of Figs. 6.12 and 6.13 show that the quantity k has a pro- 
nounced effect on the dimensions of the nozzle and thermal efficiency. 
Thus, the increase in the adiabatic exponent of the combustion products 
results in a decrease of nozzle dimensions and increase in thermal efficiency. 
Therefore, an increase in k reflects advantageously on the dimensions of 
the postcritical part of the nozzle and on the degree of recovery of heat 
energy in the motor. 


4. Tue EFFECT oF Gas RECOMBINATION AND COMPLETION OF COMBUSTION 
oN Gas FLOW ALONG THE NOZZLE OF A ROCKET MOTOR 


Until now, we have discussed the flow of gas of a constant composition 
along the nozzle. Partial recombination of dissociated gases and comple- 
tion of combustion of fuel which was not burned in the chamber take place 
within the products of combustion during their movement along the nozzle, 
and is accompanied by heat release. Besides, the composition of gases 
changes in the nozzle. This directly affects the value of the gas constant 
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R and adiabatic exponent k. The change in the composition of gas, however, 
does not have a marked effect on the process of expansion. A most substan- 
tial effect is the release of an additional quantity of heat energy in the 
nozzle. 

An energy equation derived earlier must be, in view of what has been 
said, correspondingly revised. We must introduce, into the over-all reserve 
of energy, chemical energy as well. In order to do this we must substitute 
for the heat content of the combustion products H in Eq. (6.5), the total 
heat content, E, and assume as before 


E = E, + Aw?/2g = constant 
or, for the two states corresponding to sections 0—0 and a-a 
Eo + Aw?/2g = Ee + Awe/2g. 
Equation (6.8) for calculating velocity will have the form 
We = Vwè + 2g(Ey — E.)/A. 


Specifically, for exhaust velocity from a nozzle of a rocket motor, we 
will get an expression 


We = V 2g(Eo a E.)/A. 


The total heat content of the combustion products in the chamber, 
Eo = E,, is determined from combustion process calculations. In the ab- 
sence of heat losses 


E, = Ey. 


The total heat content of the combustion products at nozzle exit plane, 
E., is determined in the same manner as the total heat content in the 
chamber, according to temperature and composition. 

It is quite obvious that the exhaust velocity, taking into account 
recombination and completion of combustion, will be greater by comparison 
with the value of we, calculated on the assumption that the gas mixture 
flowing along a nozzle is unchanged. Under these conditions the thermal 
efficiency of a motor may be found by the formula 


Ne = (Ko =~ E.)/ Ke, 
where Kg is the gravimetric heating value of the fuel. 


5. Rocker Motror NOZZLE CONFIGURATIONS 


Liquid rocket motors known at the present time are characterized by 
the following basic parameters. 
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Combustion chamber pressure is within the limits of 25-50 Ib/in.. 
Increase in the chamber pressure, all other conditions being equal, Jeads to 
the increase of the specific impulse and to the decrease in the dimensions 
of the nozzle critical section, and consequently all other dimensions, and 
to some decrease in the dimensions of the combustion chamber. Therefore, 
in the future, we may expect pressure increase In the combustion chamber 
of rocket motors. 

The pressure ratio po/Pe, in contemporary motors, equals 20-40. There- 
fore, the expansion ratio Se/Ser equals 4-7. A further increase in po/pe may 
be expected, both as a result of increased chamber pressure and of the reduc- 
tion of pressure at the exit plane of the nozzle, which is desirable for long 
range motors and certain types of aircraft engines. 

The elementary theory of the supersonic nozzle, examined above, shows 
that the velocity of the exhaust gas from the nozzle and the pressure at 
the exit plane do not depend on its configuration. It is necessary only that 
the nozzle first converge and then diverge. The change in nozzle profile, i.e., 
variation of section area with respect to length, reflects only on the change 
in velocity along the length. Pressure and velocity at exit are determined 
only by the ratio of the exit area Se to the area of the critical section Ser. 

The geometric configuration of the nozzle must be chosen so as not to 
result in large velocity losses in the nozzle. At the same time, its surface 
area must be as small as possible. With an increase in surface area there 
is an Increase both in the weight of the nozzle and the amount of heat trans- 
ferred to the cooling liquid. 

Losses in the nozzle result from friction, shocks, and vortices in the gas 
stream. Because of this, part of the kinetic energy of the directed motion is 
converted into unusable, for generation of thrust, heat energy. In this 
connection, let us examine separate sections of the nozzle profile. 

In order to decrease friction losses, the nozzle entrance must be made 
as short as possible and, at the same time, smooth, to eliminate shock and 
turbulence. 

In order to decrease surface area, the postcritical section should be 
made as short as possible, with a correspondingly increased expansion 
angle. 

However, it may so happen that, with a large expansion angle, the gas 
stream flowing with high velocity will not be able to expand and fill the 
cross section of the nozzle. Flow separation from the walls will take place, 
which will result in sharply increased losses due to formation of vortices. 
Therefore, the expansion angle of a nozzle must be limited. 

The allowable nozzle expansion angle, in liquid propelled rockets, is 
quite large. At the beginning of the nozzle divergent section it may be as 
high as 35°-40° (Fig. 6.14). Intensive reactions of combustion and recom- 
bination, which are accompanied by heat release, still take place here. For 
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this reason, the radial expansion of the stream is promoted and flow separa- 
tion from the wall will not take place even at these large expansion angles. 

At the exit of the nozzle the expansion angle must be reduced. Here, 
the reactions of recombination and combustion are considerably less 
intense, the supply of heat is reduced, and the conditions for expansion 
are less favorable. Besides, it is desirable to reduce the expansion angle at 
nozzle exit in order to reduce the radial component of stream velocity (Fig. 
6.14), which does not contribute to thrust. It is desirable that the entire 
flow have a direction parallel to the axis of the rocket. 


t 


= 
Fic. 6.14. Loss of kinetic energy due to Fre. 6.15. A reasonable nozzle configu- 
radial component of gas flow. ration with large expansion. 


For technological considerations it is more convenient to make the diver- 
gent portion of the nozzle in the shape of a cone. However, in this case, in 
order to reduce losses, the nozzle expansion angle selected must be smaller 
(20°-25°), and the resulting nozzle will be long with a large surface area. 

A more desirable nozzle is of the configuration shown in Fig. 6.15. 
Transition from a large expansion angle at the critical section to a small 
angle at the exit is accomplished by a circular arc. The nozzle length and 
losses due to radial flow are reduced to the absolute minimum. Such nozzle 
configuration, however, does not exclude the possibility of the occurrence 
of local losses. In order to eliminate losses completely, it is necessary for 
a nozzle profile to correspond to the trajectory of gas particles flowing 
near the walls. This trajectory is completely defined if the expansion ratio 
of the nozzle and distribution of velocities at the critical section and at the 
nozzle exit plane are given. Nozzles designed in this manner are called 
dynamically profiled nozzles. 

Methods of profiling, known from aerodynamics, yield very long nozzles 
and, consequently, large wall surface areas. Because of this, dynamically 
profiled nozzles are not used in rocket motors. 


D. Properties of Supersonic Nozzle. Characteristics of 
Rocket Motors 


1. PROPERTIES OF THE SUPERSONIC NOZZLE AND 
REGIMES OF ITS OPERATION 


A supersonic nozzle has a series of interesting properties which reflect 
on the operation of a rocket motor. 
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The first of these properties is the independence of the parameters of 
combustion products in the chamber and nozzle of the motor from external 
atmospheric conditions. This is explained by the fact that at supersonic 
exhaust, atmospheric disturbances cannot penetrate into the nozzle against 
the flow which moves with a velocity greater than the speed of sound. 
Therefore, with the rise of the rocket to an altitude, the pressure at the 
nozzle exit plane (during constant chamber pressure) remains unchanged. 
On the other hand, if the motor operates with a variable fuel consumption, 
the pressure at the nozzle exit plane will vary following a change in chamber 
pressure, independently of atmospheric conditions. 

The velocity at nozzle exit, with a given expansion ratio, will remain 
constant under any flow conditions if only the temperature and the gas 
composition in the combustion chamber remain constant. The exhaust 
velocity is also not affected by chamber pressure. 

Let us examine the so-called design and nondesign regimes of nozzle 
operation. 

When the conditions are such that the pressure at the nozzle exit is 
equal to atmospheric pressure, the nozzle is said to be operating in the 
design regime. If the pressure at the nozzle exit is not equal to atmospheric, 
then it is said to be operating in the nondesign regime. 

There are two distinct nondesign regimes of nozzle operation: regime 
of underexpansion and regime of overexpansion. 

When operating in the underexpansion regime, the pressure in the gas 
jet at nozzle exit is greater than atmospheric pressure. Underexpansion 
can take place, for instance, when the rocket reaches some altitude where 
the atmospheric pressure becomes less than the pressure at the nozzle exit 
plane. Gas emerging from the nozzle with an excess of pressure expands 
in the atmosphere and mixes with the external medium. During this exter- 
nal expansion, the flow velocity does not increase, for all practical purposes, 
due to formation of vortices along the flow boundary. 

In the case of overexpansion, the pressure at the nozzle exit plane is 
less than external pressure. Under these conditions there is, naturally, a 
deceleration of flow aft of the nozzle, and the velocity is reduced to sub- 
sonic. The change to subsonic velocity, as will be shown below (page 229), 
is accompanied by a sudden change in pressure. As the external pressure 
increases, the pressure shocks approach the nozzle exit and, in the presence 
of a sufficiently Jarge excess of atmospheric pressure, enter into the nozzle. 

In this case, the normal operating regime of the nozzle is disturbed, since 
the formation of shocks within the nozzle leads to the flow separation at the 
walls and gives rise to strong vortices and large losses in kinetic energy. 
According to some data, the formation of shocks within the nozzle takes 
place if the external pressure is 2.5-5.5 times the pressure at the nozzle exit. 
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The overexpansion in the motor occurs when the motor operates at such 
altitudes and such regimes that the pressure of the surrounding atmosphere 
is greater than the pressure at nozzle exit (for instance, during the decrease 
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Fig. 6.16. Comparison of operating regimes of a nozzle. 
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of fuel consumption in the motor and the corresponding reduction in 
chamber pressure). 

More interesting in the evaluation of nozzle operating regimes 1s the 
comparison of these regimes according to the magnitude of the thrust 
developed by the motor. 

Let us follow through a comparison of three nozzles shown in Fig. 6.16. 
The first nozzle is shortened and the pressure, pe, at its exit, will be greater 
then the pressure of the surrounding medium, p. The nozzle operates in a 
regime of underexpansion. The second nozzle has larger dimensions and 
operates in a design regime (pressure p. is equal to the external pressure). 
Finally, the third nozzle has even larger dimensions, and the pressure at 
nozzle exit becomes less than the pressure of the surrounding medium. 
Here we have a regime of overexpansion. 
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If we return to formula (1.6) 
P = mw ie Se(De m P), 


then it is difficult to say in which of the three cases we will get greater 
thrust. In the first case, the exhaust velocity, w, will be less, but then the 
difference pe — p will have a positive value. In the third case, the exhaust 
velocity will be greater, but the difference p, — p becomes negative. 

The posed problem, however, is solved rather simply, if it is remembered 
that the thrust is the resultant of pressure forces distributed along the 
internal and external surfaces of the motor chamber. Fig. 6.16 shows the 
distribution of pressures near the exit plane of the nozzle in the three cases 
under discussion. As we see, the enlargement of the nozzle during under- 
expansion will lead to the increase of the axial component of the pressure 
forces. These forces are shown in the first sketch by dashed lines. In the 
regime of overexpansion, the atmospheric pressure dominates over the pres- 
sure of the exhaust gases and at the right section of the nozzle yields a 
force directed opposite to thrust. Thus it becomes obvious that the optimum 
operating regime of the nozzle, with constant consumption, is the design 
regime. 


2. CHARACTERISTICS OF Rocket MOTORS 


The relation of thrust to the operating conditions of the motor is called 
the characteristic of the rocket motor. The basic characteristics are altitude 
and throttling characteristics. 

The altitude characteristic is the dependence of the motor thrust on the 
altitude at which it operates with constant fuel consumption. In order to 
obtain the altitude characteristic, let us examine the formula for the thrust 
of the motor, Eq. (1.7) 


P = Po + S.(po = Ph). 


During constant fuel consumption G, the value Pe also remains con- 
stant. Therefore, the only variable in the thrust formula is the pressure 
pn. The total variation (increase) of the motor thrust with ascent from the 
ground, where p, = 14.7 lb/in.?, to an altitude where the atmospheric 
pressure can be considered as equal to zero, is 14.7 S, lb, if the area of the 
nozzle exit plane is expressed in in.?. Under these conditions (in a vacuum), 


Pes = Po + Sepo 
and thrust at sea level, (p, = po) 
Pa = Po. 


The altitude characteristics of motors will look like the curves shown 
in Fig. 6.17. 
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It will be interesting to determine the thrusts that motors will have 
with the same fuel consumption, G, the same area of the nozzle critical 
sections, and the same chamber pressures, but different areas of the nozzle 
exit planes, Se, and, consequently, different exit pressures, pe- 

From considerations presented in the previous section, it is obvious 
that a motor with a high design altitude, h (for which pe = pa), i.e., having 
a large exit plane area S,, will develop greater thrust in the upper regions 


p 





Por 


Pot 


hear Pro h 
Fic. 6.17. Altitude characteristic of a rocket motor. 


of the atmosphere (curve B in Fig. 6.17). A motor having lower design 
altitude will develop a greater thrust in lower regions of the atmosphere 
(curve A). 

That is why many motors of the long range and anti-aircraft rockets, 
which operate for the most part at high altitudes, and also aircraft engines, 
are designed in such a way that the pressure at the nozzle exit is lower 
than atmospheric pressure at sea level (10-11 lb/in.?). Nozzles of such 
motors are called high altitude nozzles. The relative change in thrust during 
ascent to an altitude depends on the design altitude of the nozzle, and for 
contemporary motors constitutes 12-15%. 

The other characteristic of a rocket motor is the throttling charac- 
teristic. 

The dependence of thrust on the variation of fuel consumption at con- 
stant altitude is called the throttling characteristic. Generally speaking, 
variation in fuel consumption in the combustion chamber of a motor varies 
the temperature and the value of the gas constant. However, these varia- 
tions usually are not taken into account when determining the throttling 
characteristic. 

The equation for throttling characteristic can be easily obtained from 
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the equation of thrust of a rocket motor. We take into account that for a 
given motor the exhaust velocity w. does not depend (at constant R and 
To in the chamber) upon consumption G, and the pressure at nozzle exit 
drops in proportion to consumption. The Jast statement stems from the 
fact that the pressure ratio, p./po, for a nozzle of given dimensions, i.e., 
with a given ratio S./S.,, is constant. In this connection, the pressure at 
nozzle exit can be expressed as a constant percentage of chamber pressure 


De = Xo. 
The chamber pressure itself, with accepted allowable constancy of Te 
and R, is proportional to fuel consumption G. Consequently, 
Pe = LG 
from which we get 
Pe = PenomU/ Groin 


where Gnom 1s the nominal (design) fuel consumption for a given motor; and 
Penom 1S the nozzle exit pressure at this consumption. 





Gnom G(R) 






Fig. 6.18. Throttling characteristic of a rocket motor. 


Taking into account the remarks made above, expression (1.6) for the 
thrust of a motor 


P = Gw./g + Sepe — Pr) 
can be written in the following manner: 


P= G(w./g + Se Denain. G now) = SeDn- (6.40) 


The obtained equation of the throttling characteristic represents an 
equation of a straight line, intersecting the ordinate at a distance—S,p, 
(Fig. 6.18). The slope of this line is determined basically by the value w,, 
and increases with an increase in Se. 
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In certain cases, making use of the fact that the chamber pressure is 
proportional to fuel consumption, the throttling characteristic is based 
not on consumption but on the value of chamber pressure. This is particu- 
larly convenient in determining experimental characteristics obtained 
during motor testing, since the chamber pressure is much easier to measure 
than fuel consumption. 

The linear property of the throttling characteristic determined by Eq. 
(6.40) is well confirmed experimentally, with fuel consumptions greater 
than nominal, and with lower consumptions when the overexpansion is 
not pronounced enough to allow pressure shocks to enter the nozzle. There- 
fore, the throttling characteristic is linear down to consumptions equal to 
30-40% of nominal. At lower consumptions, the relationship in Eq. (6.40) 
becomes invalid and the throttling characteristic takes the approximate 
form shown in Fig. 6.18 by a dotted line. At low chamber pressure a con- 
siderable increase begins in the degree of dissociation of combustion prod- 
ucts and there is a corresponding change in chamber temperature, which 
also reflects on the throttling characteristic of a motor. 


E. Cooling of Liquid Rocket Motors 


1. Heat EXcHANGE IN Rocket Morors 


The high temperature of the combustion products and the high velocity 
of flow in the chamber, especially along the nozzle, lead to rather intensive 
heat transfer from the gas to the walls of a rocket motor. 

This would seem to contradict what has been said before with respect 
to the processes taking place in the rocket motor. Before, we have con- 
sidered gas expansion as adiabatic, i.e., taking place without heat exchange 
between the gas and the walls. Now we speak of intensive heat transfer 
to the walls. However, the crux of the matter is that heat energy, trans- 
ferred to the walls, constitutes only a small part of the total heat content 
of the gas and cannot substantially affect the character of flow. However, 
the absolute value of this small part constitutes a very large quantity of 
heat which, without suitable cooling, cannot be absorbed by the walls 
without burning through and destruction of the walls. 

The heat transfer from gas to the wall takes place in two ways; by 
convection and by radiation. 

In contemporary rocket motors heat exchange by convection plays a 
basic role. In substance, it is that particles (molecules and atoms) of gas, 
during their random motion and flow along the walls, exchange heat energy 
with the walls. In this way, only a layer of gas flowing next to the wall, the 
so-called interface layer, participates in convection heat exchange. The 
dimensions of this layer are close to those of the boundary layer. 
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The amount of heat transferred to the wall depends on the number of 
molecule collisions against the wall and the difference in temperatures 
between the gas and the wall. When the temperatures are equal the gas, on 
the average, receives as much heat from the wall as it transfers to it. 

The effect of the number of molecular collisions with the wall on heat 
exchange is determined by the coefficient of heat transfer, a, which repre- 
sents the amount of heat transferred to a unit wall surface in a unit time 
with a temperature difference between the wall and the gas of one degree. 
The number of molecular collisions of a moving gas against the wall is 
proportional to the number of molecules in a volume of gas passing by the 
wall (i.e., density of the gas p) and gas velocity w. 

Therefore, the number of collisions, which means the amount of heat 
transferred to the wall, is proportional to the product pw. This last term 
represents the specific mass discharge, which is equal to the ratio of the 
mass rate of discharge to the area of the passage section 


pw = G/gs. 


The amount of heat transferred by a gas to a unit wall surface in a unit 
time, i.e., unit heat rate, is 


q = a(T ge — Tw), (6.41) 


where 7’, is the temperature of the wall surface wetted by the gas; T'gas 18 
the gas temperature; and a, is the heat transfer coefficient from gas to the 
wall. 

Let us note here that since the gas flow near the wall is decelerated, it 
has at all sections of the nozzle temperatures close to the temperature of 
stagnation and the chamber temperature To. The wall temperature is 
considerably less than the gas temperature. As a result of this the tempera- 
ture difference Teas — Tw is considerable. On the other hand, because of 
high pressures and high velocities the specific mass discharge of gas is high 
and the coefficient of heat transfer a, becomes very high in liquid propelled 
rockets. All of these circumstances result in a condition in which the walls 
of rocket motors absorb a considerably greater quantity of heat then the 
walls of any other heat engine. 

Both coefficient a, and the temperature difference Teas — Tue vary 
along the length of the motor. The coefficient of heat transfer a, is the more 
strongly affected of the two, due to change in specific mass discharge pw 
along the length of the motor. The highest values of pw and a, are reached 
in the critical section of the nozzle. Here, also, is the maximum value of 
the heat flow q (Fig. 6.19). In contemporary motors the value of convection 
unit heat rate in the critica] section reaches 37-10% Btu/ft?/hr. 
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During the increase in the chamber pressure and temperature, the heat 
flow toward the nozzle walls increases considerably. 

The heat flow due to radiation is added to the convection heat flow. 
The amount of radiation heat flow is considerably less than heat flow by 
convection. The greatest amount of radiation takes place where gas tem- 
perature is high, i.e., in the combustion chamber. Here, radiation heat flow 
is equal to (55-75) - 104 Btu/ft?/hr. The amount of radiation is less in the 
critical section because of the lowered gas temperature. 





Fic. 6.19. Distribution of heat flow and temperature in a hot wall along the length 
of the motor: 1, heat flow across the wall; and 2, temperature of a hot wall. 


The heat flow received by the walls within the chamber must be trans- 
mitted to the surface of the wall in contact with the cooling liquid. At the 
same time, the wall surfaces absorbing and giving up heat have tempera- 
tures related by the equation of thermal conductivity 


where Twe is the temperature of the wall surface wetted by cooling liquid; 
ô is the wall thickness; and à is the coefficient of thermal conductivity of 
the wall material. 

The quantity 6/A, by its sense, is the thermal resistance of the wall. 
The thicker the wall and the poorer the heat transmission of the material, 
the higher the thermal resistance of the wall. 

The heat flow must then be transmitted to the cooling liquid, having 
a temperature T.. The amount of heat flow is related to the temperatures 
of the surfaces wetted by the coolant and the liquid by the equation 


c= ae(T we = Ts) (6.43) 
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where a, is the coefficient of heat transmission from the wall to the cooling 
liquid. 

Eliminating Tw. from Eqs. (6.24) and (6.43) we will get, for the wall 
surface temperature in contact with the gas, the following expression: 


Lo = Te + q(1/a. F 5/X). (6.44) 


Temperature variation during heat transfer from gas to wall and then 
to the cooling liquid is shown in Fig. 6.20. 





Fig. 6.20. Temperature distribution in the boundary gas layer, in the wall, and in 
the cooling liquid. 


The temperature of the wall, Te, must not exceed that upper limit at 
which the wall material noticeably loses its structural properties. In this 
sense, it is desirable for the value 7', to be as low as possible. Temperature 
Te may be limited by increasing the cofficient of heat transfer and by 
reducing the thermal resistance of the wall. 

In addition, the coefficient of heat transfer to the liquid must be large, 
because the excessive rise of Tw: above the boiling point of the cooling liquid 
cannot be allowed in the cooling system of the motor. A slight increase 
above the boiling temperature is allowable, and even desirable, since it 
increases heat transfer coefficient a,. The heat transfer is intensified as a 
result of steam bubble formation at the wall surface. The bubbles are con- 
tinually washed away by the flow of the cooling liquid and are condensed 
in the cooler liquid, giving up the heat received from the wall. 

The values of heat transfer coefficients a, and a, depend on the mass 
flow of the wall wetting liquid. However, since the density of the liquid is 
practically constant, the basic factor determining the value æ. is the flow 
velocity of the cooling liquid. The intensive cooling of chamber walls 
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and nozzle in liquid propelled rockets requires high values of coefficient 
a,. In this connection, the velocity of liquid in a cooling tract must be 
considerable. 


2. Motor CooLting DESIGN 


The problem of motor cooling consists in preventing the wall tempera- 
ture from exceeding a certain limit during the entire operation of the motor. 
This problem is solved differently in various types of rocket motors. 

In powder motors there is no liquid which could cool the motor. In 
this case, therefore, the operating time of the motor is limited in order not 
to overheat its walls to a temperature at which they lose structural integrity. 
In order to reduce the rate of wall heating and increase possible operating 
time of a powder motor, heat insulating liners, which are deposited on the 
inner surfaces of the motor, are used. A similar method is used when it 1s 
desirable to prevent overheating of the walls of short duration liquid rocket 
motors. 

In liquid rocket motors, however, reliable cooling is made more difficult 
by higher combustion chamber temperatures. On the other hand, cooling 
design is facilitated by the fact that protection of the walls may be obtained 
by means of a plentiful supply, to the internal wall surfaces, of one of the 
fuel components. It is true that such use of the fuel components reduces 
the efficiency of combustion and consequently lowers the impulse of the 
motor. However, it is still possible to use it in designs requiring the utmost 
simplicity. 

Most often the cooling of a liquid propelled rocket is accomplished by 
one of the fuel components, circulated through the cooling jacket. This 
method of cooling is quite advantageous. It does not require any additional 
liquids and location of containers for their storage in the rocket, i.e., it 
permits a motor installation of minimum weight. Besides, the heat absorbed 
by the cooling liquid is returned to the combustion chamber and is not 
uselessly lost. 

However, there are certain difficulties in designing motor cooling in 
this manner. They are determined, primarily, by the properties and the 
weight consumptions of components. It may be that one of the components, 
for instance, liquid oxygen, is unsuitable as a cooling liquid. There may 
not be enough of the second component for absorption of the total heat 
transferred by the wall of the motor during its operation. 

Another difficulty consists in the fact that in order to obtain the neces- 
sary value of heat transfer coefficient a, it is necessary to create high 
velocities of the cooling liquid in the jacket. Therefore, the gap between 
jacket walls must be very small. It is especially difficult to insure the 
required velocity of the cooling liquid in motors of small dimensions. For 
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instance, in a motor of the booster rocket shown in Fig. 5.6 the width of 
the gap at the critical section had to be made only 0.25 in. Because of this, 
a machined insert had to be added into the cooling jacket and the space 
maintained by the installation of spacers. In large motors, the size of the 
necessary gap increases and technological difficulties are partly removed. 

It must be said that the high velocities in the cooling tract are accom- 
panied by hydraulic losses, and necessitate an increase in the delivery 
pressure of the cooling component. This, in turn, complicates the opera- 
tion of the delivery system and increases its weight. 

Returning to expression (6.44), we see that the lowering of Tẹ may be 
achieved not only by increasing the velocity of the cooling liquid, with a 
corresponding increase in æ., but also by decreasing the thermal resistance 
of the wall 6/X. From this, it follows that the thickness 6 must be as small 
as possible, and the coefficient of thermal conductivity of the wall material 
A, as great as possible. 

In increasing the thickness ô, temperature 7’, is increased so sharply 
that the wall strength is reduced. With the increase of heat flow to the 
wall, the wall must be made thinner and thinner. Therefore, the environ- 
ment of the wall material is peculiar in that at high rates of heat transfer it 
is impossible to increase the strength of the wall by increasing its thickness. 

In order to decrease the temperature of a gas surface, the inner wall 
may be fabricated from a material with high heat conductivity, for instance, 
copper or aluminum. However, it must be considered that heat conductive 
materials, as a rule, are less heat resistant. 

The indicated conditions, especially the necessity for making a thin 
inner wall, lead to great difficulties in fabricating the motor chamber and 
the nozzle. 

The inner wall of the motor is subjected to considerable force due to 
the difference in pressures between the cooling liquid in the jacket cavity 
and the gases in the chamber and nozzle. This difference in pressures is 
especially large at the end of the nozzle. Therefore, in order to obtain a 
sufficiently strong inner wall, it is necessary to reinforce it by some means. 
Serving as an example may be the layout of a motor in which the inner 
jacket has been designed in the form of spirally wound copper tubes of a 
small diameter for the passage of the cooling liquid. 

Small diameter tubes can withstand the difference in liquid and gas 
pressures well. The internal chamber pressure is taken by a thick steel 
jacket, protected from the high temperature by a ring of tubes. Difficulties 
of motor design connected with the presence of large heat flows force a 
search for methods of reducing them. The basic method consists of supply- 
ing an excess of combustible into a layer adjacent to the wall. For instance, 
in long range rockets, the combustible is supplied into the interface layer 
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by several rows of ports (Fig. 5.8). Something like a curtain of combustible, 
protecting the wall, is thereby created. The action of the combustible cur- 
tain consists not only in heating and evaporating, by which it takes heat 
away from the wall, but also in lowering the temperature of that layer of 
gas which is directly in contact with the walls of the motor. The protective 
action of the layer overrich with combustible is maintained for some time. 
Gradually, however, the boundary layer is diluted by the main flow of 
the combustion products and its temperature is raised. Therefore, in 
large motors it is necessary to have several rows of ports supplying excess 
combustible. 

The introduction of combustible through radial ports is not expedi- 
tious, since the greater part of it would pass into the gas stream at some 
distance from the walls. A more reasonable method of combustible supply 
is in a motor whose chamber is shown in Fig. 5.6. Here, the combustible, 
through small openings, enters a narrow gap between the head and the 
chamber and moves along the walls without immediately mixing with the 
main gas mass. 

The expenditure of combustible for the creation of the curtain certainly 
lowers the specific impulse of the motor, since part of the combustible is 
not fully burned and carries chemical energy with it away from the motor. 
Therefore, a combined system of cooling is often used in motors in which, 
by moderate internal cooling, the heat flow to the wall is reduced, and the 
heat is removed by external cooling with comparatively small hydraulic 
losses in the cooling tract. 

A considerable reduction in the expenditure of combustible for internal 
cooling may be expected if a porous material is used for the internal cham- 
ber walls to make the so-called sweating wall,* which permits a uniform 
passage of a very small amount of liquid along its entire surface, forming 
a continuous curtain. 

There are also proposals involving the cooling of the motor by means 
of water, which would not only be boiled, but completely evaporated into 
superheated steam. This steam may be used in the turbine of the turbo- 
pump and then condensed in the heat exchanger, cooled by one of the com- 
ponents. Such cooling systems can find applications only in very large 
motors. 


* Translator’s note: evaporative cooling. 


Vil. Forces and Moments Acting on the Rocket in Flight 


A. System of Forces Acting on the Rocket in Flight and the 
Differential Equation of Motion 


1. COORDINATES DETERMINING THE POSITION OF THE ROCKET IN SPACE 


The position of a rocket in space is determined, first of all, by the three 
coordinates of its center of gravity x, y, and z, in the so-called earth coor- 
dinate system. The origin of this system is taken at the launching point of 


$ 





Fig. 7.1. Ground system of coordinates. 


the rocket. For long range rockets the x axis is considered to be a straight 
line, tangent to the arc of the great circle connecting start with the target 
(Fig. 7.1). Axis y is directed upward, and axis z is perpendicular to the first 
two axes. 

Let us also introduce the so-called coupled or ballistic coordinate system. 
We will place its origin at the center of gravity of the rocket and relate the 
axes x’, y’, and z’ to the characteristic elements of a rocket. We will direct 
axis x’ along the axis of the rocket and will call it the longitudinal axis. We 
will place axis y’ in a plane of symmetry of the rocket and will direct it 
perpendicular to axis x’, and axis z’ we will direct from the left stabilizer 
to the right. We will call axis y’ the normal and axis z’ the transverse axis 
(Fig. 7.2). 

For long range rockets, the plane x’-y’ of the rocket in a trajectory 
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coincides or, more exactly, almost coincides, with the plane zy of the earth 
system. 

In order to orient the position of a rocket fully as a solid body in space, 
in addition to coordinates x, y, and z, we will also introduce three angles 






Center 
of gravity 


Fig. 7.2. Flight system of coordinates. 


which determine the direction of the ballistic and the earth systems of 
coordinates. 

The angle between rocket axis x’ and plane x-z, i.e., the angle of inclina- 
tion of the rocket axis with respect to the initial horizon, we will designate 





Target View along axis y! View along axis x! 


View along axis Z' 
Fic. 7.3. Angles of pitch ¢, yaw y, and roll y. 
by ¢. This angle is called the pitch angle of the rocket (Fig. 7.3). The angle 
between the axis of the rocket and the plane x-y we will designate by y. 


This is called the yaw angle. It represents the deviation of the rocket axis 
from the verticle plane passing through the target. 
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Finally, let us introduce the last angle determining the rotation of the 
rocket frame about the longitudinal axis; roll angle y between the axis y’ 
and the plane x-y. 

In this manner, three coordinates, x, y, and z, and three angles, ¢, y, 
and y (see Fig. 7.3) completely define the position of a rocket in space. 


2. Forces ACTING ON A ROCKET 


Let us examine the forces acting on the rocket in flight. 

In Fig. 7.4 is shown a rocket and the trajectory of its center of gravity. 
The vector of flight velocity, v, is directed along the tangent to the tra- 
jectory. 

The rocket axis, generally speaking, does not coincide with the velocity 
vector v, but forms with it angle a, called the angle of attack. Ordinarily, 
in a guided trajectory, angle a is comparatively small (a < 6°). 







Trajectory 
` of the rocket 


Center of gravity 


Center of 
gravity 


Fie. 7.4. Rocket angle of attack a and trajectory angle with respect to the horizon 8. 


Any body moving through air is subjected to the action of a system of 
forces distributed along its surface. The resultant of these forces is called 
the total aerodynamic force. The component of aerodynamic force along 
the tangent to the trajectory of the center of gravity of the body (or the 
projection of the total force onto the direction of the velocity vector) is 
called the frontal resistance force or, simply, drag. This force is always 
opposite to the direction of motion. 

The drag force is usually designated by D. 

The component of the total aerodynamic force along the normal to the 
trajectory (or its projection on a normal to the direction of velocity) is 
called the lift force and is designated by L. Fig. 7.5 shows forces D and L; 
also the thrust, P, and gravity weight force, Mg. 
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The system of aerodynamic forces distributed along the surface of a 
rocket, just like any system of forces, may be reduced to a force (in this 
case, the total aerodynamic force) according to the laws of mechanics, and 
a moment whose magnitude depends upon the point of force application. 






Center of 
gravity 


Mg 


Fic. 7.5. Forces and moments acting on a rocket (with the aerodynamic forces 
referred to the center of gravity). 


Referring the force system to the center of gravity of the rocket we will 
get, aside from the mentioned forces D and L, the resultant moment M 
(Fig. 7.5). Ordinarily, it is considered as a summation of stabilizing and 
damping moments 


M = M ot + Ma. 


These quantities will be discussed in more detail later. 

If the rocket has steering mechanisms (aerodynamic fins, jet vanes, or 
gimbaled motor chamber) a system of forces due to steering mechanisms 
must be added to the system of forces shown in Fig. 7.5. These forces, in 
the form of two force components and a moment, are applied at the hinge 
of the steering mechanism located at a distance x’, from the center of 
gravity of the rocket. 

One of the components, D, (Fig. 7.6), is directed along the rocket axis; 
the second, L., is perpendicular to it. 

The first force, by analogy with aerodynamic forces, is called the vane 
drag and is not a steering force. The second force, L,, is called the vane lift 
and constitutes the steering force. 

Moment M, is called the vane hinge moment. 

Fig. 7.7 shows the full force system acting on the rocket. 

In order to derive the equation of motion let us use D’Alembert’s 
principle and introduce inertial forces. 
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Tangential acceleration of a rocket (acceleration along the tangent to 
the trajectory) will be 


dv/dt = v. 


The corresponding inertial force has a direction opposite to the accelera- 
tion and is equal to Mv. 


Mh 


Ay 
Fig. 7.6. Steering forces. 


The normal acceleration, depending on the curvature of the trajectory, 
is equal, as is known, to v?/R, where R is the radius of the trajectory 
curvature. 

However, 


1/R = d6é/ds = d6/dt X dt/ds = 6/0, 


where @ is the angle of inclination of the trajectory (Fig. 7.3). 

Therefore, normal acceleration directed toward the center of curvature 
becomes equal to v@. The interial force directed in an opposite direction is 
equal to Mvé. Force Mvé6 is shown in Fig. 7.7. 

Finally, let us introduce the moment of inertia, equal to the differential, 
with respect to time, of the moment of momentum of the rocket, and 
opposed to the angular acceleration ¢ i.e., clockwise (Fig. 7.7), 


[J(0 + a) = (Jey. 


The moment of inertia J of the rocket about the central axis perpendicu- 
lar to the plane of the trajectory is also a function of time. Therefore, 
moment (Jó) can be written as a sum 


(Joy = Jet do. 
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The first addend is the moment of inertia. The second addend is the 
result in inertial moment change of the rocket with time. This addend, 
proportional to the angular velocity, may be included in the aerodynamic 





Frc. 7.7. Derivation of equations of motion. 


damping moment. In derivation of the equation of motion we will add it 
to the previously examined moment M. In analyzing the damping moments 
further on, we will return to this question again. 


3. EQUATION oF MOTION 


Projecting all forces applied to the rocket onto the tangent to the flight 
trajectory we will get 


Mi + D — (P — D,) cos a + Mg sin 6+ L, sna = 0. 
Because of the small angle a, cos a ~ 1. Neglecting L, sin a, we will get 
ù = (P — D — D,)M— — g sin 8. (7.1) 
Now we will project all forces onto the normal to the flight trajectory 
L + L, cosa — Mvé + (P — D,) sin a — Mg cos 0 = 0. 
With small values of a 
6 = (1/Mv)[(P — D,)a+ L + L.) — (g/v) cos 8. (7.2) 


In these equations the term g is the acceleration due to earth’s gravity 
at flight altitude. Finally, summing up the moments about the center of 
gravity we get 


J@é+M+L,D,+ M, = 0. (7.3) 
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The system of equations (7.1), (7.2), and (7.3) describes the motion of a 
rocket in one plane only and does not account for the possibility of three- 
dimensional motion. Nevertheless, these equations can be freely used in 
any case for determining the parameters of a rocket trajectory, especially 
a ballistic rocket. 

In cases in which the flight trajectory is a spatial curve, as, for instance, 
that of an anti-aircraft guided rocket pursuing a maneuvering aircraft, the 
system of equations becomes considerably more complex. In developing 
this system it is necessary to examine forces and moments acting along 
three coordinate axes and introduce additional coordinates determining 
the position of a rocket. Instead of three equations we will get six. 

It is important to note that to the equations of motion of a guided 
rocket, equations of guidance which relate steering forces and vane angles 
to time, or motion parameters, must be added. The number of these equa- 
tions and their form depend on the steering system and the method of 
guidance and steering of the rocket. 

In an ideal case for a long range ballistic rocket, the steering equation 
looks very simple. To Eqs. (7.1) and (7.3) we must add the condition 


gp=O0+a=fit). 


where f(t) 1s the given time function. 

The steering condition is determined by the fact that during the guided 
portion of flight the rocket is oriented in a definite manner relative to a 
gyroscope whose axis is fixed in direction. We will discuss the question of 
steering more fully in Chapter IX. 

In order to be able to integrate the equations of motion it is necessary 
to know how the variables change with time and, in general, on what they 
depend. 

It is very clear that, for a given rocket, the aerodynamic forces depend 
on the velocity and flight altitude. The motor thrust does not remain con- 
stant and changes (if for no other reason than that along the flight trajec- 
tory the external atmospheric pressure, which enters into the thrust 
expression, changes). The rocket mass diminishes with time, corresponding 
to the expenditure of fuel. Finally, even the acceleration of gravity, g (if 
we are speaking of high altitudes), must be considered as a variable 
quantity. 

All of these questions we will examine below. Let us start with the prop- 
erties of the atmosphere. 


B. Earth’s Atmosphere and Its Properties 


Let us examine the properties of the atmosphere, keeping in mind that 
the magnitudes of aerodynamic forces acting on the rocket in flight depend 
on these properties. 
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The basic parameter of the atmosphere, affecting aerodynamic forces, 
is the air density, p. Temperature T also has some effect, inasmuch as a 
change in temperature, as we already know, changes the speed of sound 


c= VkgRT, 


and, depending upon the ratio of flight velocity to sound velocity (Mach 
number), changes the character of the flow, which changes the magnitude 
of aerodynamic forces, as will be shown below. 

It may also be said that the velocity of sound depends on the suspected 
variation in the chemical composition of the atmosphere with altitude. 
This will reflect on the value of the gas constant R and adiabatic exponent 

k. However, within the altitude limits for which the 
äi aerodynamic forces must be taken into considera- 
> tion in calculating rocket trajectories, i.e., altitudes 

of 60-75 miles, there is no evidence of any substan- 
tial change in the composition of the atmosphere. 

Therefore, the magnitude of aerodynamic forces 
is primarily affected by parameters of the air state— 
density p, temperature T, and, related to them, pres- 
dF sure p. Let us examine the dependence of these quan- 
tities on altitude. 

Pressure p must decrease with increase in alti- 
tude, inasmuch as the amount of pressure is deter- 

_ mined by the weight of a column of air. Let us 

ere Ohanga in develop an equilibrium equation of an elemental 
air pressure with alti- : 

tude. column with the base area equal to dF, and a 

height dh, isolated from the atmosphere (Fig. 7.8). 

Pressure p acts on the column from below. Pressure p + dp acts on the 
column from above, corresponding to the addition of height increment 
dh. Let y be the specific weight of air at altitude h. The equilibrium condi- 
tion will be as follows: 


dh 


h 


dpdF + ydhdF = 0 
from which 
dp/dh = —y 
or 
dp/dh = —gp. (7.4) 
The derivative dp/dh is negative, which indicates reduction in pressure 


with increase in altitude. It is easy to realize that in the absence of an 
external heat supply the temperature will also drop. 
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During an upward displacement of some mass of air, an expansion takes 
place with a corresponding lowering of gas temperature. Conversely, air 
being displaced downward compresses and its temperature is raised. In 
this manner a temperature equilibrium is established in the atmosphere in 
which the lower layers will have a higher temperature then the upper layers. 
Certainly, this assertion is correct only to the extent that the heat supply 
to the upper layers from without (by means of solar radiation) may be con- 
sidered small. 

The temperature distribution with altitude may be considered as 
depending on a thermodynamic process which corresponds to the above- 
mentioned expansion and compression of the gas during its slow vertical 
mixing with conservation of thermal equilibrium. Suppose that this is a 
polytropic expansion with an exponent n which is independent of the alti- 
tude h 


p/Po = (p/po)”, (7.5) 


where po and po are the pressure and density of the air at sea level. 
Having accepted this law we have the means of establishing, in the 
first approximation, the law of change of pressure, density, and tempera- 
ture of air with altitude. 
Let us find, with the aid of Eq. (7.5), an expression for density p, and 
let us substitute this expression into the equilibrium equation, Eq. (7.4). 
We will get 


dp/dh = —gpo(p/po)™. 
Separating the variables 
p-™dp = —gpopo dh 
from which, after integration, we find 
[n/(n — Dp» = — gop ™"h + C. 


The constant C is determined from the condition that when h = 0 the 
value of p = po. Consequently, 


C = [n/(n — Ilp 
and 
prin = po”) = Gpo po!” x [(n = t)/n]h 


nj(n—=1) 
_ _ , p(n—1 
p= ml Pal n Ja 


po/po = gRTo 


or 


Since 
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the final variation of pressure with velocity will be in the form 


n—- 1 n{(n—1) 
= | 7 (am) J Fe) 


Density is related to height by the relationship 


ee | 1/(n—1) 
oai- Cay ae 


According to the accepted polytropic process 
T = T(p/po) d" 





or, considering Eq. (7.6) 
T = To — [(n — 1)/nR]h. (7.8) 


Therefore, temperature drop with altitude is a linear function. Tem- 
perature drop gradient is equal to (n — 1)/nR. 

If for air we take n = k = 1.4, i.e., assuming that in a state of tem- 
perature equilibrium the expansion and compression of air, during vertical 
mixing, is adiabatic, we will get 


(n — 1)/nR = (1.4 — 1)/(1.4 X 53.3) = 0.0054 deg /ft. 


1.e., the air temperature drops approximately 5.4° for each 1,000 ft altitude. 
Actually, in the lower layers of the atmosphere the temperature drops an 
average of 3.57° per 1,000 ft, which corresponds to n = 1.23, and addition 
of heat during expansion. 

The quantity n, however, does not remain constant and changes with 
altitude. Therefore, relationships (7.6) through (7.8) must be considered 
as approximations. Investigations of the atmosphere show that up to the 
altitude h = 35,000—40,000 ft, the developed laws of air parameter varia- 
tions with altitude agree with observations sufficiently well. At high alti- 
tudes there are sharp discrepancies from the developed laws, primarily in 
temperature distribution. 

Starting with an altitude of approximately seven miles, the temperature 
no longer falls and remains approximately constant and equal, on the 
average, to —67°F. Then, starting at an altitude of 12 miles, a rise in tem- 
perature is observed, to a maximum at a 30 mile altitude. Further on, there 
is another drop in temperature. At an altitude of 50-60 miles, the tempera- 
ture is approximately equal to —76°F. Therefore, the temperature varies 
with altitude, approximately as is shown in Fig. 7.9. 

The first temperature rise at altitudes of 12-30 miles is explained, at 
present, by the presence of an ozone layer at these altitudes which intensely 
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absorbs the short wave (ultra-violet) solar radiation. The subsequent 
lowering of the temperature may be explained in the same manner as the 
first. Finally, the temperature rise in the uppermost layers is tied in with 
the bombardment of the earth’s atmosphere by cosmic particles and solar 
radiation. Therefore, the temperature of these atmospheric layers must 
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Frc. 7.9. Curve of air temperature variation with altitude. Based on NACA TN 1200. 


fluctuate within a 24-hr period. At night it is lower, and during the day, 
higher; according to hypothesis it reaches a value of the order of hundreds 
of degrees. 

It must not be supposed, however, that such high temperature will in 
any way cause complications during penetration of this region of the at- 
mosphere. The temperature of the aircraft remaining a prolonged period 
of time at these altitudes will remain considerably less then the tempera- 
ture of the surrounding medium. Due to the unusual rarefication of the gas, 
the heat transfer from the air to the aircraft will be very insignificant, and 
the balance between the heat supplied by the gas and the heat lost by the 
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body through radiation will be established at a relatively low temperature 
of the aircraft structure. In this sense an incomparably greater danger is 
represented by temperatures resulting from high velocities of rocket 
motion in the atmosphere. We shall dwell on this question later. 

The lower layers of the atmosphere, below seven miles, where the above- 
developed equations hold, are usually called the troposphere. The layers 
above seven miles are called the stratosphere. The uppermost layers, start- 
ing at 37 to 44 miles, are called the ionosphere. Here, the air is strongly 
ionized. 

In calculating rocket trajectories, it is extremely important to know 
the law of density and pressure variation of air with altitude. 


TABLE VII.1. PROPERTIES or THE ATMOSPHERE* 


Tem- Pres- f Tem- : 
Altitude pera- sure Density Altitude pera Pr eee D ensity 
h, ft ture ratio, TALIO; h, ft ture TAUG, ratio, 
T,°R. p/p p/ po T SR; p/ Po p/ po 
0 518.4 1.000 1.000 35,000 393.6 0.235 0.3098 
1,000 514.8 0.957 0.9710 | 40,000 392.4 0.185 0.2447 
2,000 511.2 0.980 0.9428 | 45,000 392.4 0.146 0.1926 
3,000 507.8 0.896 0.9151 50,000 392.4 0.115 0.1517 
4,000 504.2 0.865 0.8881 65,000 392.4 56.12 X 1075 74.14 x 1073 
5,000 500.6 0.882 0.8616 | 80,000 392.4 27.41 | 36.21 
6,000 497.0 0.801 0.8358 | 95,000 392.4 13.39 17.69 
7,000 493.4 0.772 0.8106 | 110,000 412.6 6.58 8.27 
8,000 490.0 0.744 0.7859 | 125,000 472.9 3.49 3.82 
9,000 486.4 0.714 0.7619 | 140,000 533.3 1.99 1.94 
10,000 482.8 0.687 0.7384 | 160,000 613.7 103.50 X 10% 87.42 X 10 
11,000 479.2 0.660 0.7154 | 180,000 630.0 56.98 46.89 
12,000 475.6 0.636 0.6931 | 200,000 619.4 31.40 26.28 
13,000 472.0 0.610 0.6712 | 220,000 552.4 16.56 15.54 
14,000 468.6 0.587 0.6499 | 240,000 485.3 8.03 8.58 
15,000 464.9 0.564 0.6291 | 260,000 432.0 3.51 4.21 
20,000 447.1 0.460 0.5827 | 280,000 447.4 1.51 1.66 
25,000 429.3 0.371 0.4480 | 300,000 487.4 0.72 0.69 
30,000 411.3 0.294 0.3740 | 320,000 527.5 0.38 0.32 





* Based on NACA No. 218 and TN 1200. 


In the troposphere and, to a considerable extent, in the upper atmos- 
pheric layers, the parameters of the air state vary with time of day, time 
of year, latitude, and finally with respect to the general meteorological 
conditions. However, all parameters fluctuate about certain mean values 
which were determined as a result of atmospheric observations over a period 
of some years. The average parameters are taken as the parameters of the 
so-called international standard atmosphere (ISA), i.e., a certain arbi- 
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trarily accepted atmosphere whose parameters are not considered to vary 
with time of year, day, or latitude. 

For troposphere and stratosphere, specifically, the following constants 
and relationships are accepted, tying in the air parameters with the altitude: 

1. Air—ideal gas with the gas constant R = 53.3. 

2. Pressure at sea level po = 29.92 in. Hg; density pọ = 0.002378 
slugs/ft?; temperature Tọ = 518.4°R. 

3. Temperature gradient in the troposphere is constant and is equal to 
0.003566°F /ft. 

4. Temperature in the lower layers of the stratosphere is constant and 
is equal to —67°F. 

With these standards for the troposphere, the variation in pressure, 
density and temperature with altitude according to expressions (7.6), 
(7.7), and (7.8) have the form 


p = po(l — h/147,500)5-26 
p = po(l — h/147,500)4 26 
T = Ty — 0.003566 - h. 


Ordinarily, a table of variation of air parameters with altitude, which 
was obtained by extrapolation of the international standard atmosphere 
data to high altitude region, is used for ballistic calculations; see Table 
VII.1. 


C. Aerodynamic Forces 


1. AERODYNAMIC FORCE COEFFICIENTS 


Let us suppose that the rocket has two planes of symmetry, forming 
the rocket axis by their intersection. 

If the direction of a rocket flight coincides with its axis, the lift force L 
is absent. This obviously follows from the symmetry of flow. The lift force 
will appear only as the rocket axis begins to deviate from the tangent to 
the trajectory (from the direction of velocity) and form with it an angle of 
attack a. 

Therefore, the value of the lift force L becomes zero when a = 0. As 
opposed to the lift force, the resistance force is present at all angles of 
attack. 

We will begin the examination of aerodynamic forces with the resistance 
force D, inasmuch as in calculating flight trajectories of rockets the role of 
this force, compared to other aerodynamic forces, is more substantial. The 
first and more natural attempt to determine the resistance or drag was an 
attempt to express it in terms of ram pressure. 


AERODYNAMIC FORCES 213 


If we examine the reversed motion and, as a first approximation, con- 
sider the gas incompressible, then for the simplest bodies having a shape 
approaching that of a plate, the value of resistance force D can seemingly 
be obtained in the followmg manner. 


1 





2| 


Fra. 7.10. Transverse flow around a plate. 


Designating by v the velocity of undisturbed air flow, let us write 
Bernoulli’s equation (see page 168) for a gas stream at section 1-1, suff- 
ciently removed from the body (vı = v), and at section 2-2, at the front 
surface of the body (Fig. 7.10). 


200 + pi = 3P? + po 


But ve = 0. Therefore 


Pe — Pr = FPV? = FPV. 
If we assume that the pressure at the back of the plate is equal to the 
pressure of the surrounding medium, i.e., pı, then the drag force can be 
obtained by multiplying (pə — pı) by the area of the plate S 


D = (p — pi)S = Fpv's. (7.9) 


Experience, however, does not confirm the obtained relationship even 
in this simplest case. This is understandable. A simplifying assumption 
that ve = 0 was made in the derivation. This is true only at the forward 
critical point while the air flowing away from that point has velocity other 
than zero. Further, it was assumed that directly behind the plate the pres- 
sure was equal to the pressure of the undisturbed flow. This is also incorrect. 
The pressure here would be somewhat lower. 

From experiment, the value of D for a round plate is about 11% higher 
than indicated by formula (7.9). Also, it so happens that this deviation at 
relatively low velocities remains constant independently of flow velocity 
and plate dimensions. 
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Therefore, if we introduce the correction coefficient 1.11, then for a 
round plate the drag force is 


D = 1.11 X Epv®S. 


In the general case, introducing the correction coefficient Cp, which is 
called the drag coefficient, we can write 


D = Cp X 4S (7.10) 


where Cp depends on the shape of the plate. 

The drag of any body is determined according to formula (7.10). For 
axially symmetric bodies, the area S is the projected area on the plane 
perpendicular to the axis of symmetry (the median section), or any other 
so-called characteristic area. 

The freedom in the choice of area S is reflected in the value of Cp. 
Therefore, in each case when the numerical values of Cp are applied, it is 
necessary to define the area taken as the characteristic. 

The formula for the lift force is derived in a manner analogous to Eq. 
(7.10) 


L = Cr x pv? S (7.11) 


where Cz is the coefficient of lift force and S is the same characteristic area 
as in expression (7.10). 

At relatively low flight velocities, corresponding to Mach numbers not 
exceeding 0.5—0.6, coefficients Cp and Cz may be considered as independent 
of velocity. The indicated velocity limit until recently corresponded exactly 
to velocities attained by aviation, and Cp and Čz were defined as constant 
values for a given shape of an aircraft or some of its parts. 

Formerly, these coefficients were considered as coefficients of configura- 
tion. Later, however, it was discovered that they also depend on velocity. 
As the flight velocity approaches the velocity of sound the values of Cp 
and Cz increase sharply. At high supersonic velocities the values of Cp and 
Cz decrease, asymptotically approaching some constant value. Variation 
of Cp and Cz with respect to Mach number and the angle of attack a is 
shown in Figs. 7.11 and 7.12. At the same time it must be pointed out that, 
in practical cases, Cp and Cz also depend to some degree on other factors— 
linear dimensions of the body, density of the air, ete. 

In one way or another, in rocket technology, where it is necessary to 
deal with a full range of velocities from low to high supersonic, the aero- 
dynamic coefficients Cp and Čz are considered as functions of velocity v, 
or more exactly, Mach number M. In this sense, formulas (7.11) and (7.12) 
lose their basic content, which established the proportionality between 
aerodynamic forces and the velocity head. However, due to their simplicity, 
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Fic. 7.11. Dependence of drag coefficient Cp of a ballistic rocket on Mach number 
M and angle of attack a (coefficient Cp is referred to the median section of the rocket; 
rocket jet is absent). 





Fig. 7.12. Dependence of lift coefficient Cz on Mach number M and angle of attack 
a (coefficient Cz is referred to the median section of the rocket). 


convenience, and certain traditions established in aviation, they remain 
basic design formulas regardless of whether or not low or high velocities 
are considered. 

The mechanism of formation of aerodynamic forces is quite complex. 
Determination of D and L forces by computation, even in the simplest 
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cases, often constitutes an insoluble problem. Therefore, currently, the 
determination of aerodynamic forces (values of Cp and Cz) is done by 
approximate calculations which are used as a guide and which later on are 
corrected by wind tunnel tests. 

Since wind tunnel tests of full scale rockets are possible only for small 
rockets, the tests are made with models geometrically identical to the 
rocket being investigated. Inasmuch as we are speaking about high flight 
velocities, it is necessary during tests to maintain the condition of similarity 
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Fra. 7.13. Dependence of drag and lift coefficients of a ballistic rocket on the angle 
of attack at M = 2. 
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of Mach numbers. In other words, the model must be tested in the wind 
tunnel at the same Mach number for which it is necessary to obtain the 
actual Cp. It is also desirable to maintain the simllarity in Reynold’s 
number Re, which takes into account the simultaneous action of viscosity 
and inertial forces 


Re = pvl/p 


where u is the viscosity coefficient, and l is some characteristic dimension 
of a rocket. 

The fulfillment of the last condition, i.e., the equality of the Re for the 
model and the prototype, meets with considerable difficulties and, in wind 
tunnel practice, as a rule, it is not complied with. 
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The dependence of aerodynamic forces on the angle of attack a, is 
expressed by coefficients Cp and Cz as functions of a for constant S. For a 
symmetrical rocket it is obvious that this relationship will be symmetrical 
for Cp and unsymmetrical for Cz 


Cp(—a) a Cp(a) 
Cr(—a) = — Cla). 
Dependence of Cp and Cz on angle of attack a is shown in Fig. 7.13 for 
a V-2 rocket at M = 2, and also by the curves of Figs. 7.11 and 7.12. 


It is important to note that at small angles of attack the drag coefficient 
may be considered as independent of the angle of attack. 


2. COMPONENTS OF AERODYNAMIC FORCES AND SuBSONIC FLOW 


The character of flow and the mechanism of the formation of aerody- 
namic forces differs depending on whether the flight velocity is subsonic or 





Fig. 7.14. Normal and tangential aerodynamic forces acting on the surface of a 
rocket. 


supersonic. However, in both cases, the force of the frontal resistance can 
be split up into two basic components; component of friction and com- 


ponent of pressure 
D = D; + D 


Force D, is the resultant of resistance forces, p: tangential to the sur- 
face of the body, and force D, is the resultant of forces, pa, normal to the 
surface (Fig. 7.14). 


D, = fp sin ydS 
D, = l Dn cos YdS 


where dS is an area element of a rocket and y is the local angle of the normal. 
Forces p: and p, are referred to the unit area of the surface of the body. 
Further, we can write 
D: = Co: X ovS 
D, = CDn X tov? S 


and get from this the expression for the over-all drag coefficient in the form 


Cp = Cot + Con. 
@ 
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In the simplest cases the boundary layer theory enables the calculation 
of the quantity Cp: based on Reynold’s number Re = pvl/u. For instance, 
for longitudinal flow around a plate (Fig. 7.15), coefficient C'n: for turbulent 
flow may be determined according to the formula 


Cp: = 0.074/VRe 


which holds true up to values Re = 10’. 

The plate formula for Cp: enables the evaluation of the magnitude of 
the coefficient Cp: for any body, if it is assumed that the amount of the 
“wetted” surface determines the friction force, independently of the con- 
figuration of that surface. 

The effect of viscosity on drag is not limited to simple friction, which 
was just discussed. The existence of viscosity substantially affects the dis- 
tribution of normal pressures, which form the second component of drag 
D. Depending on the value of viscosity, especially with poorly streamlined 
bodies, separation of flow occurs in the boundary layer at higher or lower 
velocities, which directly reflects on pressure distribution. 
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Fig. 7.15. Longitudinal flow around a plate. 


The theory which assumes an ideal liquid having no viscosity leads to 
unexpected results. For instance, there is the so-called Euler-D’Alembert 
paradox: a symmetrical body (sphere, cylinder) during slow motion in an 
ideal fluid should not experience drag. For instance, a cylinder with trans- 
verse flow (Fig. 7.16,a) is subjected to pressure forces symmetrical with 
respect to a plane perpendicular to the direction of flow. These forces are 
balanced and do not yield a resultant. The existence of viscosity leads not 
only to the appearance of tangential forces p, but, what is even more 
important, changes the character of flow itself. In some area the flow 
separates from the surface of the body resulting in the appearance of strong 
vortices behind the cylinder. The pressure, acting on the cylinder from 
behind, becomes less than the pressure of the surrounding medium (Fig. 
7.16,b). The pressure forces form a resultant directed against the motion 
of the cylinder. 

Generation of drag in this case, and its absence in the first case, is under- 
standable from the energetic point of view as well. During the motion of 
the cylinder through an ideal fluid, the streams separating in front of the 
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cylinder merge behind it and remain motionless. Therefore, the force 
imparting motion to the cylinder in the liquid does not perform any work 
and does not transmit any kinetic energy to the fluid. This is not true for 





-3 
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Fig. 7.16. Development of drag at subsonic velocities. 


the second case. In this case, a wake of vortices is formed behind the 
cylinder which possesses kinetic energy. This energy is imparted to the 
particles of fluid or gas by the expenditure of work of the force which 
moves the body in a viscous medium. 


(a) (b) (c) 


Fig. 7.17. Examples of streamline shapes: (a) profile of an aircraft wing for subsonic 
flight; (b) aireraft fuselage; and (c) a faired cylinder. 


From what has been said, it follows that the boundary layer separation 
and formation of vortices increase drag. In order to reduce drag forces in 
subsonic flow, it is necessary to have a profile of smooth contour which 
would facilitate flow without separation, i.e., create a so-called streamlined 
configuration. Fig. 7.17 shows some streamlined shapes. 

It is important to note that all streamlined bodies have smoothly con- 
verging contours at the aft part which do not permit flow separation, and 
rounded forward portions. The sharp-nosed shapes are not advantageous 
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in subsonic flow. The sharp edge disturbs the flow and increases the danger 
of its separation. 

At moderate flight velocities, the friction drag of a streamlined body 
comprises the main portion of the total drag; whereas the shape drag 
(pressure) plays a secondary role. For poorly streamlined bodies, such as 
a sphere or a cylinder, the form drag has considerably greater value than 
friction drag. At high velocities, however, even for well-streamlined bodies, 
the pressure drag becomes a deciding factor. 

In examining the pressure diagram shown in Fig. 7.16,b, we notice 
that the pressure on the surface of a body may be positive or negative. 
In this connection the pressure drag for bodies of revolution (rocket air 


Pressure taps 





Fig. 7.18. Rocket model for pressure measurement. 


frame) is divided into two parts: frontal drag (drag due to pressure dis- 
tributed along the air frame) and base drag arising as a result of rarefica- 
tion at the base (base suction). The drag is divided in the same manner at 
supersonic speeds. 

During flight, the air ahead of the rocket is compressed. An increased 
pressure is created here, whose magnitude determines “the pressure drag.” 
This drag increases without limit with increase of velocity. A rarefication 
is formed behind the rocket—the base vacuum, which determines base 
drag. This component of drag also increases with velocity, not without 
limit, however. It cannot be greater than the value corresponding to the 
absolute base vacuum. 

The separation of pressure drag into frontal and base is usually done on 
the basis of wind tunnel tests. The model is pressure tapped; that is, holes 
are drilled in it at many points. Through these holes the pressure is ducted 
to monometers (Fig. 7.18). As a result of wind tunnel tests a complete 
picture of pressure distribution along the surface of the model is obtained. 
After this, it is not difficult to determine the frontal and base drags. 

Fig. 7.19 shows the pressure distribution along the surface of the rocket 
body (without considering fins). In the presence of the angle of attack a, 
the character of flow around the rocket will differ from that at an angle of 
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attack a = 0. In this case there arises a side flow of air from the zone of 
high pressure into the zone of lower pressure which results in additional 






HA 


Fic. 7.19. Distribution of aerodynamic pressure along the surface of a rocket. 


formation of vortices and an additional expenditure of energy for their 
formation. The corresponding additional drag force is called the induced 
drag. The additional air flow and vortices formed by it are shown by arrows 
in Fig. 7.20. 
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Fig. 7.20. Formation of vortices resulting in induced drag at a rocket flight angle of 
attack. 


The indicated division of drag forces into separate components is valid 
both for subsonic and supersonic flow. However, during a transition to 
supersonic flow, additional components of drag forces appear which sharply 
change the flow character and correspondingly change the approach toward 
a selection of the more rational aerodynamic profiles. 


3. THe NATURE OF SUPERSONIC FLOW 


Let us return to the derivations which were evolved in Chapter VI 
[formula (6.19)] for determination of sound velocity. 
An expression was obtained for the propagation velocity of a compres- 


sion wave. 
= (2+) Ap. 
p Ap 
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For very weak (sound) disturbances, when p + Ap = p, this expression 
becomes the formula for the velocity of sound 


c = VAp/Ap 


from which, with the assumption of adiabatic process for gas compression 
during acoustical vibration, we get 


c= V kp/p = V kgRT. 


From the comparison of expressions for the propagation velocity of a 
shock wave and velocity of sound, it follows that strong disturbances always 
propagate faster than sound. 

Let us assume that in a still gas medium there is a constant source of 
weak disturbances. Let us imagine that at some time this source created 
a local compression of gas. As a result of this a spherical sound wave is 
created which begins to travel uniformly in all directions with the speed of 





Fig. 7.21. Propagation of weak wave Fic. 7.22. Propagation of weak wave 
fronts in a stationary medium fronts in a medium which moves with 
velocity v < c. 


sound c. After time t, the wave front would be at distance ct from the 
source of the disturbance. If the source of disturbances sends out periodic 
signals, the fronts of all waves will be arranged in the form of concentric 
spherical surfaces with their centers at the source of disturbances (Fig. 
7.21). The propagation of these waves may be pictured as being analogous 
to the propagation of waves on the surface of the water from a pebble 
thrown into the water. 

Let us imagine now that the source of disturbances is immobile, while 
the gas medium moves relative to it with velocity v, and v < c. During 
the time ¢, between two successive signals from the source of disturbances, 
the first wave will be displaced by quantity vt and the centers of the spheri- 
cal waves will be displaced relative to each other by the same quantity vt. 
At the same time, inasmuch as v < c the center of each subsequent sphere 
will remain within the-confines of the preceding sphere (Fig. 7.22). 
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The picture is different when v > c, i.e., in the case when the flow 
moving past the body is supersonic. In this case the center of each succeed- 
ing sphere will be outside the surface of the preceding sphere (Fig. 7.23), 
and we will get a family of spheres having a conical envelope. 

In this way the character of propagation of disturbances differs, depend- 
ing upon the velocity of flow. 

While the flow velocity was less than the velocity of sound, the disturb- 
ances propagated in all directions, even though not uniformly but with the 
flow as well as against it. However, when the flow velocity became greater 





Fig. 7.23. Propagation of weak wave fronts in a medium which moves with velocity 
v> ec. 


than the speed of sound, the disturbances began to propagate in one direc- 
tion only; with the flow. The conical envelope—the so-called Mach cone— 
constitutes a boundary separating the region of points which can be reached 
by the sound disturbances from the region which cannot be reached by 
them. 

While outside the boundary of this cone we cannot receive signals being 
sent by the source of disturbances. The angle at the apex of the cone (Mach 
angle) depends on the ratio of flow velocity to velocity of sound. From Fig. 
7.23 it may be seen that 


sin x = c/v = 1/M. 


The angle x increases with a decrease in flow velocity, and when v = c, 
becomes equal to 90°. This means that the cone of disturbances becomes a 
plane. When this happens the signals transmitted by the source of dis- 
turbances reach any point located behind the source. Let us suppose now 
that the source of disturbances is capable of continually transmitting not 
only the weak, acoustical waves, but also strong shock waves. Let us see in 
what manner these waves will be propagated in supersonic flow. 

A shock wave, as we already know, propagates with velocity c,, greater 
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than the velocity of sound. Therefore, initially, the wave traveling away 
from the source of disturbances and having sufficient initial force can propa- 
gate against supersonic flow. As it propagates, this wave will become weaker 
and its velocity will be decreasing, approaching, as the limit, the velocity 
of sound. Under these conditions, the envelope of the family of spherical 
waves will not be a simple conical surface. It will be a surface resembling 
at its head a hyperboloid, and then changing into a cone of weak disturb- 
ances. Fig. 7.24 shows such an envelope. Where the strength of the shock 
wave is greatest (directly in front of the source), the envelope is shown in 





Fic, 7.24. Envelope of a family of spherical shock wave fronts moving at a velocity 
v> ec. 


heavier outline. The concept of ‘‘disturbance source,” used until now, in- 
cludes more than a special device for transmitting sound signals. Each 
point, each projection, of a body moving in the stream may be a source 
transmitting to the stream continuous disturbances. There may be sources 
of either weak or strong disturbances. 

Let us examine a supersonic flow around some blunt-nosed body (Fig. 
7.25). Each point on the body may be considered as a source of weak dis- 
turbances. Signals from this source may be propagated only within the cone 
of weak disturbances. 

The disturbances in the vicinity of points located on the front flat part 
of the body form a zone of increased pressure in front of the body. How- 
ever, if a zone of increased pressure with a free forward surface was formed 
in front of the body, then a shock wave was also formed, which has the 
property of propagation with velocity greater than the speed of sound. A 
wave appears in this manner which moves against the flow. However, it 
does not have the opportunity to move too far forward. As soon as the dis- 
tance between the wave front and the forward edge increases, the gas is 
able to flow outward from the region of higher pressure (as is shown by 
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arrows in Fig. 7.25). During this time, the strength of the shock wave is 
reduced and the velocity of its propagation diminishes. Thus, there will 
always be a shock wave in front of a body in a supersonic stream which 
moves forward with the velocity of the supersonic flow. The distance be- 
tween the wave and the body depends on the configuration of the body, 
which determines the amount of gas spillover, and on the velocity of flow. 

Outside the boundary of the nose part of the body the shock wave front 
becomes inclined and the intensity of the wave diminishes. As a limit, the 
inclination angle of the shock wave relative to the direction of velocity 
becomes equal to Mach angle and the wave itself becomes an ordinary 
acoustical wave. 





Fig. 7.25. Shock wave ahead of a blunt-nosed shell at supersonic flight velocity. 


The presence of shock waves is a characteristic peculiarity of super- 
sonic flow. In supersonic aerodynamics, shock waves are ordinarily called 
pressure shocks. 

Shock waves develop not only at the front but also at the back of a 
body in the air flow. 

Let us examine the aft portion of a body in a supersonic stream or, what 
is the same, a body flying with supersonic velocity (Fig. 7.26). 

Let us take a corner of the body base as a source of weak disturbances. 

Let the flow, approaching point O, have the velocity vı, with vı > c. 
The approaching flow will remain undisturbed up to line OA. The location 
of line OA is determined by the angle x; = arc sin (c/v,). Going around the 
corner the flow must expand. For a supersonic flow the expansion, as is 
known, is accompanied by an increase in velocity. Therefore, beyond the 
line OA the flow velocity will be increasing and the flow will be deflected. 
Beyond a certain line OB the flow will be completely deflected and we move 
along a new direction with velocity v: Angle x: = arc sin (c/v2); since 
v > vı then x: < xı. 
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The shown schematic is somewhat simplified, inasmuch as stream 
boundaries of the turned flow will not be parallel, but an inescapable con- 
clusion is that with a flow past an external corner we will not get pressure 
shocks. Lines OA and OB represent waves of weak disturbances. 

Let us follow the subsequent behavior of the turned flow. We will note 
that the flow must turn again by the same angle but now in the opposite 
direction, having met an identical flow which has passed around the lower 
side of the symmetrical body. 

The point of merging of the flows on the axis of symmetry we will desig- 
nate by O’, and will take it as a new source of disturbances. 
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Fic. 7.26. Schematic of a formation of a tail shock wave. 


The flow approaching the angle has a velocity və. The source of dis- 
turbances O’ makes itself felt only to the right of the line O'A’, which is 
parallel to OB. Beyond the line O'A’, the flow compresses. For a supersonic 
flow, this means reduction in velocity, to value v3, for instance. The Mach 
angle corresponding to this velocity is 


X; = are sin (c/v3). 


Since vg < vz then X; > Xe. In Fig. 7.26 line O’B’ is located to the left 
of line O'A’. 

In this manner, we arrive at a contradiction: the boundary, where the 
curvature of the flow should start, is located beyond the line where this 
curvature should end. This contradiction can be explained only by the 
sudden change in gas velocity, across the shock, from v: to v3, with the 
abrupt turn of the flow. During all this, pressure and other flow parameters 
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also experience change across the shock. The shock boundary is located 
between lines O'A’ and O’B’. 

The formation of the tail shock was presented in schematic terms. In 
actuality, directly behind the flat base there is a turbulent wake, and the 
source of disturbances, O’, is purely arbitrary. The shock is formed at some 
distance from the axis. Fig. 7.27 is a schematic of the actual location of a 
shock wave. 

The appearance of shock waves during supersonic flow gives rise to an 
additional resistance, called the shock resistance. Actually, the force moving 
a body with supersonic velocity must perform additional work to maintain 
shock waves. It is important to emphasize that work is not expended on 
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Fia. 7.27. Tail shock wave and turbulent wake aft of a flat tail of a shell flying with 
supersonic speed. 


transporting the waves (they move of their own accord) but on maintaining 
the constant intensity of the shock waves. If additional energy were not 
supplied so the waves at the expense of the moving force, their velocity 
would drop below that of sound and they would be left behind a rocket 
flying at supersonic velocity. 

Inasmuch as a zone of increased pressure is formed behind a compression 
wave, the shock resistance may be treated as the consequence of additional 
pressure which: acts on the surface of the rocket in the presence of nose 
shock. Obviously, the stronger the shock the greater the losses, and the 
greater the shock resistance will be. 

We have already noted (Fig. 7.24) that the strength of the shock 
(pressure drop across the shock) depends upon the angle which the wave 
front makes with the flow vector. It is clear that for a wave to travel 
directly against the flow (the so-called normal shock), it must be stronger 
than when it moves at an angle to the flow (oblique shock). In the limiting 
case, when the wave front forms a Mach angle with the flow, its strength 
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may be negligibly small. Since the shock wave drag at supersonic speed 
represents considerable, and at high velocities even the major, part of the 
total drag, it is necessary to take special measures for its reduction. In order 
to reduce the over-all drag in supersonic flow, the rocket must be given such 
a form that shocks, since they are inescapable, be oblique and the shock 
angle approach the Mach angle as closely as possible. The least acceptable 
are shapes which give rise to normal shocks. 


Fra. 7.28. Flow around a subsonic profile in supersonic flow. 


After what has been said it becomes clear that the so-called stream- 
lined bodies, with the least frontal drag in subsonic flow, become poorly 
streamlined in supersonic flow. Fig. 7.28 shows a representation of super- 
sonic flow around a subsonic profile of a wing. A curved shock, approaching 
a normal, is formed in front of the wing as a result of which the drag on 
such a wing becomes considerable. 

Thin, elongated, pointed profiles are most advantageous at supersonic 
flight velocities. 
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Fia. 7.29. Dependence of the shock shape on the nose cone angle, at constant flight 
velocity. 


For a pointed profile, or a cone-shaped rocket head, the leading shock be- 
comes oblique, and for very acute angles of the leading edge, approaches the 
wave of weak disturbances. Energy loss in such a shock will be minimum. 

The shock angle increases with an increase in the thickness of the lead- 
ing edge, or the cone angle. When the angle of the cone or the wedge of 
the leading edge becomes equal to the so-called critical angle, the shock at 
the central partition will become normal and will detach itself from the 
nose. The losses due to this shock will be considerably greater. 
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Fig. 7.29 shows the change in the shock wave shape with respect to the 
angle of the leading edge at constant flow velocity. Therefore, the mini- 
mum included angle of the rocket head is a condition for the minimum 





Fig. 7.30. Ogival shape. 


shock resistance. However, excessive elongation of the nose leads to the 
increase in friction drag and structural difficulties. Ordinarily, the forward 
sections of rockets and shells are given the so-called ogival form which is 
described by arcs of circles or parabolas (Fig. 7.30). 


4, CRITICAL Macu NuMBER 


Shock wave drag which is the result of compression shock also occurs 
at subsonic velocities which are close to sonic velocity. Under these condi- 
tions the local stream velocity, while flowing over wings, fins, or other parts 
of the aircraft or rocket, may exceed the speed of sound and shocks will 
appear in corresponding places. 

The flight velocity at which the local velocities approach or reach the 
speed of sound is called a critical Mach number and the appearance of the 
first shocks is called the critical shock wave point. 

In order to understand the appearance of shocks in subsonic flow first 
let us convince ourselves that a supersonic flow changes into subsonic flow 
only across a shock front. Flow, in which a gradual transition from super- 
sonic to subsonic velocity occurs, is unstable.* 

Actually, in order to reduce the velocity of supersonic flow, the area of 
the section must be reduced. After the flow velocity becomes sonic, the 
area of the section must be increased in order to reduce the velocity of 
flow further. Thus we arrive at the scheme of the reversed Laval nozzle 
(Fig. 7.31). A flow through the Laval nozzle has been described earlier 
(see page 176). 


* Translators note: This is called a transsonic region. 


230 FORCES AND MOMENTS IN FLIGHT 


Let us assume that the flow velocity changes gradually from supersonic 
to subsonic (see Fig. 7.31). The flow velocity will equal the sonic velocity 
at the narrowest or the critical section of the nozzle. 





Fig. 7.31. Smooth deceleration of supersonic flow in a reversed Laval nozzle. 


Each point on the nozzle wall, and even the gas itself, may be considered 
as a source of weak disturbances. Let us take any point located in the sub- 
sonic region of flow, for instance, point O in Fig. 7.32. The disturbances 
arising from this point will be propagated in all directions and to any dis- 
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Fig. 7.32. Propagation of weak dis- Fig. 7.383. Shock in front of diffuser 
turbances in a subsonic section of flow. entrance. 


tance along the flow but only as far as the critical section against the flow. 
And this will be the case with regard to all disturbances from all points 
located in the subsonic region of flow. As a result, the disturbances from 
various sources will accumulate in the vicinity of the critical section and 
the wave so created will have the ability to move against the flow with the 
velocity greater than the speed of sound. 
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At some section, depending on the configuration of the supersonic flow, 
the velocities of flow and the wave (shock) will be equalized and the loca- 
tion of shock will be stabilized. 

In practical cases, with a rigid diffuser, the described retardation of 
supersonic flow leads to the fact that the shock emerges from the confines 
of the diffuser and is located immediately before entrance, while subsonic 
velocity is established behind the shock (Fig. 7.33). 

Therefore, the transition from supersonic to subsonic velocity takes 
place across shock and is explained by the fact that a sufficient number of 
disturbance sources in the flow is always present to form an integral wave 
capable of motion against a supersonic flow. 

Let us turn now to the question of a critical Mach number. 

During flow around a body some streamlines first converge, and then 
diverge, as is shown, for example, in Fig. 7.34. With sufficient pressure drop 
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Fic. 7.34. Change in streamline section Fig. 7.35. Shocks formed at the surface 
during flow past a body (wing profile). of a bullet at high subsonic velocities. 


across these streamlines, there is a possibility of velocity increase to super- 
sonic, as in a Laval nozzle. Further on, however, inasmuch as the free 
stream velocity is below sonic, there must be a reverse transition to a veloc- 
ity less than the speed of sound. However, we already know that this 
transition takes place through a shock. It is in this manner that shocks 
appear at subsonic flight velocities. 

Fig. 7.35 shows shocks formed at the surface of a bullet in flight. These 
shocks are formed at velocities close to the speed of sound (M = 0.900 
and M = 0.946). 

In contemporary aviation, which attains subsonic and sonic flight 
velocities, the question of critical Mach number is quite serious. At the 
appearance of the first shocks, there is a sharp increase in aerodynamic 
loading and there is qualitative change in the flow characteristic with an 
accompanying redistribution of pressure forces along the surface of the 
aircraft. All this results in high structural overloading and a change in the 
degree of stability and controllability of the aircraft. 

The difficulties thus arising are overcome by investigation and solution 
of stability and controllability problems in the region of the critical Mach 
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number, on one hand, and the application of design configuration which 
postpones the occurrence of the critical Mach number, on the other hand. 
An example of the last, for instance, is the trend toward thin, symmetrical 
wing profiles, and swept back and delta wings. 

The postponement of the critical Mach number with respect to rockets 
whose trajectory, for the most part, takes place at supersonic velocities 
(and such rockets are in the majority) does not make sense. Here it is neces- 
sary to determine the optimum configuration from the point of view of 
minimum shock drag. 


5. ROCKET Jet EFFECT oN AERODYNAMIC FORCES 


The presence of the jet stream of an operating motor introduces certain 
specific peculiarities into rocket aerodynamics. 

Let us examine the jet effect on each of the drag components. The 
motor jet cannot influence pressure drag to any appreciable extent. At 
supersonic flight velocity such an effect is simply impossible and at sub- 
sonic speed, even if any is present, it is not noticeable. 

Friction drag at supersonic flight velocity does not change in the 
presence of the jet, either. At subsonic velocity it 1s somewhat increased 
due to the fact that the jet stream aspirates (sucks in) surrounding air and 
accelerates the flow along the surface of the rocket (Fig. 7.36). The accel- 





Fig. 7.36. Aspiration of external flow by the motor jet. 


erated flow near the rocket tail is at a lower pressure, which tends to in- 
crease the frontal drag, especially in the case where the rocket has a conical 
aft section or boattail. 

The main influence of the jet is on the base drag. There is no rarefied 
zone behind the nozzle opening of an operating motor. Therefore, the total 
rocket drag at subsonic as well as supersonic flight velocities includes only 
that part of the base drag which corresponds to the base area surrounding 
the nozzle. As an example, in calculating the frontal drag for a rocket of 
the V-2 type with an operating motor, only 50% of the base drag is included. 

At low flight velocities the indicated condition lowers the frontal drag 
by a considerable amount. The change in drag due to the aspiration of the 
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surrounding flow is neglected. At high velocities, the reduction of the frontal 
drag due to the presence of the jet becomes unnoticeable, since the base 
drag plays a small part, and the shock drag becomes the main contributor. 
Therefore, in calculating rocket trajectories, the frontal drag varies 
depending on whether, all other conditions being equal, we are concerned 
with powered or coasting portions of the trajectory, 1.e., those portions of 
the trajectory with an operating or nonoperating motor. It must be noted 
that in general the thrust of a rocket exceeds many times the frontal drag 
and that in evaluating aerodynamic forces acting on the rocket there is no 
need for great refinement. In this sense the rocket differs considerably from 
an aircraft whose propeller thrust 1s equal to drag at level cruising speed. 


6. AERODYNAMIC HEATING OF Rockets AT Hicn FLIGHT VELOCITIES 


The aerodynamic effect of the medium on the air frame of a rocket in 
flight is not manifested by force alone. At high velocities there develops a 
thermal effect which manifests itself in the heating of the rocket surface. 
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Frc. 7.37. Determination of stagnation temperature. 


From the above-developed energy equation (in Chapter VI), it follows 
that during gas flow in a pipe, or around any body, there is a conversion of 
energy from one type to another. Kinetic energy is converted into heat and 
conversely. Reduction in kinetic energy (reduction in speed) necessarily 
leads to a rise in gas temperature and consequently in the temperature 
of the body. 

The gas velocity may be decreased either by contact with a body around 
which it flows, or as the result of boundary layer friction. It is obvious 
that the rise in gas temperature is independent of the source of the velocity 
loss. 

Let us examine a body in a gas flow (Fig. 7.37). 

At the leading edge of this body the flow divides into two parts. A por- 
tion of the streamlines goes up and another portion, down. It is obvious 
that at some point A at the leading edge of the body, the flow velocity is 
equal to zero. Here, there is a complete flow stagnation and all of the kinetic 
energy of the gas is converted into its heat content. This point is known as 
stagnation point, or critical point. It is not difficult to determine the gas 
temperature at this point. 

Let us separate out of the total flow a gas streamline AB, where point B 
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is located in the free stream, and let us write the energy equation (6.10) for 


this streamline. 
ve k (AÀ atk (?) 
2 +4), =- 2 Vea pja 


p/p _ gRT, 


Since 
then 
si k 
ht r URT)e = "E + RT) a. 


But vg = v and va = 0; therefore, the increase in temperature in a 
decelerating flow from point B to point A will be 


AT = T, — Tg = (k — 1)/kgR. 


The temperature at the critical point Ta represents the stagnation 
temperature. 
For air, k = 1.4, R = 53.3, and 


k — l v? 
kgR” = 12 x 10° 


for this the stagnation temperature is 
Taag = T + 0/12 X 10°. 


At low velocities, the stagnation temperature is not high. However, at 
supersonic speeds attainable by rockets it may turn out to be considerable. 
For instance, at flight velocity of 5000 ft/sec 


Tag = T + 5000?/12 X 10° = T + 2080°. 


At an ambient temperature of T = 495°R (35°F), this amounts to 2575°R 
(2115°F). 

The rocket temperature will be, naturally, noticeably below the stagna- 
tion temperature due to the partial heat transfer into the structure of the 
body and radiation to the atmosphere. However, it remains sufficiently 
high to require suitable design measures. In particular, long range ballistic 
rockets have nose tips made of heat resistant material (quartz, graphite). 

If the stagnation of the gas is not complete, then, understandably, the 
temperature of the stream will be lower than the complete stagnation 
temperature. For well-streamlined bodies, stream temperature in the 
boundary layer equals 40-70% of the full stagnation temperature. 

The heating of the rocket airframe, due to air stagnation near the rigid 
skin, is the basic obstacle in attaining supersonic flight velocities in the 
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atmosphere. This obstacle may be overcome by the development of highly 
heat resistant materials and the introduction of artificial cooling of the more 
critical places, as is done at the inner walls of the combustion chamber and 
nozzle of a liquid propelled rocket. However, the logically necessary step, 
and of the first order of importance, is, in all cases, to conform to the proper 
aerodynamic shapes which do not permit stagnation to occur. 

If we return to the problem of the ejection of combustion products 
from the nozzle of a rocket motor, it must be pointed out that stagnation 
takes place there too, with a resultant increase in temperature in the bound- 
ary layer along the nozzle walls. This, as has been mentioned earlier, leads 
to an Increased heat flow across the wall. 

In case of a complete gas stagnation in the nozzle of a liquid propelled 
rocket, and in the absence of heat losses, the kinetic energy of the flow is 
converted back into the heat content which the gas had in the combustion 
chamber. Consequently, in the absence of heat losses, the stagnation tem- 
perature of the stream in the nozzle of a liquid propelled rocket will be equal 
to the temperature in the combustion chamber. Such gas stagnation, and 
corresponding heating, take place, for example, at the leading edges of jet 
vanes. 


D. Stabilizing and Damping Moments 


1. STABILIZING AERODYNAMIC MOMENT 


It already has been mentioned above that a system of aerodynamic 
forces, distributed along the surface of a rocket, can be referenced to any 
point on the rocket in the form of a resultant force and a moment whose 
magnitude depends on the location of the point to which a system of forces 
is referred. 

If the system of aerodynamic forces is referred to the center of gravity 
of a rocket, then at an angle of attack other than zero, we will get, besides 
the L and D forces discussed above, a resulting moment—the so-called 
restoring or stabilizing moment M. acting to reduce the angle of attack 
a (Fig. 7.38). 

The magnitude of the stabilizing moment is determined by a formula 
analogous to the formulas used to determine drag force D and lift force L 


M st = Cm x 1 ov2Sl (7.12) 


where p, v, and S are air density, flight velocity, and characteristic area 
respectively; l is some characteristic linear dimension (ordinarily rocket 
length); and C,, is the dimensionless moment coefficient. 

At low subsonic flight velocities, the term Cm is independent of velocity, 
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At high velocities Cm, like coefficients Cp and Cz, is considered as a function 


of velocity. 
Coefficient Cm, in a first approximation (at small a), is proportional to 


the angle of attack 
Cae = Cra: 
The value of Cm“ unlike other aerodynamic characteristics, does not 
depend on flight velocity alone. Inasmuch as aerodynamic forces are 


referred to the center of gravity of the rocket, and the center of gravity 
itself moves along the axis as the fuel is expended, the value of C,,% also 
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Fic. 7.88. Reference of aerodynamic Fic. 7.39. Reference of aerodynamic 
forces to the rocket’s center of gravity. forces to the rocket’s center of pressure. 


depends on the internal distribution of the rocket mass. Because of this, 
it is more convenient to select some other characteristic instead of Cm” as a 
measure of the stabilizing moment. 

Along the rocket axis a point exists so located that, when the system of 
aerodynamic forces is referred to it, we obtain a moment equal to zero. This 
point is called the center of pressure. 

The center of pressure may be considered as a point of intersection of 
resultant aerodynamic forces with the rocket axis (Fig. 7.39). The location 
of the center of pressure depends exclusively on the external aerodynamics 
of the rocket. 

The relative location of the center of pressure and the center of gravity 
is important for stabilizing the rocket in flight. 

If the center of pressure is located behind the center of gravity, then, 
during deviation of the rocket axis from direction of flight, the aerodynamic 
forces create a moment which reduces the angle of attack and restores the 


STABILIZING AND DAMPING MOMENTS 237 


initial direction of the axis. If the center of pressure is located ahead of the 
center of gravity, the aerodynamic forces, during axis deviation from direc- 
tion of flight, will create a moment increasing the angle of attack. In the 
second case it is said that the rocket is not stabilized or is statically un- 
stable. In order to displace the center of pressure aft, the rocket is equipped 
with fins. An unfinned rocket, as a rule, is statically unstable. 

Fig. 7.40 shows the relationship of the center of pressure location of a 





Fig. 7.40. Shift of center of pressure location on the rocket with flight Mach number. 


long range ballistic rocket to flight velocity (Mach number). This graph, 
however, is not very informative, since it does not give a concept of rocket 
stabilization. It is preferable to show the relationship of the locations of 
center of gravity and center of pressure together with the function of time 
of rocket flight. Fig. 7.41 shows this relationship for the same rocket. 

It is seen from the included curves that the center of pressure always 
remains behind the center of gravity of the rocket during the guided flight 
(t = 64 sec). Therefore, the rocket is statically stable during its entire 
guided flight. 

The stability of a rocket is evaluated by the degree of stabilization, or, 
as it is called, the margin of stability—the ratio of the distance between 
the center of pressure and the center of gravity to the total length of the 
rocket. This value for finned rockets varies within the limits of 5 to 15%. 

It must be noted, however, that stabilization by means of fins cannot 
be effective during flights beyond the atmosphere. Fins retain their useful- 
ness only during the atmospheric portion of the rocket trajectory, and, in 
some cases, may be altogether discarded. 

In constructing a graph similar to the one shown in Fig. 7.41, the loca- 
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tion of the center of gravity is determined by the simple calculations of 
weight and distribution of a fuel remaining in the rocket at a given moment. 
The location of the center of pressure is determined by approximate 
calculations with the aid of model tests in wind tunnels. 
It is possible to determine the distribution of pressure forces along the 
surface of a rocket with the aid of a model equipped with pressure probes. 
Referring the system of forces thus obtained to some point, it is possible to 
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Fra. 7.41. Shift in location of center of gravity and center of pressure of a rocket with 
time during powered phase. 


find the moment and the total aerodynamic force, and then the location 
of the center of pressure. The moment and the aerodynamic forces may also 
be determined directly with the aid of an aerodynamic balance. Then the 
location of the center of pressure may be found by so selecting the point 
of model support as to have a zero overturning moment. 


2. DAMPING MOMENT 


Damping moments are moments whose magnitude depends on the 
angular velocity of rocket rotation about some axis. Damping moments may 
be considered relative to the longitudinal and the two transverse axes. 
Damping moments are opposed to the direction of rotation, and, in the 
first approximation, are proportional to the angular velocity. 

The development and magnitude of the damping moments depend on 
the conditions of flow around the rocket by an external stream, and on the 
conditions of flow of liquids and gases within the rocket and motor. Corre- 
spondingly, the moments are differentiated as external aerodynamic, and 
internal damping moments. 
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The external damping moment is obviously the result of a simple resist- 
ance of the air to the turning of a rocket. 

During the rotation of a rocket with an angular velocity w, the flow 
past a point located at a distance x from the center of gravity will take 
place with a change in the local angle of attack by the amount 


ar = wxr/v 


(Fig. 7.42). Because of this, a local aerodynamic moment develops in the 
opposite direction from rotation. The damping moment is determined by 
the summation of the elemental moments along the surface of the rocket. 
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Fig. 7.42. Development of an additional local angle of attack during rotation of a 
rocket. 


An expression for calculation of the damping moment can be derived 
by analogy with the expressions which were derived earlier for the aero- 
dynamic forces D and L and the moment M,. We will consider the total 
damping moment as being proportional to the velocity head 3pv?, charac- 
teristic area S, and ratio w/v. Also, introducing from dimensional consider- 
ations the square of the rocket length l, we get 


2 
Ma = Ca s x ~ SP, (7.13) 


The dimensionless damping moment coefficient Cz, like all the other 
aerodynamic coefficients so far encountered, is independent of flight veloc- 
ity at low speeds. At high velocities this quantity is considered as a function 
of velocity. Fig. 7.43 shows, as an example, a graph of the relationship of 
coefficient Cz to Mach number M for a V-2 rocket during rotation about its 
transverse axis. 

The internal damping moment is due to the Coriolis acceleration which 
occurs during the change in direction of the liquid flow through tanks and 
tubing of the rocket, and gas flow through the chamber and the nozzle of 
the motor. These moments are easily determined if we assume that the 
above-mentioned flows completely follow the turning of the rocket. Let us 
take an element of flow with the length dz’ and cross section S,, at a dis- 
tance from the center of gravity x’, and, consequently, with a mass pS,dz’. 

As is known, the Coriolis acceleration equals 2ww, where w is the flow 
velocity in the rocket; and w, the angular velocity of the rocket. 
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The elemental moment of the Coriolis force is 
aM = 2wwpS,dx’ex. 


The damping moment is determined by integrating this expression 
between the limits of 2’. 


Ma = 2w f. j pwS,x'dz’. (7.14) 
But, from the conditions of constant discharge 


pw, = m = constant. 


Therefore 
a yy , , 
Ma = 20m w T dz = wm(x3"2 — Xo’). 


For this expression, limits of integration x2’ and 2;’ in a liquid fuel rocket 
must be established between the surface of the component in the lower 


ad 
E 
N 
Zi 
Æ 
E 
BE 
= 
E 
HE 


ptt tt | Yt Et 





Fig. 7.43. The relationship of damping moment coefficient to Mach number during 
rotation of a rocket about its transverse axis. 


tank and the exit plane of the nozzle (see Fig. 7.44). In this interval the 
sum of the discharges of both components is equal to m. To the obtained 
expression it is necessary to add a similar one in which the total discharge 
m is replaced by the mass discharge of the upper component m, and the 
integration performed between z; and 2,’ (see Fig. 7.44), i.e., between the 
surfaces of the upper and the lower components. Finally, we will get 


Ma = wlm(x3” a £z?) -+ Mlz? = x17)]. 
Part of the mass moves toward the center of gravity and not away from 


it. Because of this the damping moment will be partially reduced. 
To the obtained expression may be added moment Jọ, proportional to 
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the angular velocity ġ, which is isolated from the moment of inertia (see 
page 204). 

The rate of change of the moment of inertia J of a rocket is determined 
from simple considerations. 





Center of gravity 


Fig. 7.44. Determination of the damping moment developed by the motor jet. 


During time At, the distribution of mass in a rocket will change only in 
the vicinity of the surfaces of the first and second fuel components (see 
Fig. 7.44). Consequently 


AJ = (— myx,” = Mtz ?)At 
and 
J= —m ny? = Mok", 


where m, and m: are the discharge rates of the first and second components. 
Adding Jw to the above-obtained expression for the moment, and taking 
into consideration that m = mı + me, we obtain 


Ma = wm (z? = 2x1?) + Malz? = 22'2)]. 


During rocket flight in dense regions of the atmosphere, and with the 
usual relationship between ¢ and ¢, the damping moment of a finned 
rocket constitutes approximately 10% of the stabilizing moment. Further- 
more, the greater part of this 10% is due to the damping moment of the 
external aerodynamic forces and only a small part is due to the internal 
moment determined by the last equation. Therefore, in calculating rocket 
trajectories in the dense atmospheric regions the damping moment due to 
internal forces may be neglected. Only during rocket flight beyond the 
limits of the atmosphere, where aerodynamic forces are nonexistent, does 
this moment become independently significant. 


E. Steering Forces 


Let us examine the gas (jet) vanes as the basic functional steering 
members. 
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The jet vane in the stream of a motor is acted upon by the gas dynamic 
forces of the flow in the form of the frontal force D,, the lift force L,, and 
the moment M, (Fig. 7.45). 
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Fig. 7.45. Jet vane in the stream of a motor. 


A vane in the gas stream behaves in a manner similar to an aircraft 
rudder in the airflow. 

The frontal and lift forces on the vane, and also moment M,, may be 
determined by the usual aerodynamic formulas 


D, = Cp» x tow? S; 
Ly, = Crs X tpu? S 
M, = Cun X pw? Sl 


where w is the gas flow velocity past the vane; p, the density of gas in the 
stream; S, the characteristic vane area; and l, the characteristic vane 
length. 

Coefficients Cpe, Cra, and Cy, in this case may be considered as values 
independent of the flow velocity w, inasmuch as the flow is considerably 
supersonic (M = 2.4), and at these velocities aerodynamic coefficients 
vary insignificantly with change in velocity. Finally, the exhaust gas 
velocity through the nozzle for all liquid fuel motors is approximately 
the same, and usually equals 6500-7200 ft/sec. 

Coefficients Cp,, Cz», and Cy», as do all aerodynamic coefficients, 
depend on the shape of the vanes and the angle of attack, i.e., on the vane 
turn angle 8 in the stream (see Fig. 7.45). It must be said that this angle 
varies through a considerable range. For instance, jet vanes in long range 
ballistic rockets are capable of turning in the stream up to +25°, while the 
working range of the vanes is within the interval +15°. This leads to the 
fact that the forces acting on the vane during the guided portion of rocket 
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flight vary between considerable limits. The typical dependence of coeffi- 
cients Cp», Cz», and Cy, on the vane angle of attack 8 is shown in Fig. 7.46. 
The magnitude of the operating angle 8 in flight may be considered as 
consisting of two addends. The first addend is determined by the pro- 
grammed flight of the rocket to the target, i.e., by that trajectory which 
must be imparted to the rocket under the ideal conditions of undisturbed 
flight. This vane angle of attack may be calculated beforehand. The 
necessary steering forces may be calculated according to the shape of the 
trajectory and, from them, the necessary angles may be calculated. 
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Fig. 7.46. Relationship of coefficients Cpr, Ct», and Cu» with the vane angle of 
rotation 8. 


The second addend is an angle, resulting from the reaction of the 
stabilizing mechanism and, consequently, of the vanes, to random dis- 
turbances. They include factors which constantly act during the trajectory, 
such as noncoaxiality of the aerodynamic and reactive forces and variable 
wind gusts. Their exact magnitude cannot be determined beforehand. It is 
determined on the basis of allowances for fabricating malalignment and 
from statistical data of observed atmospheric conditions. Allowance is 
made for additional vane motion, above that which is determined from 
programmed conditions, to compensate for possible wind gusts. 

The first of the two mentioned addends of the operating angle 8 is basic, 
and results in systematically acting forces taken into consideration in 
ballistic calculations. Therefore, the steering forces D,, L,, and M, are 
examined, first of all, as definite functions of time, inasmuch as angle 8, 
as a function of time, is known. Concerning the second addend of angle 
8, it leads to random, unsystematically acting forces, which are not con- 
sidered in ballistic calculations. 

Let us examine the forces D, and L,, and the moment M, separately. 

The total jet vane force D, is not a steering force and is considered as 
thrust loss due to jet vanes. This loss is considerable. For instance, for the 
V-2 rocket, when £ = 0, it constitutes approximately 1400 lb (350 Ib per 
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vane) and increases with vane rotation. A factor which somewhat reduces 
the magnitude of D,, for graphite vanes, is the erosion of vanes at the sur- 
face during the operation of the motor. 

The moment M, consists of the sum of the stabilizing and damping 
moments proportional to the vane angle of rotation 8 and angular velocity 
B, respectively. If, to this moment, we also add the moment of inertia of 
the vane J,@, then we will get the magnitude of the hinge moment M,, 
(see page 205): 


M, = aß + aß T aĝ, (7.15) 


where ao, ai, and as are the coefficients of proportionality. 

For graphite vanes, these coefficients change somewhat during the 
operation of the motor due to erosion of graphite at the surface. The most 
substantial change is in the coefficient a, since the graphite erosion takes 
place essentially at the leading edge, which is subjected to the strongest. 
flow (Fig. 7.47,a). This displaces the vane center of pressure considerably 
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Frc. 7.47. Variation in the graphite vane configuration during the operating time of 
the motor. 


aft (coefficient ao increases). In order to avoid this, the graphite vane is 
given the shape shown in Fig. 7.47,b. Here, the erosion takes place on both 
sides of the vane hinge line; therefore, the center of pressure is not displaced. 

The hinge moment included in Eq. (7.8) is considerably smaller than 
the magnitude of L,2,’, and it may be neglected in the equation of motion. 
The magnitude of the hinge moment is important only for evaluating the 
necessary power requirements of the steering mechanisms which rotate 
the vanes. 

The jet vane lift force L, creates the main steering moment L,2,’, 
which turns the rocket in the required direction; see Eq. (7.3). The lift 
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force L, itself is usually included in the aerodynamic lift force L; see 
Eq. (7.2). This is done in the following manner. 

Let us assume the rocket flies with a constant programmed angle g. 
Then, from Eq. (7.3), neglecting the term M,, we get 


M + Lz,’ = 0. 
But the moment 
M = Le 


where e is the distance from the center of pressure to the rocket center of 
gravity. Therefore 
L, = —Le/z,’. 


In this manner Eq. (7.2) is converted to the form 
ô = (1/Mv)[(P — D,)a + L(1 — e/x,’)] — (g/v) cos 6. (7.16) 


If the steering of the rocket is performed by the rotation of the combus- 
tion chamber (Fig. 7.48), then the thrust loss D, will be very small, inas- 





Fie. 7.48. Development of the steering force with the canting of the motor chamber. 
much as the angle of chamber rotation £ is small while the magnitude 
D, = P(1 — cos B) = P(p?/2). 
The force L, in this case is determined according to the formula 
L, = P sin B = PB. 


F. Change in Thrust and Rocket Weight Along a Trajectory 


Of all the forces entering into the equations of motion (7.1), (7.2), and 
(7.3), we have left only two which were not examined: the thrust of the 
motor, and the weight of the rocket. The main reason for the change of 
thrust during flight is the pressure change in the surrounding medium. In 
gaining altitude, the atmospheric pressure changes and, together with it, 
the axial component of the external static pressure. This follows from 
expression (1.6) 


P=mw+S8S.(pe — p). 
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In other words, as the rocket gains altitude, the thrust increases to corre- 
spond to the drop in atmospheric pressure. 

For liquid fuel rockets, a second reason exists for some change in thrust 
during flight. This is the change in the discharge of fuel components due to 
variable conditions of delivery. 


Fig. 7.49, Pressure determination at 
the intake of the turbopump assembly. 
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The pressure of fuel components at the entrance to the turbopump 
assembly is made up of pressures pı and pz, above the liquid surfaces in the 
tanks (due to pressurization), and the pressure head of the liquid columns 
H, and H: (Fig. 7.49). 

In a stationary rocket, the pressure at the entrance to the turbopump 
assembly for a component designated by subscript “1” will be 


Pt = Pi + yH: 


where y; is the specific weight of the component. 
If the rocket moves with an axial acceleration a, and its axis forms an 
angle y with the horizon (Fig. 7.50), then 


p, = Pit nHila + gar sin ¢)/g 


where g» is the acceleration of gravity at an altitude h. 


CHANGE IN THRUST AND WEIGHT 247 


In this expression, the terms pı, a, ga, Hı, and ¢ are variables along the 
trajectory with time; a increases, ga, Hı, and ¢ diminish, while pı varies 
depending on pressurizing conditions. Sometimes, for instance, pressuriza- 
tion is accomplished by ram pressure. In this case, pressure p, equals the 
atmospheric pressure at a given altitude plus the ram pressure, with a cor- 
rection for compressibility of air. With this system of pressurization, the 
pressure at the entrance to the turbopump first rises and then drops. Analo- 
gous reasoning applies to the second fuel component. 





Fig. 7.50. Development of axial overloading resulting from axial acceleration a. 


Increase in pressure p; will cause approximately the same increase in 
component pressures at the entrance to the combustion chamber, inasmuch 
as the centrifugal pump creates approximately constant pressure drop. The 
picture is the same for the displacement delivery, for which p; is simply 
the pressure at the entrance to the jets. 

The thrust of the motor is a complex function of pressures p: for each 
of the components. During change in the delivery pressure, the total fuel 
consumption varies, as well as the percentage ratio of the components. 

It must be said, however, that the described complex dependence of 
thrust on conditions of delivery results qualitatively in comparatively 
small variations, considerably smaller, in any case, than those which arise 
in the static component of thrust as a result of the changes in the external 
pressure. In long range ballistic rockets the change in atmospheric pressure 
results in thrust variation of approximately 15%, whereas delivery condi- 
tions vary thrust by only 1-2%. 

In speaking of thrust, it is impossible not to mention its variation during 
starting and stopping of the motor. 

Some motors develop their full power not immediately upon fuel igni- 
tion, but after a passage of a more or less prolonged period of time, which 
constitutes some small fraction of the total operating time. Similarly, 
during cut-off of a liquid motor the thrust does not disappear immediately, 
but the so-called aftereffect is observed. After cut-off, there is still some 
thrust, due to the burning of the remains of the fuel. The magnitude of this 
thrust remains somewhat indefinite, which results in a noticeable range 
dispersion for ballistic long range rockets. In order to reduce the indicated 
dispersion, it is possible to shift the operation of the motor to a preliminary 
stage with reduced fuel supply before cut-off. Fig. 7.51 shows the typical 
graph of thrust variation for a motor of a long range ballistic rocket. 
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The weight of the rocket Mg, also changes along the trajectory. This 
occurs because of the change in mass M and the variation in acceleration 
of gravity gn. 
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Fia. 7.51. Variation of thrust with time for a long range ballistic rocket. 


The mass of a rocket at a given time during flight equals the initial 
mass Mo minus the mass of the burned fuel. 


M = My — Í, ' mdt. 
If the consumption remains constant, then 
M = My — mi. 
The acceleration of gravity at an altitude h above sea level is 
ga = gh?/(R + h? 


where g is the acceleration of gravity at sea level, and R is the earth’s 
radius. 


VIH. Rocket Flight Trajectory 


A. Trajectory Phases 


From the instant of initial motion, and ending at the time of impact or 
explosion, the flight trajectory of a rocket is usually divided into two phases: 
the powered phase,fi.e., flight with an operating power plant, and a free, 
coasting phase with'inoperative power plant (Fig. 8.1). Therefore, from the 


y Coasting phase 





Powered phase 


Fig. 8.1. Powered and coasting phases of a rocket trajectory. 


point of view of flight dynamics the difference between trajectory phases is 
that, during the coasting phase, the thrust and the steering forces are 
excluded from the forces acting on the rocket. 

The powered phase in its turn may be divided into phases. 

In those cases in which the rocket is launched from trainable launchers, 
the motion of the launcher is examined separately. This region is primarily 
characteristic of unguided rockets; in particular, short range powder 
rockets. During motion along the launcher, coupling forces caused by 
launcher guide rails are added to the external forces. 

For guided anti-aircraft rockets, the powered phase can be divided, 
depending on the system of guidance, into unguided and guided phases, 
and in the case of a rocket equipped with a seeking or homing system, into 
a homing phase also. 

Long range ballistic missiles are launched vertically from a launching 
pad, and during the first few seconds rise straight up. This phase of the 
flight may be called the launching phase. Further on, the rocket starts its 
entry into the trajectory. The rocket deviates from the vertical and, 
describing an arc in the entry stage, enters the last, inclined phase (cut-off 
phase), where the power plant is switched off. Further on, the coasting 
phase of flight begins (free flight). 
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The coasting phase of the long range ballistic rocket is divided into two 
parts. The first part of the coasting phase takes place in such a rarefied 
region of the atmosphere that the aerodynamic forces have an insignificant 
effect, compared to the force of the rocket weight. During the second part 
of the coasting phase, in moving through denser atmospheric strata, a 
noticeable deceleration and stabilization of the rocket takes place, which 
up to that time moved without a definite angular orientation (Fig. 8.2). 


y ‘Free flight region 





Stabilized region 


Fic. 8.2. Trajectory of a long range ballistic rocket. 


In summarizing the above, it must be noted that the division of the 
trajectory into phases is arbitrary. It is done, primarily, according to the 
presence of various forces determining the motion of a rocket which dictate, 
to a considerable degree, the approach to the integration of differential 
equations of motion. 

Now, let us return to the equations of motion (7.1), (7.2), and (7.3). 
Let us rewrite them, keeping in mind that Eq. (7.2) 1s reduced to the form 
of Eq. (7.16) and that in Eq. (7.3) the quantity M, amy be omitted and 
moment M replaced by the sum Ms: + Ma 


v = (P - D—D,)M— — g sin 6 (8.1) 
ô = (1/Mv)[(P — D,)a + L(1 — e/2,’)] — (g/v) cos 0 (8.2) 
J0 +a + Mat Ma t+ Lz,’ = 0. (8.3) 


Analytical integration of these equations in the general form is impossi- 
ble. Therefore, we will dwell on the simplest specific cases. 
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B. Powered Phase 


1. PowERED PHASE oF A LONG RANGE Bauuistic ROCKET 


It is simplest to determine the law of motion of a long range ballistic 
rocket during the powered phase under the conditions of ideal steering, 
i.e., steering which completely assures conformity to the given regime of 
entry into the free flight trajectory. 

For a given set of coordinates, at the end of powered phase, a definite 
velocity and flight direction must be imparted to the rocket at the termina- 
tion of the powered phase in order for it to reach a given point. This is pro- 
vided, first of all, by programmed steering. 

The rocket, as has been mentioned above, 1s launched vertically. Then 
its axis deviates from the vertical. This deviation takes place gradually, 
and stops shortly before the end of power plant operation. The turning of 
the rocket is performed by steering mechanisms which receive their signal 
from the programming device. The rocket turns relative to a stable vertical 
reference, fixed by a gyroscope, according to a given function of time. The 
angle of inclination of the rocket axis to the horizon g = @ -+ a (see page 
201) is called the programmed angle, and its variation with respect to time, 
the rocket program 


p =0+a = fil). 


The long range ballistic rocket program has the shape of a curve, as 
shown in Fig. 8.3. 

At the end of powered phase, angle ¢ is maintained constant, in order 
to reduce as much as possible the deviation from a given flight direction, 
in case the cut-off of the power plant takes place a bit sooner or later. The 
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Fig. 8.3. The program curve of a V-2 rocket. 


attainment of a given velocity is controlled by a special device, the so-called 
velocity integrator, which at the time of reaching the regimed velocity 
generates a signal to switch off the power plant. 

Let us return to the equations of motion. Suppose that the rocket 
program has been established and that the steering is ideal. 
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For trajectory calculations in the first approximation, it may be assumed 
that the angle a is small, and that the lift force L, because of this, is also 
small compared to the projection of the weight force vector onto the normal 
to the trajectory, Mg cos @. At the same time 


=o = fit) 


i.e., angle @ turns out to be the given function of time. In this case Eq.(8.1) 
can be quite simply integrated numerically, independently of Eqs. (8.2) 
and (8.3). Let us rewrite this equation in the following form: 


Av = [(P — D — D,)M— — g sin 6]At (8.4) 
Two more obvious ones are added to this equation in order to determine 


the coordinates of a moving rocket 


Az = v cos 6At 
Ay = v sin bAt 


(8.5) 


where Az and Ay are displacements of the rocket center of gravity along 
the axes x and y during the time At (Fig. 8.4). The variable coordinates x 


Fıc. 8.4. Calculation of a rocket trajectory. 


and y determine the position of the rocket along the trajectory at any 
instant. 

During integration it is necessary, first of all, to select an integration 
interval Aż, for instance At = 1 sec. The smaller the interval the higher the 
precision of integration will be. 

The integration is done in the following manner. 

We assume that for the initial instant, the rocket velocity is equal to 
zero, and @ = 90°. According to formula (8.4) we find the velocity after a 
time interval At after launching 


o (Py—D, _ 
y= Av = (72 s) At. 
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Here, the drag D is assumed to be zero. Furthermore, it is obvious that 


zx=0 
and 
yi = Ay = vAt. 


Therefore, after time At after launching, the rocket is at an altitude yı. 
In the standard atmosphere table we find the air density p corresponding 
to this altitude; from the graph Cp = f(M), we find the coefficient Cp; 
and from velocity vı, we calculate the drag 


Dı = Cpi X pws 


During time At the thrust P also changes: with increase in altitude, the 
pressure of the surrounding medium is reduced, and the thrust is increased. 
Besides, the thrust may change because of variation in the delivery 
conditions. 

The mass of the rocket will decrease by the amount of fuel expended 
in time Aż, and will have a value M. 

Angle 6, according to the program, may have a value other than 90° 
and will be equal to 61. 

Therefore, at the end of the second step of integration, with the aid 
of formula (8.4) we get 


n= n + (PBR Dot sin a) ai 


and then according to formula (8.5) we find 


Yo = £ı + vı cos 0At 
Yo = Yı + v sin OA. 


Furthermore, we determine parameters at the end of the third step 
P, — D: — D, i 

v + a — g sin a) At 

T3 = Tə —+ Vo COS AAt 

Yz = Yo + v Sin bAt. 


V3 


Similar calculations are performed for all of the given trajectory phase. 

As always, in numerical integration of equations, the results of calcula- 
tions should be tabulated and reduced to a system requiring a minimum 
expenditure of effort. 

The relationship of velocity, rocket acceleration, aerodynamic forces, 
and inertial forces to time and to rocket coordinates x and y along the 
trajectory may be obtained from these calculations. 
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Fig. 8.5 shows the trajectory of the V-2 rocket in a powered phase, and 
Fig. 8.6, variation of velocity, acceleration, and drag with respect to time 
for the same rocket. As can be seen from the curves, the drag initially 
increases sharply, due to increase in rocket velocity. Further on, the drag 
decreases due to reduction of air density p along the trajectory. 

Rocket acceleration increases during the powered phase. This occurs 
because the mass of the rocket decreases and the thrust increases. Only in 
a small region (in a region of velocities close to the ve'ocity of sound) does 
acceleration decrease somewhat, since the drag here increases sharply. 
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Fria. 8.5. V-2 rocket trajectory during Fig. 8.6. Relationships of velocity, ac- 
the powered phase. celeration, and drag with respect to time 
of a V-2 rocket in a powered phase of the 
trajectory. 


The angle of attack a can be determined as a function of time from 
Eq. (8.2), knowing the variation of the angle g. 

Fig. 8.7 shows, as a sample, the change in the angle of attack for a long 
range ballistic rocket. The shape of this curve is almost entirely determined 
by the rocket program. 

Initially, while the rocket, according to the program, rises vertically, 
the angle of attack is equal to zero. With deviation of the rocket from the 
vertical, the angle of attack necessarily assumes a negative value (Fig. 
8.8). If the rocket deviates sharply from the vertical, then the absolute 
value of the angle of attack will be large. Ordinarily, the absolute value of 
a does not exceed 2.5-3°. 

The rocket is programmed in such a way that, at the time when the ram 
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pressure is maximum, the angle of attack would be equal to zero. This is 
necessary for the elimination of transverse aerodynamic forces, under whose 
action the rocket would fail structurally. 

Further on, the angle of attack becomes positive id at the end of the 
powered phase, rapidly increases, inasmuch as the programmed angle 
p = 0 + a remains constant and angle @ decreases, due to the action of the 
acceleration of gravity. 





Fig. 8.7. An example of the relation- Fra. 8.8. Negative angle of attack at 
ship of the angle of attack to time of along the start, and positive angle at the end, of 
range ballistic rocket in a powered phase of the powered phase of the trajectory. 
the trajectory. 


After the angle a has been determined as a function of time, the above- 
described ballistic calculation may be repeated with refined data, consider- 
ing that 


0= p-a. 


The trajectory calculations are not performed for range determination 
alone. In integrating the equations of motion, the forces acting on the rocket 
with respect to time may be determined. This will enable the performance 
of structural analysis of rocket assemblies based on ballistic calculations. 
Specifically, knowing the values of acceleration and drag, it is possible to 
analyze the rocket air frame under the action of axial compressive forces. 
After determining the angle of attack a, transverse aerodynamic forces 
may be calculated, and the structural analysis may be performed to deter- 
mine the significant bending moments on the air frame of the rocket. 


2. THE UNGUIDED, IDEALLY STABILIZED ROCKET 


As a second example of the integration of the equation of motion during 
the powered phase, let us examine the flight of an unguided, ideally stabi- 
lized rocket missile, launched from guide rails. 

We will neglect the moments of inertia which develop during the rota- 
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tion of the rocket about transverse axis. This is equivalent to the assump- 
tions that the stabilizing moment, if the rocket is statically stable (see 
page 237) will always so turn the rocket that its axis will coincide with 
the velocity vector, and the angle of attack a will become zero, and that the 
angular velocity of the velocity vector rotation 6 is sufficiently small. 

In this case, Eq. (8.3) will be satisfied for all points along the trajectory 
since, according to assumption, 6 is vanishingly small, M,, and Mz become 
zero, and force L,, for an unguided rocket, is absent. In Eq. (8.2) the force 
L will become zero, since the angle of attack is equal to zero. 

Equations (8.1) and (8.2) are rewritten in the form 








. P-D ; 
v= a 79a 0 
(8.6) 
ĝ = — : g cos 6. 
V 
These equations are most rationally integrated by the numerical 
method. 
Av = (r — g sin a) At 
1 (8.7) 
A§ = — z9 cos ĝAt. 


Equations (8.5) are added to Eqs. (8.7). 


Ax = v cos bAt 
Ay = v sin bAt. 


During the motion of the rocket along the rails, the value of angle @ is 
constant and is equal to the launcher angle of elevation with respect to the 
horizon @. In this case, the first Eq. (8.7) is integrated independently of 
the second. 

If, while still on the launching rails, the rocket reaches such a velocity 
that the frictional resistance on the rails and drag D become significant, the 
integration of the first Eq. (8.7) should be done numerically. Assuming an 
integration interval At, we find 


P f 
v = Av= (i — g sm a) Al. 


Then, from the obtained velocity vı, we determine drag D, find the new 
value of rocket mass Mı = My — mAt (m is the fuel consumption in a unit 
time), and take the next interval: 


Vo = V, + (252 — g sin a) At. 
1 
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The value of thrust Pı may differ from the value of Po. It is in this manner 
that the graph of change in rocket velocity on the launching rails is cal- 
culated as a function of time. 

If the force D is negligible in comparison with the thrust P, Eqs. (8.6) 


are very simply integrated. 
Assuming the mass discharge m and the value of thrust P as constant, 


we get with the aid of the first of Eqs. (8.6) 
P R 
v= aan g “n 0) A 


v = —(P/m) In (Mo — mt) — gt sin % + C. 


The constant C = (P/m) In Ms, since only in this case is the initial condi- 
tion, v = 0 at ¢ = 0, satisfied. Therefore, 


v = —(P/m) In (1 — mt/Mo) — gt sin bo. 


Therefore, the velocity of the rocket leaving the launcher rails may be 
calculated in one way or another. Let us designate this velocity by vo. 

Subsequently, Eqs. (8.7) are integrated simultaneously. Again, the 
integration interval At is selected; and the velocity vı, which the rocket will 
have in a time interval At after leaving the rails, and the angle @,, are 
determined 


from which 


Py — Do ) 
vv, = 0 —_—— — g sin 4] Af 
1 0 + ( My g 0 
0i = þh — L g COS At. 
Vo 


It is understood that here Po, Do, and Mo are values corresponding to 
the instant of leaving the rails, and not to the initial motion on the rails. 
Furthermore, the values of xı and yı are determined by formulas (8.5) 


zı = To + vo COS HAL 
Yi = Yo + vo Sin At 


and the second interval is performed. 


mt (PP g sin a) as 
1 


1 
= 6, — — g COs 0At 
vi 


Yq = Lı + vı cos 6,AE 
Y2 = Yı + vi sin 0At. 
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The behavior of rocket motion and the shape of its trajectory are deter- 
mined in this manner. 

The indicated integration of equations of motion may also be continued 
through the coasting phase of the trajectory, simply by assuming that at 
the time of motor cut-off P = 0 and M = constant. 

The form of the trajectory of a powder rocket missile is shown in Fig. 8.9. 

The disregard of the inertial moment during turn, to simplify the equa- 
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Fic. 8.9. Trajectory shape of an unguided powder rocket. 
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Fig. 8.10. Motion of a rocket upon leaving the launcher rails. 
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Fig. 8.11. Trajectory of the rocket center of gravity: a, in the absence of angular 
perturbations; b, in the presence of angular perturbations. 


tions of motion, results primarily in the elimination from the obtained solu- 
tion of the oscillatory motions of the rocket. Actually, the initial rotational 
impulse is imparted to the rocket immediately upon leaving the rails (Fig. 
8.10). The rocket rotates in a direction reducing angle 6, and leaves the 
rails with a negative angle of attack. Further on, the rocket is stabilized, 
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performing an angular oscillatory motion. The center of gravity of the 
rocket also receives transverse oscillatory displacements due to the periodic 
change in the lift force L. The trajectory of the rocket center of gravity, 
taking into account angular oscillations, appears as shown in Fig. 8.11. As 
a first approximation, these oscillatory motions may be neglected, which 
has been done above. 


C. Flight Beyond the Limits of the Atmosphere Within the 
Earth's Gravitational Field 


1. EQUATIONS oF MOTION 


The greater part of the trajectory of a long range rocket, or a high 
altitude meteorological rocket, occurs in such rarefied regions of the atmos- 
phere that during this flight phase it 1s possible to neglect aerodynamic 
forces acting on the rocket. If, in addition to this, the flight is in the coasting 
phase, the equations of motion of a rocket, as a mass point, may be inte- 
grated for any rocket elevation above the earth, no matter how great. 


gt? 





Fig. 8.12. Trajectory of a body projected at an angle to the horizon (parabola). 


It is known in physics that any body thrown at an angle with respect 
to the horizontal with velocity vo, moves (if the air resistance forces are not 
considered) along a parabola. Indeed, at time t, after initial motion, the 
coordinates of a thrown body will be 


Xx = vol COS bo 
y = vol sin bo — $a. 


By eliminating t, we find the equation of the parabola shown in Fig. 
8.12: 


g 


m — 2 e—a 
y x tan bo x Sng cos? 8, 


The obtained expression, however, is true only within limited bounda- 
ries. During its derivation it was assumed that the gravity acceleration 
vector g remains parallel to the axis y for all points on the trajectory, and 
does not change in value. Actually, this is not so. Acceleration due to gravity 
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g decreases with increase in altitude, proportionally to the square of the 
distance from the center of the earth. Besides, the force of earth’s gravity is 
constantly directed toward its center, and g vectors at different points on 
the trajectory will not be parallel to each other. These properties of earth’s 
gravitation cannot have significance if we are talking about low altitudes 
and short distances of the thrown body. However, if we are interested in 
long distances and high altitudes peculiar to long range ballistic and many 
meteorological rockets, then the indicated condition must necessarily be 
kept in mind. 

The problem of determining the trajectory of a rocket flight under 
these conditions corresponds to the problem of determining the trajectory 
of an astronomical body in free flight (Kepler’s problem). The theory of 


ds g 
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Fig. 8.13. Derivation of the equation of rocket motion in a polar system of coordinates. 


body motion under these conditions is called “elliptical,” as differentiated 
from the above-mentioned elementary problem of body motion along the 
parabola. The greatest interest of the elliptical theory lies in solving basic 
problems of astronautics, and other problems, such as the determination 
of the trajectories of a rocket earth satellite and of an interplanetary ship 
after its escape from the field of earth’s gravitation. 

Let us examine the motion of a rocket in a polar system of coordinates 
whose origin is at the center of the earth. 

Let a rocket, represented in Fig. 8.13 by point A, move by inertia along 
some trajectory, and at a given time be located at a distance r from the 
center of the earth. During transition from point A to point B, the kinetic 
energy of a rocket Mv?/2 changes by the amount d(Mv?/2). The potential 
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energy will change by the amount d(Mg,r). The acceleration of gravity g 
is not a constant quantity. 

The change in kinetic energy is equal to the change in the potential 
energy, inasmuch as the rocket motor 1s not operating: 


d(4Mv? — Mg,r) = 0. 


Since the mass M remains constant 


łu — g,r = constant. 


The acceleration of gravity 1s inversely proportional to the square of 
the distance to the center of the earth 


g = ght 


where R is the radius of the earth; and g, the acceleration of gravity at sea 
level. 
Therefore 


sv? — gR?/r = constant. 
At the initial point on the trajectory, where r = ro, v = vo. From this 
ły? — gR?/r = $v — gR?/ro. (8.8) 
but the velocity 


4 ds = / dr? + rèd? ae 7. 1. odd 
Mig orgy ee Ae ee 
Then, instead of expression (8.8), we will get 
4(r + rg?) — gR?/r = wè — gR?/ro. (8.9) 


Now let us derive the second equation. We can obtain it from the con- 
dition that the momentum of the rocket, relative to the center of the earth 
during coasting, is constant 


Mrv cos 6 = constant. (8.10) 


Angle @ is shown in Fig. 8.13; v cos 6 is a velocity component perpendicu- 
lar to the radius. 
It follows from expression (8.10) that 


rv COS Ô = TWWo COS Oo. 
but it can be seen in Fig. 8.13 that 
cos 0 = rdy/ds. 
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Therefore 
rv cos 0 = r(ds/dt)(rdg/ds) = rọ 
and 
T? = voro COS Ap. (8.11) 


If Eqs. (8.9) and (8.11) are solved simultaneously, and if time ¢ is elim- 
inated, we get the relationship of r with respect to ¢, i.e., the equation of 
rocket trajectory. 


2. TRAJECTORY OF FLIGHT 
According to Eq. (8.11) 


e _ Volo COS bo 
y2 


but 
o = dg/dt = (de/dr)(dr/dt) = (dẹ/dryr. 
Taking this into consideration, we will get 


NE dr Varo Cos Oo 
de r? 
or 


r = (dr/dẹ)(k/r°) 


k = Voro COS Ao. 


where 


Substituting this expression for r, and using the expression for ¢ obtained 
from (8.11) into Eq. (8.9), then 
(dr/dy)*k?/r4 + k?/r? — 2gk?/r = v? — 2gR?/ro. 


From this 


(dr/dg)k/1? = Vu? — 2gR?/ro — k?/r + 2gR?/r 
or 
—d(k/r) l 
V vè — 2gR?/ro — k?/ + 2gR?/r 
Under the square root sign let us add and subtract the constant term 
g?R4/k?, and also introduce the term —gR?/k under the differential sign in 
the numerator, after which we will get 
ma —d(k/r — gR*/k) | 
V (wè — 2gR?/7 + g?R4/k?) — (k/r — gR?/k)? 


de = 
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Integration of the last expression yields 
k/r — gh?/k 
V 0c? — 2gR?/ro + GP RA/K 
from which the required trajectory equation is obtained 
Pea k?/gR? l 
1 + [(k/gR?) Vv — 2gR?/ro + gRt/k?)] cos (p — go) 


p — po = arc cos 


(8.12) 


The constant of integration go depends upon the reference angle ¢ (see 
Fig. 8.13). 
Let us introduce the following notation 


k? Vy" COS” Oo 


P= eg (8.13) 
-E fpe 2E , oR! 
—— gR? ajo To + k? 
or 
- Opi pa A pira one? D. 
e= J 1 — “ite cos ae ae (8.14) 


The final equation of the flight trajectory (8.12) will take on the form 
r = pil + e cos (p — p). (8.15) 


Let us analyze the obtained expression. 

It is known from a course in analytical geometry that Eq. (8.15) is an 
equation of the second order in a polar system of coordinates with the origin 
at one of the foci of the curve. Coefficient e is the eccentricity of the curve. 
When e < 1, Eq. (8.15) represents the equation of an ellipse; when e > 1, 
the equation of a hyperbola; and, finally, when e = 1, the equation of a 
parabola. 

The parabola obtained from Eq. (8.15), when e = 1, should not be con- 
fused with that which was obtained earlier for a freely projected body 
moving in the earth’s gravitational field under the influence of a gravita- 
tional vector, whose absolute value and direction are constant. 

First, let us set e = 0. Equation (8.15) will take the form 


r = p = constant. 


In this case the body will move in a circle. 
From Eq. (8.14) we will get 


ve? = (gR?2/r)(1 & V1 — (1/cos? 6). 
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It is not difficult to see that the velocity v determined by the last 
formula has significance only in a case where cos? 6) = 1 and ĝo = 0. This 
means that the motion of a body along the circumference is possible only 
when the direction of the initial velocity vector is parallel to the horizon. 

We get an expression for flight velocity 


vo = V/ gR?/r9 = Vj. 


This is an expression for the orbital velocity with which we are ac- 
quainted from Chapter I. When ro = R, the numerical value of the orbital 
velocity: 

vı = 26,000 ft/sec. 


Now, let us determine under what conditions the eccentricity e becomes 
greater than unity. 

From expression (8.14) it can be established, without difficulty, that 
independent of the elevation angle @ of the initial velocity vector, when 


Vo > ~ 2g R?/ro 


the value 
e >l. 


In this case the trajectory of motion will be a hyperbola. 
With 
Vo = / 2gR?/ To 
the motion is along a parabola, and with 
Vo < «v/V 2gR?/ To 


along an ellipse. 
The velocity 


ve = / 2gR?/r9 


is called the parabolic velocity or, sometimes, the escape velocity. 
Assuming that ro = R we will get for the escape velocity the following 
numerical value: 


v = V2gR = 37,000 ft/sec. 


Figure 8.14 shows the variation in the shape of the trajectory of a body 
with respect to velocity at a constant projection angle 0) = 30°. 

At low velocity, a body released at an angle 6) with respect to the horizon 
will return to earth, having described an elliptical arc. The range of flight 
will be extended with an increase in initial velocity. The are described by 
the flying body will also increase. 
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The trajectory becomes a parabola at velocity vo = v2. A body projected 
from the earth with velocity vo > v will not return to the earth. Hence, the 
name for v—escape velocity. | 

Up to the present time, such high velocities were not attained. However, 
with the velocities of the long range ballistic rockets which have been 
attained, the trajectory calculations according to the elementary ‘‘para- 





Fig. 8.14. Trajectories of a body related to the projection velocity at 0 = 30°. 


bolic” theory cannot yield sufficiently accurate results and it is necessary 
to perform the calculations according to the formulas of the “elliptical” 
theory. 

The following procedure may be followed. 

Up to a certain altitude, while the effect of the aerodynamic forces is 
still significant, the trajectory calculation may be carried out by the 
method of numerical integration described above. Starting with the time 
when it becomes possible to neglect drag, the calculation should be per- 
formed according to formulas (8.13), (8.14), and (8.15). In this way, it is 
possible to get the full picture of the variations in velocity and acceleration 
for the entire trajectory, depending on the duration of flight. 

Figure 8.15 shows the graph of velocity and tangential acceleration of 
the ballistic rocket for its entire flight trajectory. 
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The velocity increases up to the time of motor cut-off. After cutoff, the 
velocity decreases, while the rocket gains altitude due to its momentum. 
After the rocket passes the trajectory apex, the velocity again starts to 
increase. Afterward, due to the action of air resistance forces, the velocity 
decreases. 

Initially, acceleration increases in the powered phase, shown in Fig. 8.6. 
During the switch over of the motor to the preliminary stage, the accelera- 
tion sharply decreases. Following this, at full cut-off of the motor, the accel- 


V ft/sec V ft/sec? 
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Fig. 8.15. Velocity and tangential acceleration along the trajectory of a ballistic 
rocket. 


eration becomes negative. Obviously it equals —g sin 6. Further on, the 
tangential acceleration varies as the result of variation in @ and g. At the 
apex of the trajectory it equals zero. During the descending phase of the 
trajectory the force of earth’s gravity accelerates the motion of the rocket. 
At reentry into the upper regions of the atmosphere, there begins a retarda- 
tion of the rocket due to air resistance. At the very end of the trajectory 
the rocket velocity drops noticeably and the drag force lessens. The abso- 
lute value of the acceleration also falls at that time. 


3. RANGE 


At low velocities the range of a body projected at an angle 4 to the 
horizon is determined by an elementary formula derived from expressions 
on page 259 

L 


(v/g) sin 20o. 


The maximum range is at 6) = 45°. This result is well known. 
Things are more complicated at high projection velocities. 
In this case, it is more expedient to express the range of the rocket by 
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the angle 28, formed by the polar radius which passes through points A 
and C, the initial and final points, respectively, on the trajectory (Fig. 
8.16). Due to the symmetry of the trajectory, it is sufficient to determine 
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Fic. 8.16. Range determination of a ballistic rocket. 


half of this angle—angle 8 (from point A to point B on the trajectory). 
Let us find angle po in Eq. (8.15), setting ọ = 0 and r = 7% 


p — To 
yo = — arc Cos s 
ero 





Equation (8.15) now takes the form 


4—1 
r= p | 1 + e cos (o + are cos = r) | . 


0 





Let us consider the eccentricity e < 1, which corresponds to the present 
velocity capabilities. 


When the expression contained in brackets in the last equation 
becomes equal to r, the radius r will, obviously, reach its greatest value 


Tmax = p/(1 — e). 
At the same time, angle ¢ will be equal to 8. Therefore, 
B + arc cos [(p — 1o)/ero] = r. 
Substituting expressions (8.13) and (8.14) for p and e we get 
oe (vo?ro/gR?) cos? bo — 1 
E et GOR [1 — (2vro/ gR?) cos? 60 + (vor?/g R) cos? 0o)!” 


If this expression is differentiated with respect to cos? bo, and the deriva- 
tive is equated to zero, then, without undue difficulty, it is possible to find 
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that angle of projection @ at which the value 8 will be maximum for a given 
velocity vo 


cos? @ = a PY 5 = 7k? 
or 
sais I 
Cos’ Oo = 2 — v¢/ve (8.16) 


where v; is the orbital velocity. 
If velocity vo is low 


cos? Oy = 4. 


This result has already been shown above; in order to obtain maximum 
range, the vector of the initial velocity must have an angle of 45° with the 
horizon. 

When v = v, as has been noted, cos? — = 1 and 6 = 0. This case 
corresponds to the constant rotation of a body about the earth with an 
orbital velocity. 

Figure 8.17 shows a graph of the optimum projection angle as a function 
of the ratio v/v. It is clear from the curve that the optimum angle de- 
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Fic. 8.17. Relationship of the optimum projection angle 6 to the initial velocity v 
(vı is the orbital velocity). 


creases with an increase in the projection velocity. At velocities of the order 
of 6500 ft/sec, the angle 4 only differs slightly from 45° (it is equal to 
approximately 44°). At high velocities it diminishes more and more 
noticeably. 

At an optimum angle 6 


—2 4/1 — vo?/ 01" 


ea (8.17) 


B = T — arc cos 


If the square of the ratio v/v, is small compared to unity, then it is 
possible to obtain 


28 a v/ve. 
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For instance, at vo = 6500 ft/sec, 


6500 \) 1 
ie (E) = 6 
The range along the arc AC (see Fig. 8.16) 


2Bro = 26R = z 4000 = 250 mi. 

The presented range calculation is applied to the region of the trajec- 
tory which lies above the dense atmospheric layers (region ABC in Fig. 
8.16). The actual range will be somewhat greater due to trajectory sectors 
AA’ and CC’. The first of these (AA’) is determined, as has already been 
mentioned, by the programming of the rocket into its basic trajectory. The 
region CC’ may be considered, because of its shape, as the continuation of 
the basic elliptic curve. Here, the shape of the trajectory cannot noticeably 
vary from an ellipse. Only the rocket velocity changes noticeably, due to 
air resistance. The error in range determination due to substitution of an 
elliptical arc for region CC” is small, since the region CC” itself is small 
compared to total range. Therefore, the additional range in region CC’ may 
be approximately calculated by formulas (8.13), (8.14), and (8.15). 

The existence of section CC’ affects the optimum angle of projection 8o. 
The condition of maximum range from point A to point C” will differ from 
condition (8.16) postulated for the maximum range from A to C. Due to 
the difference in altitudes of points A and C’, the optimum angle @ is 
actually less than the one determined by formula (8.16). For instance, for 
a long range ballistic rocket, traveling at the end of powered phase with 
a velocity of 5000 ft/sec, the optimum projection angle bo, calculated accord- 
ing to formula (8.16), is equal to 44°30’. Actually, due to the existence of 
section CC’, this angle must be approximately 42°. 


4. Some Basic PROBLEMS IN MASTERING INTERPLANETARY SPACE 


In concluding this chapter, devoted specifically to the problems of 
flight with velocities commensurate with interplanetary velocities, let us 
pause briefly on the means of mastering interplanetary space. First of all 
we must agree which mastery we mean: rocket flights into interplanetary 
space with people as passengers, or without them. The first problem is 
immeasurably more complex than the second and cannot be the order of the 
day. At the same time, rocket flights into interplanetary space without 
passengers seem to be entirely possible in the immediate future. 

The first, basic, and still unsolved, problem in this direction is the con- 
struction of an artificial earth satellite in the form of a rocket, equipped with 
transmitting instruments. The solution of this problem would mean, first 
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of all, an attainment of a definite and quite high degree in the development 
of rocket technology. At the same time, with the aid of an artificial satellite, 
it would be possible to collect extensive scientific data which would enable 
us to check and refine our concepts of the physical conditions beyond the 
limits of the earth’s atmosphere. 





Fig. 8.18. Direction of a rational rocket path for launching an artificial earth satellite. 


In order to put an artificial rocket satellite beyond the limits of the 
atmosphere it is necessary to impart to it a velocity of the order of 26,000 
ft/sec. The velocity of the earth’s rotation may be partially used for that. 
At the equator the earth’s surface has a circumferential velocity of about 
1500 ft/sec, directed from west to east. Therefore, in putting an artificial 
satellite into orbit, it is more reasonable to launch the rocket from the 
equator in an easterly direction (Fig. 8.18). 





Fic. 8.19. Symmetrical trajectory around the moon. 


The second interesting problem in mastering interplanetary space is 
the launching of a rocket to circle the moon, with return to earth (Fig. 
8.19). With this flight it would be possible to obtain a photograph of the 
part of the moon invisible from the earth. In order to solve this problem, 
attainment of a considerably greater velocity is necessary; specifically, 
velocity close to escape velocity. But this is not the only difficulty. Calcula- 
tions show that for a rocket to return to earth after flying around the moon 
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it is necessary to maintain the magnitude and direction of rocket velocity 
very strictly at the end of powered phase. 

With a symmetrical trajectory shape as shown in Fig. 8.19, it is neces- 
sary that the angle 6) be within limits of from 81° to 89°. It is not difficult 
to satisfy this condition. At the same time it is necessary that the velocity 
vo be within limits of from 0.99 v to 1.0001 v2, where v2 is the escape velocity. 
If velocity vo is less than 0.99 ve, the rocket will fall on the moon. Should the 
velocity be greater than 1.0001 v: the rocket, having flown around the 
moon, will not return to the earth. To maintain such a narrow velocity 
margin represents quite a difficult problem at the present time. It seems 





Fig. 8.20. Trajectory of a rocket flight to the planet Mars with a minimum expendi- 
ture of energy. 


that for its solution it may not be necessary to increase the precision of 
steering mechanisms, but to introduce auxiliary motors which, when 
required, will correct the rocket flight along a trajectory. 

In launching a rocket to the moon, just as in launching a satellite, the 
rotational velocity of the earth may be partially exploited. 

Problems even more complicated than the two described are involved in 
reaching the nearest planets, Mars and Venus. 

If it were possible to place a rocket on these planets, even without 
passengers but equipped with powerful transmitting apparatus, we could 
expect to obtain considerable data on the physical conditions existing on 
the surface of these planets. 

In launching a rocket to Mars or Venus along an optimum trajectory, 
not only the earth’s rotation about its axis may be used, as was shown 
above, but the earth’s velocity around the sun as well. The optimum flight 
trajectory to Mars from the standpoint of energy expenditure (but not of 
time) is shown in Fig. 8.20. It represents a curve, approaching an ellipse, 
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with the sun at one of its foci. For free flight along such a trajectory, the 
rocket, upon leaving the earth, must have a velocity of about 8.7 mi/sec 
imparted to it. 

The problem of mastering interplanetary space represents a very com- 
plex but not an insoluable problem. In 1929 C. E. Tsiolkovski wrote on 
this subject: 

“Interplanetary travel cannot even be compared with flying through 
the air. The last is a toy compared to the first. 

Undoubtedly, success will come, but the time necessary for this achieve- 
ment cannot be predicted by me. 

The notion of the ease of its solution is only a temporary fallacy. Of 
course, it is beneficial since it adds courage and strength. 

If the difficulties of the problem were known, many who are now work- 
ing with enthusiasm would be dismayed. 

But, how beautiful the achieved goal will be.” 

The time that has passed since these lines were written shows that 
C. E. Tsiolkovski was very right in his judgements. This problem is very 
complex indeed. At the present time, rocket technology is on the threshold 
of solving the first problem of interplanetary flight—the creation of the 
artificial earth satellite. 


IX. The Basic Principles of Stabilization and Steering 


A. Methods of Stabilizing and Steering a Rocket 


1. STABILITY AND STABILIZATION 


The motion of a rocket under the influence of regular, systematically 
acting forces (regardless of whether or not these forces are considered in 
the equations of motions) is customarily called the undisturbed motion of 
a rocket. For instance, the motion of a long range ballistic rocket, subjected 
to regular thrust and aerodynamic forces, will be considered undisturbed. 
The motion of an anti-aircraft rocket seeking a maneuvering airplane, in 
spite of the indeterminate nature of this maneuvering, must also be con- 
sidered as undisturbed. Along the pursuit trajectory the rocket is acted 
upon by steering forces, which are subject to the prescribed conditions, 
depending upon the method of guidance. 

Aside from the regular forces, random, short-duration, external influ- 
ences may act on the rocket in flight, such as wind gusts, short-duration 
small variation in the motor thust, etc. Under the action of these forces the 
rocket will take on a motion which is called the disturbed motion. 

Let us assume that the action of random forces is small, and that the 
disturbed motion differs only slightly from the undisturbed. If, after the 
discontinuance of the action of the disturbing forces, the oscillatory motion 
of the rocket again becomes nonoscillatory, this nonoscillatory motion of 
the rocket should be called stable. If, however, after the discontinuance of 
the action of the external disturbance the motion of the rocket does not 
approach nonoscillatory motion, then such motion will be considered 
unstable (Fig. 9.1).* The concept of stability and instability of the motion 
is applicable, naturally, not only to rockets, but, in general, to the motion 
of any body. 

It is understood that in all cases, if we wish to control motion, it is neces- 
sary, first of all, to eliminate undesirable influences of random external 
disturbances and thus assure the stability of the motion. 

In actual cases, the rocket in flight is continually subjected to the action 
of various random disturbing factors. In the time interval between two 
consecutive external disturbances, the oscillatory motion of a rocket does 


* In practical investigations of rocket stability, the motion is considered to be stable 
even in the case where the rocket, after discontinuance of action of the disturbance, 
returns not to the nonoscillatory, but to some other, slightly different motion. 
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not have time to return to the parameters of the initial nonoscillatory 
motion. Therefore, under actual conditions of flight, the motion of the 
rocket only approximately corresponds to that required. The rocket con- 
tinually follows an oscillatory motion whose parameters fluctuate about 
the parameters of the stable nonoscillatory motion. If the deviations from 
the desired motion are not large, then, for practical purposes, the motion 
of the rocket may be considered as being stable. - 

In connection with what has been said, the above definition of stability 
of motion may be changed. If, in the final analysis, the effect of the dis- 
turbing factors on the rocket varies, within allowable limits, the motion 
may be considered stable. Or, briefly put, in stable motion the external 
random disturbances do not cause large variations in motion parameters. 


Trajectory of 


stable undisturbed motion Trajectory of unstable motion 








Trajectory of of denied 
disturbed motion motion 


Fig. 9.1. a, Stable; and b, unstable rocket motions. 


The boundaries of allowable variations in the prescribed motion of a 
rocket are determined by a whole series of considerations. Major among 
them is the requirement for accuracy of aiming at a target and the struc- 
tural integrity of the rocket since, at large angles of attack and transverse 
accelerations, the rocket is subjected to considerable side forces. 

It must be said that motion stable relative to one set of parameters 
may be unstable relative to another set. Therefore, when the stability of 
motion is discussed, especially stability along allowable deviations, it is 
necessary to specify the type of motion which is being considered. The 
motion of a rocket is defined by the motion of its center of gravity and the 
angular displacement of the whole rocket as a rigid body about the center 
of gravity. 

The concept of stability and instability, in the sense of angular displace- 
ment, is usually replaced by the concept of angular stabilization or unstabil- 
ization of a rocket. If a rocket to which angular deviation has been imparted 
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returns to a nonoscillatory motion after the disturbing influence ceases to 
act, the rocket is called stabilized. If this does not happen, then the rocket 
is called unstabilized. 

In many cases the questions of angular stabilization may be examined 
independently of the stability of rocket motion as a mass point, even 
though, in the example of an unguided powder rocket (see page 258), we 
saw that the angular disturbance is related to the disturbances in the mo- 
tion of the center of gravity. The rocket may be found to be unstabilized 
even when deviation from the prescribed motion of its center of gravity 
may not fall outside the allowable limits. 


2. PROPERTIES OF OSCILLATING MOTION AND STABILIZATION 
oF UNGUIDED ROCKETS 


Stabilizing aerodynamic moment and aerodynamic stabilization have 
already been mentioned in Chapter VII. If, in referring all aerodynamic 
forces to the rocket center of gravity, the stabilizing moment acts to de- 
crease the angle of attack, the rocket is called statically stable. If the mo- 
ment is directed in an opposite direction, the rocket will be statically 
unstable.* 

In Chapter VII the concept of the center of pressure as an intersection 
point of the resultant of aerodynamic forces with the axis of the rocket was 
also developed. In the case in which the center or pressure 1s located aft of 
the center of gravity the rocket is statically stable. If the center of pressure 
is located forward of the center of gravity the rocket is statically unstable. 

The ratio of the distance between the center of pressure and the center 
of gravity to the length of the rocket is called the margin of stability. 

With a positive margin of stability the rocket having an angular devia- 
tion is turned by aerodynamic forces to its initial position. When the angle 
of attack becomes zero the restoring force will become zero. But, since the 
rocket maintains angular velocity, the motion will continue in an opposite 
direction. An angle of attack with an opposite sign will develop and angular 
oscillations of the rocket will arise. Due to the presence of damping forces 
the oscillations will decay (Fig. 9.2). 

Oscillations of this type are called short period oscillations. Their fre- 
quency depends on the margin of stability and the moment of inertia of a 
rocket. The higher the margin of stability and the lower the moment of 
inertia, the higher the frequency of these oscillations, the greater the damp- 


* In the aircraft industry the frequently used expression ‘‘static stability” is equiva- 
lent to the concept of “aerodynamic stability.” Static stability is ordinarily considered 
to be stability without interference of steering devices; as opposed to dynamic stability, 
stability of motion insured by the introduction of a steering system, it being immaterial 
whether the steering system is automatic or nonautomatic. 
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ing, the quicker the oscillations will decay, and the faster the rocket is 
restored to the steady state motion. 

However, excessive increase in the margin of static stability may yield 
negative results. A guided rocket with a high margin of stability becomes 
unresponsive. The unguided rocket with a high margin of stability will give 
large values of dispersion (see page 318). 


49 


Ag t 


Fig. 9.2. Decay of a short period oscillation. 


The short period angular oscillations described above may occur, as is 
easily understood, about the three central axes. 

The presence of the so-called long period, or phugoid oscillations, is 
peculiar to a rocket, especially a winged rocket, of the pilotless aircraft 
type. Their origin is most conveniently explained by examining a uniform 
horizontal flight of a pilotless aircraft whose thrust, on the average, is 
balanced by a drag force at a constant angle of attack. 


soes Ta 


Fra. 9.3. Phugoid, or long period, oscillations of a pilotless aircraft. 


Let us assume that, for some reason, the thrust has increased for a 
short period of time. In this case the flight velocity increases, and lift force 
also increases, becoming greater than the weight force. The pilotless air- 
craft starts to climb and because of this loses speed. There comes a time 
when the lift force becomes equal to the weight force, but the aircraft, 
having vertical velocity, continued to climb because of momentum. Con- 
tinuing climb will result in the reduction of the lift force. The aircraft will 
begin to descend. The phugoid oscillations develop in this manner. The air- 
craft trajectory for this case is shown in Fig. 9.3. 

The phugoid oscillations have a long period, frequently considerably 
exceeding the operating time of the rocket motor, i.e., exceeding the 
time of flight in the powered phase. Therefore, phugoid oscillations need 
not be considered for many types of rockets. 

The various kinds of rocket oscillations are not independent. During 
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transverse angular oscillations, due to periodic development of lift force, 
transverse oscillation of the center of gravity appears about the steady 
state trajectory. 

If the rocket moves with an angle of attack, then, with the appearance 
of the yaw angle, the rocket will roll. This is clearly evident in Fig. 9.4, 
which shows a rocket with a yaw angle y. The right stabilizer has addi- 
tional velocity, the lift force on it increasing, whereas on the left it decreases. 
A moment develops, rotating the rocket about its longitudinal axis and 
giving rise to rocket roll. 

The problem of stabilization consists of limiting and damping out as 
quickly as possible all oscillations of the types described. 

With relatively simple static stabilization, the effects of short duration, 
random disturbing influences acting on the rocket during flight are satis- 





Frc. 9.4. Development of roll during rocket yaw. 


factorily excluded, but the influence of systematic, constantly acting fac- 
tors which are not considered beforehand remains. To their number belong, 
first of all, technological errors. 

The rocket cannot be fabricated exactly. The rocket motor always has 
a certain eccentricity, and the line of thrust does not pass through the center 
of gravity of the rocket. This creates a moment, turning the rocket side- 
ways. The external contours of the rocket cannot be fabricated ideally. 
Therefore, aerodynamic forces also have a certain eccentricity. Stabilizers 
themselves, whose purpose it is to insure stability of motion, always have 
inaccuracies due to fabrication which lead to the deviation of the rocket 
from the required trajectory. The aerodynamic stabilization cannot be 
effective during flight beyond the limits of the atmosphere. 

There is another possible way of assuring stability of motion of an un- 
guided rocket. It consists of imparting to the rocket rapid rotation about 
its longitudinal axis. This method, as is known, is widely used for stabilizing 
artillery shells. In this case, the rocket maintains the direction of its axis 
due to the gyroscopic effect. 

The method of spin stabilization is used primarily in small powder 
field rockets. These rockets, as we already know, are called turbo-reactive 
shells. Large rockets cannot be stabilized by these means. With rapid rota- 
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tion it is difficult to insure the necessary strength of a large rocket, let alone 
the difficulties in making such a rocket guided. 

A more refined, but immeasurably more complex, method than the two 
described is the installation in the rocket of an automatic stabilizer, which 
will at the same time guide the motion of the rocket according to a pre- 
determined program. 


3. STABILIZATION BY MEANS OF AN AUTOMATIC PILOT 


Automatic stabilizers were first produced to insure the constant direc- 
tion of motion of marine torpedoes. Later, with the development of aviation 
and a need for automatic steering of aircraft, autopilots were produced 
whose basic principles were applied to automatic stabilizing apparatus for 
rockets and pilotless aircraft during the Second World War. 

The basic operating principle of such robots is most easily understood 
from an example of the simplest steering device (Fig. 9.5). 
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Fia. 9.5. Schematic of the simplest automatic pilot. 


The actuating member of the simplest steering device is a cylinder with 
a movable piston which actuates the vane of the aircraft. The pressure on 
both sides of the piston in the cylinder is regulated by a distributing valve 
which is actuated by a sensitive diaphragm. In order to supply oil to the 
right or left cylinder cavities the piston valve must be moved correspond- 
ingly to the right or to the left. The displacement of the regulating valve 
does not require appreciable force and, therefore, the thin sensitive dia- 
phragm is able to control the motion of the rudder. 
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The pressure is admitted into both chambers of the diaphragm housing 
from the pneumatic system of the aircraft through ports half closed by 
slides. A yoke, rigidly attached to the gyroscope, and therefore remaining 
motionless during the turning of the aircraft, serves as the slides. While 
following the course, both ports are closed to the same degree and the pres- 
sure in both chambers of the diaphragm housing is equal. The rudder 
remains in neutral position. 

Let us suppose now that the aircraft, for some reason, begins to deviate 
from the established course, for instance, to the left. Port A will be par- 
tially closed and port B opened. The pressure in the right chamber of the 
diaphragm housing will rise and the piston of the valve will be displaced 
to the left. This will result in admission of oil into the forward chamber of 
the actuating mechanism (steering device). This will turn the vane counter- 
clockwise and the aircraft will begin a turn to the right, 1.e., resume its 
course. 

It would seem that the described robot fully solves the problem of steer- 
ing. However, this system has a major deficiency. 

Let us see how the aircraft will move later on. 

After the aircraft has assumed the desired direction the slides, dia- 
phragm, and valve will return to neutral position, but the piston of the 
actuating mechanism and the vane will remain in the same position and the 
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Fie. 9.6. Motion of an aircraft while controlled by the simplest automatic pilot. 


aircraft will consequently continue in a right turn. Only when the aircraft 
overrides to the right, the slides close port B and open port A, and the valve 
moves in the opposite direction, will the vane force the aircraft to turn left. 
This gives rise to oscillations about a course—yawing. The aircraft will 
travel in a “snake” path (Fig. 9.6). 

A similar pattern could be observed if the driver of a car in turning to 
the right began to straighten out his steering wheel only after the car had 
turned fully to the right. If the driver continued to repeat this error, i.e., 
was late in straightening out the wheel to a neutral position, the car would 
also travel like a snake. In order to eliminate this motion the driver must 
begin straightening out the wheel prior to the time when the car has 
assumed the proper direction. 

Obviously, it is necessary that the steering device, in controlling the 
motion of the aircraft, return the vane from its deflected position at the 
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time the aircraft is back on course. This is the purpose of the feedback. In 
the broad sense of the word, the feedback is an element of a regulating sys- 
tem tying in the following link with one of the preceding ones in such a way 
that the signals transmitted to the regulating system conform with the 
results ‘performed by the system. 

Figure 9.7 shows the same steering mechanism, but with the vane feed- 
back. Here the air distributor AB is no longer rigidly connected with the 
aircraft, and is capable of rotation about a vertical axis. This rotation is 
linked to vane movement. 





Fic. 9.7. Schematic of an automatic pilot with vane feedback. 


Let us now see how the aircraft will get back on course. 

Again, let us assume that for some reason the aircraft deviated from 
the course to the left. The elements of the mechanism will then operate in 
the above-described sequence and the vane will turn counterclockwise. 
Simultaneously, by means of feedback, the vane will turn the air distributor 
and will return it to the neutral position. While the aircraft is turning to 
the right, port A will begin to open and port B to close, with the result 
that the vane will turn in the opposite direction. When the aircraft resumes 
its course the vane and, with it, the slide of the pressure distributor, will 
assume a neutral position relative to ports A and B. The resumption of 
course by the aireraft will take place without yawing (Fig. 9.8). 

The above-described oscillating instability (see Fig. 9.6) may be avoided 
by methods other than the introduction of a rigid feedback from the rudder. 
It is possible, for instance, to give the actuating device—the steering 
mechanism—signals not only proportional to the angle of rocket deviation 
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from the true direction, but also proportional to the velocity and accelera- 
tion of this deviation, i.e., by introducing the so-called proportional feed- 
back. In this case, during rapidly increasing rocket deviation the stabilizing 
devices react to the deviation more actively and return the rocket on course 
more rapidly. This type of stabilizing system is used in ballistic rockets. 

The operation of the automatic steering in yaw was described above. 
It is understandable that for a pilotless missile or a rocket this automation 
must be realized not only in yaw but also in pitch and roll. 

The function of automatic stabilization in a rocket is not confined to 
pure stabilization alone. Inasmuch as the stabilizing mechanism controls 
the motion of the rocket according to a definite program, it also becomes a 
steering device. 

The rocket equipped with a stabilizing mechanism is called a guided 
rocket. If its guidance, from the moment of launch, is not tied in to the 
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Fic. 9.8. Motion of an aircraft while controlled by automatic pilot with feedback. 


earth, then such a guidance system is called self-contained. However, if the 
link with the ground is maintained after launch, for instance, by a radar 
beam, then the guidance system is called remote. 

The self-contained system has one obvious deficiency, namely, after 
the moment of launch any further interference with the motion of the rocket 
is precluded. In particular, it is impossible to correct the motion of the 
rocket if it begins to deviate from its course. 

Transmission of a signal along a beam enables introduction of correc- 
tions into the motion of a rocket and attainment of a considerably greater 
accuracy of hits. For anti-aircraft rockets beam guidance allows the aiming 
of the rocket directly at the enemy aircraft and sharply increases the per- 
centage of hits as compared to unguided rockets. 

However, beam guidance has a major drawback in that the enemy may, 
by means of creating sufficiently powerful jamming, nullify the advantage 
of remote guidance. 

The problem of beam guidance is one of the basic problems of contem- 
porary rocket technology. Concurrent with the development of beam guid- 
ance, means of jamming as countermeasures are also developed. The devel- 
opment of the last necessitates the development of countermeasure-proof 
guidance apparatus. All of these questions are problems of a special branch 
of radio technology (electronics). 
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B. The Gyroscope and Its Application 


1. PROPERTIES OF THE GYROSCOPE 


A massive, precisely balanced flywheel, rotating with high angular 
velocity, is called a gyroscope. 

Everyone is well aware of the peculiar behavior of the gyroscope; the 
rapidly rotating flywheel has the ability more or less to maintain the 
direction of its axis stably, and displays an unexpected “disobedience” 
when an attempt is made to rotate this axis. 

The gyroscope is the basic element of contemporary autopilots and the 
majority of stabilizing devices. Therefore, it is reasonable to begin the exam- 
ination of the operation of stabilizing devices with the gyroscope. 

The theory of the gyroscope is a specific problem in the general theory 
of a moving body with one stationary point. Equations of motion may be 
derived for such a body (Euler’s equations) which in a general form cannot 





Fig. 9.9. Gyroscope in a gimbal suspen- Fic. 9.10. Orientation of fixed gyro- 
sion. scope axes. 


be solved, but which enable us to obtain answers to certain specific prob- 
lems of motion, including the problem of motion of a symmetrical, rapidly 
rotating disk. The last problem constitutes the substance of the so-called 
theory of the gyroscope. We will not delve into this theory and will restrict 
ourselves primarily to the examination of its qualitative side. 

First of all, we will assume that the gyroscope consists of a symmetrical 
body with one stationary point, and that this point is its center of gravity. 
A suitable mount may be made using a Cardan suspension (Fig. 9.9). Under 
these conditions, the motion of the gyroscope may be represented as the 
rotation, with angular velocity Q, about some instant axis passing through 
the center of gravity. 

We} Let us take three mutually perpendicular axes, x, y, z, rigidly connected 
to the gyroscope (Fig. 9.10). Let us direct axis z along the main spin axis 
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of the gyroscope. The other two axes will pass through the center of gravity 
of the gyroscope perpendicular to axis z. 

The vector of the instant angular velocity Q may be resolved along 
these axes into components wz, @,, @,. For a rapidly rotating gyroscope the 
absolute value w: is immeasurably greater than the values of w: and wy. 
These last for a balanced gyroscope are generally equal to zero and appear, 
as we shall see later, only when external forces begin to act on the gyroscope. 
Inasmuch as w: and w are small, w: = Q. 

The components of the main angular momentum vector N will be Aoz, 
Boy, Coz, where A, B, C are the moments of inertia of the gyroscope about 
axes x, y, z. From the conditions of symmetry A = B. 

An obvious fact must be noted: vector N, generally speaking, does not 
coincide in direction with vector Q. Only because, for a gyroscope, wz and 
wy are very small compared to wz is it possible to consider that vectors N 
and Q are directionally quite close to each other. This considerably facili- 
tates the problem and enables the immediate determination of the direc- 
tion of the vector Q by observing the behavior of vector N (being coincident 
with it). 

The basic relationship determining the rotational motion of a body is 
the law of the change of the angular momentum 


dN/di = M (9.1) 


where M is the vector of the moment acting on the body; and N is the vector 
of the angular momentum. 

Let us introduce a system of semi-fixed central axis xv’y’ (see Fig. 9.10). 
These axes are rigidly fixed to the axis z and turn with it, since they are not 
fixed to the rotor of the gyroscope, i.e., they do not rotate with the spin 
velocity Q. For the gyroscope shown in Fig. 9.9 axes z’ and y’ coincide with 
the gimbal axes. 

Let us now see what motion will be performed by the gyroscope axis z 
if we try to turn the gyroscope by applying to it, through the gimbal ring, 
a moment M about axis x’ (Fig. 9.11). 

Up to the time of the application of moment M, the gyroscope has an 
angular momentum 


N = CQ. 


According to expression (9.1) the change in the angular momentum in 
time Af is 


AN = MAt. 


Vector AN coincides in direction with vector M. Adding N and AN we 
obtain a new angular momentum vector N’. Due to the smallness of AN 
vectors N and N’ differ from each other only in direction. Therefore, we 
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see that as a result of the action of moment M, the angular momentum 
vector N, without changing its magnitude, rotates through an angle AN/N 
about axis y’. 

It has already been mentioned above that, for a gyroscope, the angular 
momentum vector almost coincides with the instant axis of rotation. Con- 
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Fig. 9.11. Derivation of an expression for angular precession velocity. 


sequently, we can say that under the influence of moment M the axis of 
gyroscope rotation turns through the same angle AN/N = MAt/N about 
axis y’ during time At. Therefore, the angular velocity of the gyroscope axis 


w = M/N = M/CQ. (9.2) 


The described motion of the gyroscope about axis y’ is called precession, 
and angular velocity w is called the precessional rate. 

The principle obeyed by the motion of the gyroscope during the appli- 
cation of an external moment may be formulated in the following manner: 
If a gyroscope, having three degrees of freedom, has a moment applied to 
it about an axis perpendicular to the axis of rotation, the gyroscope will 
turn In such a manner as to bring the vector of the main rotation Q into 
coincidence with the moment vector M along the shortest path (see Fig. 
9.11). 

Strictly speaking, this law contains all of the “‘unusual’’ behavior of 
the gyroscope. The gyroscope axis turns not in the plane of the applied 
moment but in a plane perpendicular to it. Besides, the precessional mo- 
tion is not accelerated. Precessional rate w increases only while the applied 
moment increases. With a constant moment the angular velocity of pre- 
cession remains unchanged and, after the action of the moment is stopped, 
precession ceases. 

If initial rotation with velocity 2 were not imparted to the gyroscope, 
it would behave as an ordinary body with the application of an external 
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moment. From Fig. 9.11 it is seen that in that case (when N = 0), vector 
AN would entirely determine the motion of the flywheel, and the last, due 
to the action of the moment, would be rotating with an acceleration about 
axis x’, 1.e., in a plane of the action of the moment. 
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Fig. 9.12. Explanation of precession by means of Coriolis acceleration. 


The behavior of the gyroscope may also be quite simply explained by 
examining accelerations which develop during the precession. 

For simplicity’s sake let us examine a flat disk rotating with velocity 
Q (Fig. 9.12). Let us suppose that by some means precessional motion about 
axis y’ has already been imparted to the disk. 

Due to the smallness of w, the velocity of an arbitrarily chosen point A 


v = Qr. 
The component of this velocity perpendicular to axis y’ 
Vs = Qr cos go. 


The turning of the disk about axis y’ gives rise to Coriolis acceleration 
at point A, whose absolute value is 


7 = Qwrz = 2wOr COS g. (9.3) 
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The vector of Coriolis acceleration is perpendicularly directed to the 
plane of the rotating disk. By varying r and ¢ it is possible to find accelera- 
tions 7 at each point in the disk. The distribution of these accelerations on 
the surface of the disk (acceleration diagram) is presented in Fig. 9.13. 
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Fig. 9.13. Distribution of Coriolis acceleration in the plane of a disk. 


In order to impart this system of accelerations to the disk it is neces- 
sary, as we see, to apply to it a moment about axis x’ perpendicular to the 
precession axis y’. 

The magnitude of the moment 


M= h, jy'dm 


where y’ = r cos g is the moment arm of the element, and dm is the element 
of the mass. 

Integration is done for the entire volume of the disk V. Considering 
Fq. (9.3) we obtain 


M = WQ Í; r? cos? gdm = 2wQ h ydm. 


But the last integral is the moment of inertia of the disk about axis x’. 
For a flat disk it is equal to half of the moment of inertia about axis z. 


h y?dm = A = C/2. 
In this way we arrive at expression (9.2) 
M = Cw. 


Knowing the law of precession it is possible to explain and foresee many 
manifestations connected with gyroscopic effect. 
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For instance, let us trace the behavior of a spin stabilized missile in 
flight. 

The missile leaves the launching rails with zero, or almost zero, angle 
of attack. Having high angular velocity during its further motion, it 
behaves like a gyroscope and strives to maintain, unchanged, the direction 
of the longitudinal axis. Due to the curvature of the trajectory an angle of 
attack appears, and with it an aerodynamic moment in the vertical plane 
(Fig. 9.14). 

If the missile is statically unstable (which is usually the case with the 
nonfinned missiles), the aerodynamic moment will be directed so as to 
increase the angle of attack, and the vector of the moment will be perpen- 
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Fic. 9.14. Spin stabilized missile in a Fig. 9.15. Precession of a spin stabilized 
trajectory. missile. 


dicular to the plane of the drawing and directed toward the nose. With 
axial rotation of the missile to the right, vector Q is directed forward. 
Therefore, according to the principle of precession, the missile begins to 
precess to the right. During further motion, the axis of the missile will be 
describing a cone about the trajectory of the center of gravity with the 
motion to the right (Fig. 9.15). For a statically stable missile, the preces- 
sion will be to the left. 


2. APPLICATION OF THE GYROSCOPE 


The property of the gyroscope to maintain the direction of the spin axis 
is widely used in the field of navigation and stabilizing devices. 

This widely accepted “property,” however, is conditional. We already 
know that with the application of an external moment the gyroscope pre- 
cesses and does not maintain unchanged direction of the spin axis. 

A gyroscope installed in a gimbal mount is subject to the moments of 
friction forces in the bearings of the gimbal. These moments develop from 
the rotation of the gyroscopic rotor itself, and also from the rotation of the 
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external gimbals of the gyroscope relative to the rotor during oscillating 
motion of the flying apparatus. 

After the passage of more or less considerable time, the axis of a free 
gyroscope deviates noticeably from the original direction. Even with well- 
designed gimbal bearings, unpermissible deviation of the axis takes place 
within a few minutes of the unrestrained operation of the rotor. Therefore, 
it is clear that in the operating cycle of a gyroscopic device the position of 
the gyroscopic axis must be continuously corrected. 

As a simpler, and at the same time quite ingenious, method of gyro- 
scope correction, the operational scheme of the simplest artificial horizon 
should be examined. 





Fig. 9.16. Schematic of the artificial horizon device: 7, air supply tube; 2, air exit 
ports; 3, 4, 5, 6, dampers; and 7, rotor. 

The artificial horizon instrument is used to show the pilot, while flying 
blind, the position of the horizon. The sensitive element of this device, in 
this case the gyroscope, must constantly maintain its orientation in space. 

Figure 9.16 shows schematically the construction of an artificial horizon. 

The device consists of a gyroscope with a vertical axis mounted in a 
gimbal. 

The rotor of the gyroscope is set in motion by a stream of air entering 
through the axis tube 7, and impinging on the ribbed surface of the rotor. 
From the rotor box the air escapes through four openings 2, halfway closed 
by pendulum dampers 3, 4, J, and 6. 

Let us see how the correction is performed. Let us assume that, due to 
the action of friction forces in the bearings, the gyroscope axis turns about 
axis x’ through an angle A0 (Fig. 9.16). During this motion damper 8 will 
open up and damper 6 will close. The escaping stream of air will create a 
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reacting moment M about axis y’. According to the law of precession, the 
gyroscope will start to turn in order to bring vector Q closer to vector M 
along the shortest path, and the vertical position of the gyroscope axis will 
be restored. An analogous sequence will take place if the gyroscope, due to 
the action of friction forces, turns about axis y’. 

There are other schemes of correction besides the one described. In par- 
ticular, the magnetic needle of an ordinary compass is usually used for the 
correction of a directional gyroscope, and the gyroscope is oriented along 
its direction. 

In examining all of these schemes a natural question comes up. Is the 
gyroscope so necessary if it has to be corrected anyway with the aid of an 
ordinary pendulum or an ordinary magnetic needle? Is it not simpler to 
use the pendulum directly as an indicator of the vertical, and the magnetic 
needle as an indicator of direction? 

The fact is, however, that the pendulum and the needle of a compass 
have very low masses and low inertias, and are very much subject to various 
spurious effects. For instance, should an aircraft go into a turn the pendu- 
lum would orient itself in the direction of total acceleration and the device, 
without a gyroscope, would not indicate the true horizon. The gyroscope 
has considerable inertia and a prolonged systematic disturbance is required, 
for instance, a prolonged bank, for the artificial horizon equipped with a 
gyroscope to show a noticeable error. The greater the kinetic moment of 
the gyroscope CQ, the more inert is the gyroscope, and the more effective 
is its operation. 

Therefore, the gyroscope finds application in apparatus due to its high 
inertia, the low frequency of its own oscillations, and low susceptibility 
to the action of spurious disturbing forces. As opposed to the magnetic 
needle and the pendulum, the gyroscope may be used in instruments as a 
power source setting in motion certain mechanisms. 


C. Gyroscopic Devices and Actuating Components of the 
Automatic Stabilizer of a Long Range Rocket 


1. ARRANGEMENT AND FUNCTION OF THE RUDDERS 
oF A Lona RANGE ROCKET 


We will examine the construction and operation of an automatic stabil- 
izer of a long range rocket, using the self-contained steering apparatus of 
the V-2 ballistic rocket as an example. 

The rocket is equipped with four pairs of rudders. The arrangement of 
rudders is shown in Fig. 9.17. Rudders 1, II, III, and IV are jet vanes, 
while I’, II’, III’, and IV’ are fins. 

Fins and vanes are numbered clockwise, looking forward. The forward 
vane, which is in the plane of the trajectory, is designated by the number I. 
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Fig. 9.17. Schematic arrangement of jet vanes and fins on a V-2 long range ballistic 
rocket. 


Since in guided flight the rocket does not rotate about its longitudinal axis, 
this orientation of vanes is maintained during the entire powered phase of 
the trajectory. 

Vanes II and IV stabilize the rocket in pitch. By turning vanes II and 
IV in synchronism, according to a definite plan, we are able to program a 
rocket in any manner. 
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The function of vanes I and III is to maintain the rocket on a given 
course, l.e., not to permit it to deviate from the plane of the trajectory. 
This is accomplished by synchronous turning of vanes I and III. Roll 
stabilization of the rocket is performed by the same vanes, I and III. If 
the rocket has a small angle of roll, vanes I and III turn in opposite direc- 
tions, thereby establishing a restoring moment about the longitudinal axis. 

Fins I’, II’, III’, and IV’ have an auxiliary purpose and perform the 
same functions as the jet vanes. 

Their installation on the rocket is not essential in principle. Therefore, 
we will not discuss the operation of the fins. 

In examining the operation of the automatic stabilizer we will begin 
with a description of the construction and operation of the gyroscopic 
devices. 


2. THE HORIZONTAL FREE GYROSCOPE 


The first device, called the horizontal free gyroscope, is designed for 
stabilizing the rocket in the plane of the trajectory. This same device also 
programs the rocket. Therefore, the signals from this gyroscope must be 
transmitted to vanes II and IV. 

The construction of the horizontal free gyroscope is shown in Fig. 9.18. 

The gyroscope is mounted on a gimbal in such a way that the spin axis 
of the rotor is horizontal and lies in the plane of the trajectory. 

Rotor 1 also serves as the armature of a nonsynchronous motor, whose 
stator winding is supplied by a 500-cycle alternating current. The rotor is 
brought up to speed a few minutes before rocket launching and automati- 
cally assumes the necessary orientation in space with the aid of a pendulum- 
actuated correcting device. 

If the axis x’ (see Fig. 9.18) inclines to the vertical, the pendulum, 3, 
will close a contact, and the electromagnet, 4, will receive a properly 
directed signal. The electromagnet will create a moment about the vertical 
axis x’ which, as we already know, will cause precession of the gyroscope 
about the horizontal axis y’. This precession will continue until the moment 
about axis x’ becomes zero, i.e., until the pendulum contact, 3, is opened. 

With the deflection of the axis x’ in the other direction pendulum 8 will 
close a contact on the other side. In this case the electromagnet will create 
a moment of an opposite sign. 

Orientation of the rotor relative to the gimbal frame, 5, is accomplished 
with the aid of an analogous device. If the rotor turns about axis x’, then 
one of the contacts, 6, will be closed and a signal will be transmitted to the 
electromagnet, 7. Due to the action of this electromagnet, precession will 
start about axis x’ and the necessary orientation will be restored. 

The system of gyroscopic correction operates up to the moment of 
launch. After launching the correction system is disengaged and the gyro- 
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scope becomes free. During the time of the guided flight (64 sec), the axis 
of the uncorrected gyroscope does not have a change to deviate noticeably 
from the given direction. 

Should the rocket in flight deviate from the required direction in the 
plane of firing in a given trajectory plane, the potentiometer, 8, connected 
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Fig. 9.18. Construction of a horizontal free gyroscope: 1, gyroscope rotor; 2, stator 
winding; 3, pendulum; 4, electromagnet; 5, gyroscope gimbal; 6, contact arrangement; 
7, electromagnet; 8, potentiometer; 9, potentiometer wheel; 10, metal band of the 
programming mechanism; 77, eccentric of the programming mechanism; 12, pawl of 
the step relay; 13, ratchet wheel; and 14, latch. 


to wheel 9, and through it to the frame of the rocket, will rotate together 
with the rocket about a stationary gyroscope, thereby generating an elec- 
trical signal which, having passed through the amplifier-transformer sys- 
tem, will be applied to the actuating components—the steering mecha- 
nisms. Vanes II and IV will turn, and the required rocket direction will be 
reestablished. 
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If, during flight, the potentiometer is turned through some angle Ag 
relative to the frame of the rocket, then obviously the vanes will respond in 
the same manner as if the rocket itself deviated by that angle and will 
turn the rocket through an angle Ag. Therefore, by turning the potentiom- 
eter, 8, according to a given plan, we will cause the turning of a rocket 
according to the same plan in the plane of the trajectory, i.e., we will 
impart to the rocket its program. The turning of the potentiometer is per- 
formed by the so-called programming mechanism (see Fig. 9.18). 

Potentiometer 8, and the wheel, 9, constitute a rigid unit. A thin metal 
band, 10, is passed over the wheel, 9. This metal band is connected on the 
opposite side to the eccentric, 11, which is profiled to correspond to a given 
program. Eccentric 1/1 is rotated through a worm gear transmission by a 
step relay, 12. The last consists of an electromagnet with an armature. 
When the electromagnet receives an impulse, the armature is attracted 
to the magnet and by its pawl moves the ratchet wheel, 13, by one tooth. 
Therefore, the speed of rotation of the wheel, 73, depends on the impulse 
frequency received by the electromagnet. 

The leaf spring, 74, shown in Fig. 9.18, is a latch for the ratchet wheel 
which prevents its rotation in the reverse direction. 


. 3. THe VERTICAL FREE GYROSCOPE 


The second gyroscopic device of the automatic stabilizer is the vertical 
free gyroscope. This apparatus is installed in the rocket to insure stabiliza- 
tion in yaw and roll and controls vanes I and III. 

The construction of the device is shown in Fig. 9.19. 

The axis of the rotor is perpendicular to the plane of the trajectory. 
Therefore, the gyroscope is insensitive to the rocket maneuvers in the plane 
of the trajectory (in pitch), but reacts to roll and yaw. 

As in the case of the horizontal free gyroscope, alternating current is 
supplied to the coil, 2, for the rotation of the rotor, 1. 

The correction of the device is the same as for the horizontal free gyro- 
scope and is also performed only up to the time of launching. 

With the deviation of axis x’ from the vertical, the pendulum, 3, closes 
one of the contacts, resulting in a signal being supplied to the electromag- 
net, 4, and axis x’ returns to the vertical position. If the gyroscope turns 
about axis x’, a signal is supplied by potentiometer 5 to electromagnet 6, 
and the normal position of the gyroscope is again restored. 

After launching, potentiometer 5 transmits a signal, not to the correct- 
ing magnet 6, but to vanes I and III. Since the axis 2’ coincides with the 
longitudinal axis of the rocket, it is obvious that, in flight, potentiometer 5 
will respond to the roll of the rocket and will transmit signals to move vanes 
I and III in opposite directions. 
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Fig. 9.19. Construction of a verticle free gyroscope: 1, gyroscope rotor; 2, stator 
winding; 3, pendulum; 4, electromagnet; 5, potentiometer; 6, electromagnet; and 7, 
potentiometer. 


Potentiometer 7 will transmit signals for correcting direction. These 
signals must be transmitted to vanes I and IIT in a manner that will cause 
them to turn synchronously. 


4. STEERING MECHANISMS 


The steering mechanisms are the actuating components which turn the 
vanes in the necessary direction and by a necessary amount depending on 
the received signal. 
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Figure 9.20 shows a schematic, and Fig. 9.21 the construction of a 
hydraulic steering mechanism. 

The signal obtained from the gyroscope potentiometer, transformed 
and amplified in the intermediate devices of the automatic stabilizer, is 
received by the polarized relay 1 of the steering mechanism. Depending on 
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Fig. 9.20. Schematic of a rocket steering mechanism (see also Fig. 9.21). 


the sign of the signal, the polarized relay moves yoke 2 in the proper direc- 
tion and piston valve distributor 3, connected to it and mounted on a flat 
spring 4, which replaces the mechanical hinge. 

With the rotation of the distributor, one of the pistons (6 or 6) will 
close the bypass, and the oil, supplied by the gear pump, 7, will enter the 
chamber of the power cylinder, displacing the power piston, 8. The force 
will be transmitted to the vane shaft, 11, by the piston rod, 9, and the crank, 
10. The vane, 12, attached to the shaft will turn in the required direction. 
The gear pump is rotated by motor 18 by means of shaft 14. 
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Frc. 9.21. Construction of a steering mechanism: /, polarized relay; 2, yoke; 3, slide 
valve; 4, spring; 5, 6, slide valve pistons; 7, gear pump; 8, power piston; 9, rocker arm; 
10, crank; 11, vane shaft; 12, vane; 13, gear pump motor; and 14, motor shaft. 


The described steering mechanism operated only one vane. There must 
be four such mechanisms installed in the rocket. 


5. CONDITIONS FOR STABLE FLIGHT 


In the example of the simplest automatic steering (see page 279), we 
have already seen that with the deviation of the flying apparatus off course 
the simple tilting of the rudder in the required direction is insufficient for 
the attainment of stable equilibrium. If the rudder tilt is in the proper 
direction but not in the required degree, or is lagging, we always risk having 
an unstable regime of rocket steering—divergent steering or steering with 
full departure from the given course. 

Let us determine the conditions which must be satisfied by a stabilizing 
system in order to insure stable motion. With this in mind, let us turn to 
the equation of rocket motion, Eq. (7.38). 


J +M + LD: + M, = 0 
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and solve for the type of motion that the rocket will have after it has 
received some disturbance, even in the guise of an instantaneous angular 
impulse. 

We will simplify the problem by eliminating from the equation the hinge 
moment M, as a quantity small compared to L,D,. Besides, we will consider 
that the remaining parameters of rocket motion, after imparting to it the 
disturbance corresponding to angle ¢, remain unchanged. So, for example, 
the rocket velocity v and aerodynamic coefficients will be considered as con- 
stant, retaining the value which they had prior to the time of disturbance. 

In the first approximation of stability evaluation this assumption is 
quite acceptable, inasmuch as we are interested only in the deviation of the 
rocket from the steady state motion and not in the motion of the rocket 
itself. 

The initial equation of motion, Eq. (7.3), being analyzed, takes the form 


Jo ELD: =: (9.4) 


With a disturbance, the quantities entering into this equation change, 
and we get 


where Ag is the angle of deviation from the programmed angle resulting 
from disturbance; AM, the corresponding change in aerodynamic and damp- 
ing moments; and AL,, the change in vane lift force. 

Subtracting Eq. (9.5) from Eq. (9.4) we find 


JA + AM + AL,D, = 0. (9.6) 


Moment M is composed, as we know, of stabilizing and damping mo- 
ments. The first of these is proportional in the first approximation to the 
angle of attack, and the second to the angular velocity of the rocket rota- 
tion. Consequently, the change in the moment AM will be composed of two 
addends, one of which is proportional to the angle of rotation Ag, and the 
second to the velocity of rotation A¢g 


AM = adg + bAg (9.7) 
where a and b, due to our assumption of the constancy of flight parameters 
of disturbed flight, are constant values. 

As to the lift force Z,, we already know that, in the first approximation, 
it 1s proportional to the vane rotation angle 8. Consequently 
AL, = cAg 


where A@ is the additional angle to vane rotation due to the reaction of the 
automatic stabilizer to the disturbed motion of the rocket; the quantity 
AB is a function of time. 
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Therefore, Eq. (9.6) takes the form 
JA¢d + bAg + aAgy + cD,AB = 0. (9.8) 


Now the problem consists in finding out in what way the angles of vane 
rotation Af and rocket disturbance Ag are interrelated. If they are so related 
that the solution of Eq. (9.8) will yield decaying or, even better, the rapidly 
decaying solution for Ag, the rocket steering will be stable. If in solving 
Eq. (9.8) we discover that function Ag does not decay with time, it will 
indicate that automatic stabilization does not insure stable motion. 

The angle of vane rotation and the hinge moment M, are related by 
the following equation 


M, = 8 + ab + ağ 


where Qo, Qı, @2 are constant quantities. 
Consequently 


AM, = aAB + a,AB + œA. (9.9) 


However, the variation in the hinge moment AM, depends on the sig- 
nals which are transmitted to the steering mechanism from guidance 
devices. Depending on the method of reception and transmission of signals 
from the guidance gyroscope, the addition of the hinge moment AM, may 
be, for example, proportional to the angle of deviation of the rocket axis 
Ag from the program. It is possible, for instance, to arrange for the value 
AM, to be proportional to the velocity or the acceleration of the angular 
disturbance of the rocket motion (A¢g or Ag). 

Let us assume, for instance, that 


Coefficients bo, b1, and be are constant quantities and, for the time being, 
undetermined, but it is they which characterize the behavior of the steering 
system. The right hand side of Eq. (9.10) is called the time constant. Intro- 
ducing design changes into the steering system, we can change values bo, 
bı, and bz and achieve stable motion of a rocket by the stabilizing system. 

Comparing expressions (9.9) and (9.10) we find 


AB + AB + AB = bAy + DAG + bA. 


Let us substitute into this the expression for A8 from Eq. (9.8). As a 
result we will get one equation in terms of Ag. 


Ag’ + d:Ag’ + dA + dAg + doâ = 0 (9.11) 
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where 


do = aoa bocD, 
od od 

me bao + aa , bieD, 

dg dad Aad 
_ Ad + ab + ana Pa bcD, 
7 dod “dad 
= Ay + ab 
7 dad f 





(9.12) 
də 


ds 


Thus, the principle of the disturbed rocket motion will be determined 
by the solution of a linear differential equation, Eq. (9.11). 
Solution of Eq. (9.11) we will seek as usual in the form. 


Ag = Aet. 
Substituting this expression into Eq. (9.11), we will get the typical 
equation 
kt + dak? + dak? + dik + do = 0. (9.13) 


Equation (9.13) has four roots: ky, ke, ks, and k4. As a result, the solution 
of Eq. (9.11) will take the form 


Ag = Ayeh*t + A*t + Age? + Asem, (9.14) 


Let us see now under what conditions the motion will be decaying. 

The roots of Eq. (9.13) may be either real or complex. 

For real roots of k, the function e*t will be divergent for positive values 
of k and convergent for negative values of k. 

If the root k is complex 


k=ptu 
where u and v are real numbers. Then, according to Euler’s identity 
ekt = eutvdt = ext(cos pt + 7 sin vt). 


Each algebraic equation with real coefficients yields only an even num- 
ber of complex roots, whereupon these roots, as is known, are conjugate 
pairs. 

Suppose, for instance, that the roots kı and ke are conjugate 


k= u+ 
ka = u — vt. 


300 PRINCIPLES OF STABILIZATION AND STEERING 


Then, for the first two addends in the right hand side of Eq. (9.14) we will 
find 


Aet + Ase! = Aje##(cos vt + 72 sin vt) + Are#**(cos vt — 7 sin vt) 


or 
Aye! + Ase? = (Ay + Aert cos vt + (A — As)e#! sin vt. 
Denoting 
A, + A» T Ci i(Áı a A?) za Co 
we get 
Ayes! + Ave! = Cet! cos vt + Cæ! sin vt. (9.15) 


If the other two roots, k; and ky, in Eq. (9.13) are also complex, then the 
sum A3e*! + A,e*! on the right side of Eq. (9.14) will be transformed 
analogously. 

The obtained function (9.15) will be damped if u is negative. The case 
of u = 0 indicates oscillating undamped motion; the case of u > 0 indicates 
oscillations with increasing amplitude. 

Therefore, the conditions for stable steering may be formulated in the 
following form: For stable steering it is necessary for the real parts of all 
roots of Eq. (9.13) to be negative. This condition may be fulfilled by choos- 
ing, in a proper manner, Eq. (9.13) and coefficients do, di, dz, and dz, which 
depend on the parameters of the steering system. 

The answer to the question of whether the steering system is stable or 
unstable may be obtained without solving Eq. (9.13). Using methods of 
higher algebra, it is possible to determine relationships between coefficients 
do, dı, dz, and d; for which the real parts of all roots of the algebraic equation 
will be negative. For the equation of the fourth degree, Eq. (9.13), it is 
necessary and sufficient to satisfy the condition of Routh according to 
which the real parts of all roots will be negative if simultaneously 


dy > 0, dı > 0, d > 0, d > 0 
and 
didods — dod — d,? > 0. 


It is not difficult to be convinced of the fact that the motion of a guided 
rocket may be either stable or unstable. 

For instance, let a signal be transmitted to a steering mechanism which 
is proportional to the deviation of the rocket, and which is independent of 
the angular velocity or acceleration of the rocket. Then, as it follows from 
expression (9.10), bı = b: = 0. Further, let us suppose for the sake of sim- 
plicity that the vane masses, and moments of inertia during their rotation, 
are small, that the vanes themselves are totally compensated in the flow, 
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and that the hinge moment is independent of the angle of rotation. This 
means that, in Eq. (9.9), ao = az = 0. Besides, if it is assumed that the 
rocket has left the region of the atmosphere and the role of aerodynamic 
forces is negligible, then, in expression (9.7), a = b = 0. 

Multiplying all terms of Eq. (9.11) by &J, and keeping in mind the 
notation of Eq. (9.12), we get 


adJAg + becD,Ag = 0 


or 
Ag’ + Ago = 0 
where 
@ = bocD,/aJ. 
The typical equation will have the form 
B+ 2 = 0 
whence 


k= —d, keg = 4d + tid V3. 


The real part of roots kz and ks turns out to be positive. Consequently, 
the motion of the rocket will be unstable. This came about not because we 
have neglected vane inertias and aerodynamic forces, but primarily because 
the chosen principle of regulation, with the transmittal of signals to the 
vanes proportional only to the angle Ag, was incorrect. 

In order to obtain stable motion, it is sufficient to transmit to the steer- 
ing mechanisms signals which are proportional to Ag and, for the improve- 
ment of the control system, also proportional to Ag. Then, in expression 
(9.11), bı and be will not be equal to zero. By the proper selection of these 
coefficients it is easy to obtain a stable system of stabilization. 

A stabilizing system of the V-2 long range ballistic missile is based 
precisely on this principle. Signals obtained from the potentiometers of the 
horizontal and vertical free gyroscopes are differentiated by special devices 
and are transmitted to the vanes in the form of functions of deviation angle 
Ag, Eq. (9.10). 

Let us turn to the examination of intermediate devices which perform 
the transmission of signals from the gyroscopes to the steering mechanisms. 


D. Intermediate Devices of the Long Range Rocket 
Automatic Stabilizer 


1. THE GENERAL SCHEMATIC OF THE INTERMEDIATE DEVICES 


The intermediate devices of the automatic stabilizer are intended for 
transformation and amplification of signals obtained from gyroscopic 
potentiometers, and their transmission to the steering mechanisms. 
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The basic schematic of the intermediate devices of the automatic 
stabilizer is shown in Fig. 9.22. 

Stabilization in the plane of the trajectory with the aid of the hori- 
zontal free gyroscope is performed independently of stabilization in yaw 
and roll. Therefore, signals transmitted by the horizontal and the vertical 
free gyroscopes have independent systems of transformation and trans- 
mission to the vanes. 

The signal in the form of de voltage obtained by the horizontal free 
gyroscopic potentiometer is transmitted, first of all, to a differentiating 
circuit where it is transformed into the form 


vy = oAyg + (Ag + Ag. (9.16) 


The function of the differentiating circuit consists, therefore, in adding 
to the voltage, which is proportional to the angle of rocket deviation, com- 
ponents which are proportional to the angular velocity and angular accel- 
eration of the rocket. These additional signals, as we already know, are 
necessary in order to insure steering stability. 

The potentiometer is supplied by a source of direct current, since it is 
very difficult to design a differentiating circuit which would operate on 
alternating current. 

After the differentiating circuit, the signal is transmitted to a converter, 
where it is converted into an amplitude-modulated alternating current sig- 
nal with a frequency of 500 cycles. This conversion is performed in order 
to simplify the subsequent amplification, since amplification of a direct 
current signal poses great difficulties. 

The converted signal is transmitted through a transformer to the ampli- 
fiers, which control the steering mechanisms of vanes IT and IV. 

Since vanes IT and IV are not mechanically interconnected, a synchron- 
izing device is introduced into the pitch system of stabilization in order to 
insure synchronous operation of vanes II and IV, excluding, in this manner, 
the occurrence of undesirable perturbations in roll. When vanes II and IV, 
together with the synchronizing potentiometers, are turned to the same 
angle, no signal is transmitted to the primary coil of the transformer 7'p,. 
When the vanes are turned to a different angle, however, the primary coil 
of the transformer Tp, will receive a signal proportional to the difference 
of the vane angles. This signal will then be transmitted to the amplifiers 
and, through them, to the relay of the steering mechanisms. As a result, the 
coordination of the vanes will be reestablished. 

The system of control for vanes I and III is analogous to the control 
system of vanes II and IV. Signals from two potentiometers are separately 
differentiated, converted, and transmitted through transformers to the 
amplifiers of the steering mechanisms of vanes I and III. The signal from 
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the roll potentiometer forces the vanes to defiect in opposite directions, 
and a signal from the yaw potentiometer is transmitted in such a way that 
vanes turn in the same direction simultaneously. 

In the schematic of Fig. 9.22 an extra spare coil is shown in transformer 
Tp:. Through it may be impressed a signal from radio control. This signal, 
obviously, as the signal from the yaw potentiometer, will correct direction 
of rocket flight. 

Let us examine separate elements of the intermediate devices. 


2. Tue DIFFERENTIATING CIRCUIT 


Voltage v from the gyroscopic potentiometer is proportional to the angle 
of rocket deviation Ag. Since Ag is a function of time, v is also a function of 
time. 

Differentiation of the signal v may be performed with the aid of the 
simplest RC circuit, 1.e., by introducing in sequence capacitance and ohmic 
resistance (Fig. 9.23a). 


R, 





(a) (b) 


Fic. 9.23. Schematic of the simplest differentiating circuit. 


The signal being differentiated is applied across input terminals and its 
derivative is picked off at output terminals. Obviously, 


Vin — Ve = Vout 


where v. is the voltage drop across the condenser. 
The current which passes through the capacitance 


te = C(dv,/dt). 


If no load is connected to the terminals, i.e., if the circuit operates with 
an amplifier, then 


Vout = aft. 
Eliminating from the three obtained equations v. and te we get 


Vout + RC (dvour/dt) = RC(dvin/dt). 
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If the selected value of RC is sufficiently small, it is possible to neglect 
the second addend in the left part of this equation, and we will then have 
almost pure differentiation 


Vout © RC (dvin/dt) = RCvin 


i.e., the output will be a voltage almost proportional to the derivative of the 
input signal. 

It is essential to note that, inasmuch as RC must be very small, the 
converted signal, as it follows from the last expression, is sharply reduced 
in strength by the differentiating circuit. Therefore, the introduction of 
amplifiers further on becomes necessary. 

If we wish to obtain an output signal which depends not only on the 
derivative dv;,/dt, but also on the function v;, itself, then it is necessary to 
add to the simplest differentiating circuit the resistance A, in parallel with 
capacitance C (Fig. 9.23b). The output signal will then have the form 


Vout = MWin + Movin 
where m, and Mm: are constants depending on parameters R, A, and C. 


Therefore, the output voltage in a circuit becomes proportional to volt- 
age Vin and its first derivative, and the loop shown in Fig. 9.23b performs a 





Fig. 9.24. Differentiating circuit schematic of the automatic stabilizer. 


single differential transformation of the signal. A scheme shown in Fig. 
9.24 is used for double differentiation of a signal. The output voltage of 
this circuit is 

Your = Nin + Nin + NaVin 


where no, ni, and ne depend on the parameters of the circuit. 
This signal, converted and amplified, finally reaches the steering 
apparatus in the form of a demand to vary the hinge moment, Eq. (9.10) 


AM; = biAg + bidg + beAg. 
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The relationships among coefficients bo, 61, and be are determined by the 
relationships among coefficients no, nı, and ne. Therefore, the relationships 
among coefficients bo, b1, and be necessary for stable steering are obtained 
by suitable selection of resistances and capacitances of the differentiating 
circuit. 


3. SIGNAL CONVERTER (MODULATOR) 


In passing through the differentiating circuit the signal is drastically 
weakened and must be amplified. Since amplification of a direct current 
signal poses considerable difficulties, in the system of the automatic stabil- 
izer under discussion the signal is converted to an alternating current 
modulated in amplitude. A scheme shown in Fig. 9.25 is used for this. 


5007 
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Fic. 9.25. Converter schematic. 


The voltage from the differentiating circuit is applied to the center taps 
of the coils of the two transformers: the first, Tp, input, excited by 500- 
cycle alternating voltage, and the second, Tpz, output, from which the 
converted signal is picked off. 

The basic element of this scheme is a hrid of four selenium rectifiers. 

A selenium rectifier consists of a stack of iron or aluminum disks covered 
on one side with a layer of selenium. In its turn a layer of some fusible metal 
is deposited on the selenium. A thin layer at the boundary between the 
selenium and the fusible metal has the property of being a good conductor 
in one direction and a very poor conductor in the opposite direction. The 
graph in Fig. 9.26 shows the amount of current passing through a selenium 
rectifier with respect to the applied voltage. 
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The selenium rectifier bridge ABCD (see Fig. 9.25) is made strictly 
symmetrical. In the absence of an input signal at terminals EF it is bal- 
anced. Notwithstanding the fact that the transformer Tp, impresses a 
voltage across the diagonal AC, the voltage in the second diagonal BD con- 
nected to the primary coil of the output transformer Tp: is equal to zero. 
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Fra. 9.26. Characteristic of a selenium Fig. 9.27. Change in the output signal 
rectifier. during one cycle of the alternating current 
carrier depending on the polarity of the 
signal. 


Upon transmittal of a signal the symmetry of the bridge, because of the 
nonlinear characteristic of selenium rectifiers, is disturbed and there is a 
voltage across the diagonal BD proportional to the magnitude of the input 
signal. Voltage also appears across the secondary coil of the transformer 
Tpz. This will be an alternating voltage whose amplitude is proportional 
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Fig. 9.28. Modulated signal. 


to the magnitude of the direct current input signal. If the direction of the 
input signal current changes then the phase of the output alternating volt- 
age will change by 180° (Fig. 9.27). 

As a result we get in the secondary coil of the output transformer an 
alternating voltage with a frequency of 500 cycles, and an amplitude 
variable with the signal arriving from the differentiating circuit (Fig. 9.28). 
The phase of the output alternating voltage depends on the polarity of 
the input signal. 
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4, AMPLIFIER AND DEMODULATOR 


A schematic of the amplifier and demodulator is shown in Fig. 9.29. 

The basic element of the amplifier is an electronic tube (pentode), 
widely used in the most diversified devices and systems. The tube is sup- 
plied from a voltage source with 200-250 volts. 

Voltage from the secondary coil of the output converter transformer is 
impressed across the grid of the tube. This voltage, amplified by the tube, 
is supplied to the coil, A, of a three-cored transformer (coils of the trans- 
former are wound on three cores). 


iach 
N) 27.7: 


+v 


g l l ea 


modulator 


Y 
-V4 Fic. 4.29 


Fic. 9.29. Amplifier schematic. 


The amplified voltage is the voltage of the alternating current with a 
frequency of 500 cycles, whose amplitude is proportional to the output 
signal of the differentiating circuit and whose phase depends on the polarity 
of the signal. 

Since, in the final analysis, the signal must control a polarized relay of 
the steering mechanism which reacts only to the direct current signals, it 
is necessary for the amplified signal to be reconverted into direct current 
signal. This reverse, second, conversion of the signal (or demodulation) is 
performed by the three-cored transformer and selenium rectifiers connected 
as shown in Fig. 9.29. 

The three-cored transformer, besides the amplified signal voltage (coil 
A), is supplied by the alternating voltage with a frequency of 500 cycles 
from the power source. The last is called a carrier, since the phase of the 
amplified signal is determined by comparison with it, i.e., the sign of the 
original signal obtained from the gyroscopic potentiometer. 

If the steering signal is absent in coil A, the voltages across coils E and 
D are equal to each other. 
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Voltages across coils Æ and D are rectified by selenium rectifiers and 
are supplied to the two coils of a polarized relay. The rectifiers are so con- 
nected that the currents in each coil oppose each other. For a zero signal, 
both currents 2, and t: (see Fig. 9.29) are equal in magnitude and the total 
net effect is equal to zero. Consequently, the armature of the relay remains 
at rest. 

The three-cored transformer operates in such a manner that, depending 
on whether the currents in coils BC and A correspond or are opposite in 
phase, the voltage in coil Æ correspondingly increases or decreases (while 
the voltage across coil D correspondingly decreases or increases). 

As the signal is transmitted, one of the currents (2; or t2) will increase 
and the other one will decrease, depending on the sign of the signal. The 
polarized relay will then operate either in one direction or the other. The 
net result is the turning of the vane, depending upon the angle, angular 
velocity, and angular acceleration of the rocket. 


E. Range Control 


1. METHODS oF RANGE CONTROL 


Range is one of the elements of flight trajectory and is determined, 
generally speaking, by the parameters of rocket motion. 

If it is possible to interfere with the motion of the rocket along the tra- 
jectory, the reaching of a given target and the necessary range may be 
obtained with various relationships of flight parameters. To put it another 
way, the rocket may be guided to a given target by various trajectories. 

However, starting with the time when the rocket ceases to be guided 
and continues its motion subject to forces which cannot be changed by us, 
the flight trajectory no longer depends on our wishes and is almost entirely 
determined -by those initial conditions of the coasting flight which were 
imparted to the rocket at the end of the powered phase. Therefore, for a 
rocket to reach a given target it is necessary that the parameters at the 
beginning of coasting flight have definite relationships. 

The initial parameters determining the motion of a ballistic body are, 
obviously, velocity vo, the direction of velocity vector (for a plane motion, 
angle ĝo), and the coordinates at the end of the powered phase x» and yo 
(Fig. 9.30). Besides, the rocket flight during the coasting phase may be 
affected by such factors, variable with time, as wind direction and meteoro- 
logical environment for the entire coasting phase or during that part of the 
trajectory which takes place in sufficiently dense regions of the atmosphere. 
The trajectory of flight is also subjected to systematic secondary factors, 
for instance, the velocity of earth’s rotation and the variation in the magni- 
tude of gravity acceleration vector, depending on the latitude of the locale. 
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Therefore, the rocket range is a function of several variables: 
L = f (Xo, Yo, Vo, bo, c) (9.17) 


where parameter c is a correction which takes into account the effect of the 
enumerated secondary factors. 

The range control of powder rockets is performed simply by means of 
changing the elevation angle of the launcher during launch (Fig. 9.31). 





Fic. 9.30. Rocket trajectory and range. 


With this method of range control it is obviously necessary that the 
thrust and operating time of the motor do not vary from rocket to rocket. 
Otherwise, for a constant angle ĝo, there will be a large variation in range. 
In order to avoid this dispersion, increased requirements are imposed upon 
powder grains in dimensional stability and weight and also in accuracy of 
fabricating powder chambers and nozzles. 





Fic. 9.31. Variation in range with the projection angle 6o. 


However, even with full conformity to the indicated conditions there is 
always range dispersion, inasmuch as for individual rockets, variations in 
ignition, combustion regime, and afterburning of a powder charge are 
always possible. With small rockets, this dispersion may be compensated 
by the number of rockets released to cover a given area. For long range 
rockets, because of their high cost, such compensation cannot be allowed. 
For these rockets it is necessary to introduce methods of range control which 
would insure substantially greater accuracy. 
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In examining expression (9.17), it is not difficult to note that the range 
control of a rocket may be attained only by variation in quantities 4%) and 
vo, Since the parameter c is a quantity not subject to arbitrary change and 
variation in coordinates zo and yo may influence the range only within 
restricted limits. 

If we are to control the range of a ballistic rocket by changing angle 4, 
assuming that the constancy of velocity vo is insured by the operating 
time of the motor with a given thrust, then obviously we will get a range 
dispersion of the same order as for nonguided powder rockets. At the same 
time, the necessity for changing the program with range will be an inescap- 
able inconvenience. 

It is more rational to control range by means of changing velocity vo 
at the end of the powered phase. 

The first, somewhat coarse, method of control consists of switching off 
the motor (cut-off of fuel components) at an appropriate time, depending 
on the required range. This method of control, based on the assumption of 
invariable parameters for a group of rockets, does not, however, yield 
satisfactory accuracy. The thrust variation of actual samples of the same 
motor is noticeably greater than the allowable dispersion which is required 
here. The thrust regulation, in order to maintain strictly its constant value 
for a group of rockets, involves considerable technical difficulties. 

Another, more acceptable, method of range control is the fuel component 
cut-off according to actual attained velocity (and not according to calcu- 
lated time function). Thus, the second quantity (vo) of the two (6) and vo), 
which more strongly affects range, is under direct control. As to the second- 
ary factors zo and yo, the influence of the change of these quantities on range 
need not be taken into account in the first approximation. This method of 
control yields considerably greater range accuracy compared to the motor 
cut-off with time. 

The velocity of the rocket may be measured with the aid of radar, 
utilizing the Doppler effect as its operating principle. The device measures 
the difference in frequency of radio waves transmitted to, and reflected 
from, the rocket. This frequency difference, as is known, depends on the 
velocity of the reflecting object, in this case a rocket. Motor cut-off is per- 
formed at the required time by radio beam. 

This method of range control yields high accuracy. However, it requires 
a complicated apparatus, and like any other method of radio control is 
subject to the action of countermeasures (interference). 

A simpler method and one free of the indicated drawback is the method 
of mechanical integration of axial rocket accelerations during the powered 
phase. The device which performs the acceleration integration and thus 
determines the velocity of a rocket is the so-called integrator. 
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2. CONSTRUCTION AND OPERATION OF A GYROSCOPIC INTEGRATOR 


The most feasible, for application, gyroscopic integrator is shown in 
Fig. 9.32. 





Fig. 9.32. Construction of an integrator: 1, gyroscope rotor; 2, stator winding; 3, 
slip rings and slide contacts; 4, rotor housing; 8, gyroscope suspension axis; 6, external 
yoke; 7, gear; 8, motor cut-off cogs; 9, latch; 10, range setting; 11, gear; 12, 13, zeroing 
contacts; 14, erection contacts; and 15, erection motor. 


The integrator has a gyroscopic rotor 1, which is actuated by alternating 
current supplied to the stator winding 2. The current and signals from the 
integrator are conducted through rings 3 with sliding contacts. 
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The rotor is located in housing 4, which is free to turn about axis 6. The 
weight of the rotor and housing creates a moment about this axis. 

The axis 6 is connected with the external gyroscopic gimbal 6, which 
may turn about axis xz’. Through a train of gears the rotation is transmitted 
to disk 7, with cogs 8. 

A few minutes before launching, the gyroscope is brought up to speed. 
During this time the gyroscope axis is Immovably locked by means of a 
detent 9. In Fig. 9.32 the detent is brought out to the side. It consists of 
an electromagnet and two plates. When a signal is transmitted to the 
electromagnet (this takes place at the instant of rocket take-off from the 
pad), the plate of the latch frees the axis of the gyroscope and the last, 
together with the housing, becomes supported by axis 5. 

Precession about axis x’ begins due to the action of weight forces. 

For a stationary rocket, according to Eq. (9.2), the angular velocity of 
precession is 


w = Pa/CQ 


where P is the weight of the rotor and housing; and a is the moment arm of 
force P about axis 5. 

For a rocket moving straight up (as we will assume for the time being 
for the sake of simplicity), with acceleration 7, the apparent force of the 
weight will increase by comparison with P by an amount proportional to 
a/g, and then 


w = (Pa/CQ)(1 + 3/9). 


The angle of gyroscopic rotation about axis x’ in time ¢ 


v= Í, i wdt. 
If the kinetic moment of the gyroscope remains constant, then 
Pa Pa v 
v= CO (1 -+ i) at = Pafi -+ 2) (9.18) 


i.e., the angle of gyroscopic rotation about axis x’ has a linear relationship 
dependent on the flight velocity of the rocket at a given instant of time. 

If the time ¢ of rocket flight is known, then by eliminating addend 
(Pa/CQ)é in the last expression, it is possible to find the rocket velocity by 
means of angle # after a period of time t after launching. 

The rocket range setting is performed by a knob, 10, on the integrator. 
The disk, 17, with a contact arrangement, rotates through the required 
angle relative to cogs 8. This setting is performed, taking into consideration 
the correction due to the first addend in expression (9.18). This correction 
is easily introduced, since the rocket program is known in advance. 
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When the rocket attains the required velocity, one of the cogs, 8, will 
initiate a signal to switch the motor into a preliminary stage, and the second 
cog will initiate total cut-off. 

During operation, the spin axis of the gyroscope must always remain 
perpendicular to axis x’. In order to fulfill this requirement a correcting 
device is made part of the integrator. 

If the spin axis of the gyroscope deviates either up or down, either the 
upper or lower contacts, 74, will be closed and the motor, 18, will receive 
a signal of the corresponding sign. Thus, a moment will act on the gyroscope 





Fic. 9.383. Determination of the apparent axial acceleration. 


about axis x’ through a gear train. Such a moment, as we already know, 
will neither accelerate nor decelerate the rotation of the gyroscope relative 
to axis x’, 1.e., it will not vary the conditions of integration, but will only 
force the gyroscope to turn about axis 5 and resume the proper position. 

The integration of accelerations in the manner described has its defects. 
One of these is the necessity of taking into account the rocket program in 
determining range, which leads to the so-called precedural error of the 
integrator. 

Axis x’ is parallel to the rocket axis. Consequently, the integrator will 
react, not to the accelerations tangential to the trajectory, but to the axial 
rocket accelerations. 

As is seen in Fig. 9.33 the apparent weight force acting on an integrator 
rotor in the axial position is 


P’ = mlj cosa + vó sina + g sin (0 + @)|] 


where m = P/g is the mass of the gyroscopic rotor and housing. 
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In this case the precessional velocity of the gyroscope 
w = (Pa/gCQ)[j cos a + vå sin a + g sin (6 + @)]. 


The angle of rotation of the gyroscope relative to axis x’ 
3 f ibe Í Y 6 si in (0 dt 
= fo = CG lod cosa TP sina + g sin (@ + a)]dé. 


Therefore, we see that the integrator does not yield, in a pure form, the 
rocket flight velocity, and its indications require corrections on the basis of 
the required program 6 + a = f(t) and the calculated change in the angle 
of attack a along the trajectory. Since, in flight, the program is executed 
with certain errors, these errors are reflected, to some degree, in the deter- 
mination of the rocket velocity v by the integrator. 


F. Rocket Dispersion 


1. THe MEASURE OF DISPERSION 


The same indices as used in artillery practice may be used to evaluate 
the degree of accuracy of rocket firing. 

Let us assume that we are firing at a certain target with an unlimited 
number of similar rockets with constant launching conditions (Fig. 9.34). 





Fic. 9.34. Dispersion of rocket fire. 


In this case, when we speak of ‘similar rockets” and “constant condi- 
tions”? we understand by this only “similarity” and “constancy” within 
those limits which can be insured by technological considerations and by 
practical conditions of conducting the firings. 

In actuality, the geometric dimensions of the rocket will vary somewhat 
from round to round. The conditions of powder grain burning in one rocket 
will differ from conditions of burning in another, if only slightly. The con- 
dition of launching itself will not remain invariant. Accidental, even 
though completely unnoticeable, displacement of the launcher from shot 
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to shot and the random wind gusts during rocket flight will change in some 
degree the conditions of launching. 

However, in speaking of similar rockets and constant conditions we 
at the same time keep in mind that systematic errors and systematic devia- 
tions from certain average standards are excluded from the technical 
rocket equipment, their geometry, and conditions of launch. 

For instance, the inaccuracy of the master gage, which might lead to 
the systematic deviation in the geometric dimension of a rocket or its 
charge, is eliminated. The constant wind velocity, which affects the range, 
is taken into account, etc. Therefore, in firing similar rockets with constant 
launching conditions, the deviation of the rocket from the target will have 
a random character and the points of impact shown in Fig. 9.34 will be, to 
a greater or lesser extent, scattered relative to the target. For a large num- 
ber of shots, the random deviation from the target or, more accurately, 
the distribution of missiles on a plane, will follow the Gauss random dis- 
tribution law. 

If, in plane z-y (see Fig. 9.34), a rectangle is isolated with sides Ax and 
Ay, then the number of missiles AN hitting this rectangle, with a large 
number of firings, will be proportional to the area Az - Ay 


AN = Ardy. 
The proportionality coefficient ® is a function of coordinates x and y 
$ = (x, y). 
This is the distribution function. According to the Gauss law 
& = Agren (9.19) 


where A, k, and k are certain constants and where the quantities k and h 
determine the degree of dispersion along the axes x and y, respectively. 
The surface defined by Eq. (9.19) is shown in Fig. 9.35. 


$ 
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Frc. 9.35. Surface described by the dis- Fig. 9.36. Family of dispersion ellipses. 
tribution function (zy). 
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The function © is a maximum at the origin of the coordinates when 
x = y = 0, and becomes zero when z = + © or y = + œ. 

Let us determine the geometric locus of points of equal distribution. 
Assuming ® = C = constant, from Eq. (9.19) we will find 


ka? + h2y? = In (A/C) 


i.e., the equation of an ellipse with axes coinciding with the axes of the 
coordinates. The dimensions of the semi-axes of the ellipse depend on the 
value of C. Therefore, the family of contours of the surface ®(z, y) consists 
of a family of ellipses of dispersion (Fig. 9.36). 





50 % of shells 


Fig. 9.37. Determination of B,. 


As a measure of missile dispersion during firing, the axis dimensions of 
one of these ellipses may be used, for instance, an ellipse within whose 
boundaries half of a large number of missiles impacts, 1.e., an ellipse of the 
best half. 

In artillery practice the measure of dispersion is not the dimension of 
an ellipse, but the width of a strip within whose bounds half of a large 
number of shells impact. The dimensions of this stip in the direction of firing 
(down range) are designated by 2B, (Fig. 9.37) and, sideways, by 2B,. 


2. DISPERSION OF UNGUIDED ROCKETS - 


At the present time, the dispersion of unguided rockets is many times 
greater than the dispersion of ordinary artillery shells. This is explained, 
primarily, by the fact that the powered phase of the unguided rocket flight, 
i.e., that phase during which it accelerates up to speed, is not sufficiently 
stable, while the artillery shell, during acceleration, i.e., in the gun barrel, 
has a strictly determined direction of flight. The above is also confirmed 
by the fact that the unguided rockets, whose burning time ends while still 
on the launching rails, have dispersion approximately of the same order of 
magnitude as artillery shells. 

Therefore, in’discussing the sources causing rocket dispersion it is neces- 
sary to differentiate between those which act during the powered and during 
the coasting phases. 
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The basic causes of rocket dispersion during the powered phase are the 
eccentricity of application of the motive and aerodynamic forces and 
external disturbances which force the rocket to turn in flight about a trans- 
verse axis. 

With an eccentrically applied thrust force, there develops a constantly 
acting moment, forcing the rocket to move continuously along a curved 





Fig. 9.38. Rocket dispersion due to eccentrically applied thrust (thrust misalign- 
ment). 


trajectory and deflecting the rocket from the target (Fig. 9.38). In other 
words, during the powered phase, it is essential that the direction of a 
motive force constantly follow the direction of the turning rocket axis. If 
this force deflects the axis of the rocket from a given direction then it itself 
follows the turn and deflects the rocket from its target even further. 


Direction of 
rotation 
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Fic. 9.39. Effect of cross wind on the flight of an unguided rocket. 


Aerodynamic eccentricity acts during the powered phase in an analogous 
manner, inasmuch as it changs the direction of the rocket axis and, at the 
same time, the direction of the thrust force. 

The action of external disturbing factors is quite singular. A cross wind 
gust, for instance, forces the rocket to deviate from its course, as opposed 
to an artillery shell, not with the gust but against it. Figure 9.39 shows a 
rocket in flight with a cross wind gust acting upon it. If the rocket is stati- 
cally stable (and it is specifically made so), then it will turn into the result- 
ant flow like a weather vane and will continue its flight, due to the action 
of thrust, in a new direction. In this sense, it is not advantageous to have a 
rocket with a large margin of static stability. Such a rocket will be overly 
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sensitive to random disturbances of this type and the dispersion will be 
large. 

In order to reduce the effect of the eccentricity of thrust and aero- 
dynamic forces, the unguided rocket is rotated in flight, i.e., it is made to 
turn slowly relative to the longitudinal axis. This rotation should not be 
confused with the rapid rotation which is imparted to spin stabilized 
missiles. Here, the turning of the rocket is necessary only as a means of 
canceling the moment of the eccentrically applied forces. 

The rotation is imparted to the rocket either by canted rails or by means 
of inelined fins installed on the rocket. In some cases the rotation is imparted 





Fra. 9.40. Rocket spin during powered phase. 


by the reactive force of jets emerging from small auxiliary nozzles drilled 
in transverse rods of the rocket (Fig. 9.40). The rocket slides along the rails 
by means of these rods during launching. 

Aside from the enumerated basic causes which affect the dispersion of a 
rocket there are also other causes which we have touched upon earlier. They 
are fabricating variation in the dimensions of the rocket, the size of the 
charge, and variation in the combustion process itself, whose uniformity 
from rocket to rocket is quite difficult to maintain. 

During flight in the coasting phase, the character of rocket dispersion 
is approximately the same as for ordinary artillery shells. The dispersion 
is basically conditioned by variation in aerodynamic properties from rocket 
to rocket and by eccentricity of aerodynamic forces. In the presence of 
aerodynamic eccentricity the rocket flies with some angle of attack. This 
leads to the appearance of a transverse force deflecting the rocket to one 
side. 

From themathematical point of view, the dispersion can be considered 
as the result of changes introduced into the trajectory parameters by vari- 
ous sources. 

For instance, in the simplest case of a rocket, or an artillery shell, 
projected at an angle to the horizon, in the absence of resistance the range is 


L = (v?/g) sin 24. 


In this case, the range dispersion may be considered the result of scatter 
of Av and Ab in the values of vp and @. Assuming that the changes in the 


320 PRINCIPLES OF STABILIZATION AND STEERING 


values of v and 6 are small, we will get 


AL = (ƏL/ðv)Avo + (8L/86)A% 
or 
AL = (2u0/g) (sin 24,Av0 + Vo COS 20 A). 


It may be seen from the last expression that the deviation in the values 
of vo and 6) have various effects on range dispersion depending on velocity 
v and projection angle ĝo. Specifically, if the firing is conducted under con- 
ditions of maximum range, i.e., when 6) = 45°, the range dispersion is 
determined only by the value Avo. | 

Relative range dispersion for the example under discussion 


AL/L = (2/vo) (Avo + v ctn 20A bo). 


The value AL/L is a minimum when 6 = 45° and becomes infinite 
when f= 0 and @ = 90°, since for these values the range L becomes zero. 

A similar analysis with respect to changing parameters may be per- 
formed for rockets, taking into account the air resistance and certain other 
factors. While doing this it 1s possible to determine not only the qualitative 
but, in certain cases, the quantitative, portion of the problem and approxi- 
mately evaluate the magnitude of dispersion for proposed rockets. 


3. DISPERSION OF GUIDED BALLISTIC ROCKETS 


As has been indicated above, the range dispersion of guided long range 
ballistic rockets is tied in primarily with the errors in the integrator. An- 
other cause of dispersion of rockets with automatic stabilization is the 
transverse drift to which the described gyroscopic devices do not react. 
If the flow around a rocket is asymmetrical, the development of lateral 





Launch To target 


Fic. 9.41. Drift of a rocket equipped with an automatic stabilizer. 


aerodynamic forces is possible during its motion in the dense atmospheric 
region. Lateral force also develops if the line of action of thrust is eccentric 
and is not directed along the axis of the rocket (along the gyroscopic axis). 

Due to the action of transverse force, the rocket is laterally displaced 
(Fig. 9.41). This displacement or, as it is known, drift, will take place 
without rotation of the rocket axis (inasmuch as this rotation is not per- 
mitted by the automatic stabilizer). The skidding may take place in a 


ROCKET DISPERSION 321 


plane of the trajectory as well as in the direction perpendicular to this 
plane. 

Gyroscopic apparatus react only to the angular displacement of the 
rocket and drift is not corrected by them. As a result, noticeable dispersion 
develops. 

The drift may be eliminated only with the aid of special apparatus. 

One of the possible means is the design of an independent integrator of 
lateral accelerations. By means of second order integration of accelerations 
it is possible to determine drift and to transmit an appropriate signal to 
vanes I and III (see Fig. 9.17). Practical implementation of such a measure, 
however, requires somewhat cumbersome gyroscopic apparatus necessary 
primarily to eliminate the effect of earth’s gravity. 


X. Ground Equipment and Launching Devices 


A. Launching Devices of the Field, Anti-aircraft, 
and Aircraft Rockets 


The complex of ground equipment is an independent, and quite impor- 
tant, part of rocket technology. The success or failure of launchings, the 
rate of design development, and the practical application of rockets depend, 
to a considerable degree, on the arrangement of ground servicing, design 
philosophy, and simplicity. 

The ground complex includes all physical and organizational provisions 
for rocket launching, subsequent observation of flight, and guidance, if it 
is done from the ground. The physical provision for launching consists of 
means of transporting the rocket to the launching area, fuel supply (if 
fueling is done prior to launching), fueling, and auxiliary facilities. 

In practice, the problems of physical and organizational provisions for 
rocket launching are solved variously, depending on the arrangement, 
purpose, and design of the rockets. For large, long range rockets, the 
implementation of launching represents a rather complicated problem. Con- 
versely, for small, field rockets, this problem is extremely simple and does 
not present any greater complexity than the firing of an artillery shell. 
For a start, we will dwell on these simpler problems of launching short 
range anti-aircraft and aircraft rockets. 

Transportation of small rockets (weighing up to 1 to 1.5 tons) does not 
present special difficulties. Special crating and ordinary means of trans- 
portation are used to deliver the rocket to the launching area. Larger 
rockets may require special dollies and hoists for their transportation. 
Understandably, special hoisting means are necessary for lifting such a 
rocket. 

The launching of a large anti-aircraft rocket is also done in a different 
manner from the launching of medium-sized rockets. The method of launch- 
ing large anti-aircraft rockets links them with the long range rockets. In 
particular, the launching area of one such anti-aircraft rocket will be exam- 
ined below, in a section dealing with the launching of long range rockets. 

Guides of various designs are used as launchers for powder rockets. 

In firing powder rockets, the guides are assembled in blocks and, for 
mobility, are installed on some means of transportation. The guides are in 
the form of profiled rails, along which the rocket slides on special “‘shoes.”’ 
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The firing is controlled from within the cab of the vehicle. The controller 
has a panel on which he can set any order of firing: one, two, three shots 
or salvo firing. The ignition of the rocket charge is done by means of an 
electrical squib. 

Minethrowers are also transported by automotive means. One of these 
is the German six-barrel minethrower of the Second World War. The mine 
used with this minethrower was shown in Fig. 2.4. Another minethrower is 
characterized by the use of crates as launchers of its jet spin stabilized 
rockets. The crate is attached to the body of the tank. The missiles are 
related to a class of rockets, one of which was shown in Fig. 2.7. 

For aiming powder rockets, as is necessary, for instance, in anti-tank 
warfare, long tubes are used as launchers (bazookas). The rocket is inserted 
into the tube from the rear. 
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Fıc. 10.1. Construction of an anti-tank gun. 


In a short barrel design, the powder charge does not have time to burn 
completely in the tube, and completes combustion during free flight. For 
the protection of the gunner’s face from the gas jet, a jet deflector, in the 
form of a conical mushroom made of fine mesh, which allows seeing the 
target but at the same time impedes the jet stream, is-installed on the for- 
ward part of the barrel. In some cases the protection of the gunner is pro- 
vided by a special mask worn over the face. 

When the gun barrel is somewhat longer, the round is so designed that 
the moving powder charge burns out while the rocket is still in the launch- 
ing tube. The gas jet escapes through the rear end of the tube and the 
necessity for protection of the gunner’s face is obviated. 

Figure 10.1 shows a transverse section of an anti-tank gun with an 
inserted round. This gun is characterized by the fact that the caliber of 
the head section of the round considerably exceeds the caliber of the barrel. 
This allows considerable increase in the power of the round. 

The shell is inserted into the barrel from the front. The empinage, as 
is seen in Fig. 10.1, is made of four thin, flexible plates which are bent by 
hand around the tail of the rocket at the time of loading the gun. 

The launchers of powder aircraft rockets differ very little from launch- 
ers mounted on surface means of transportation. 

Figure 10.2 shows a typical launcher installation located under the 
wings of an aircraft. 
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A more contemporary design is one in which the cell with the powder 
rockets is extended prior to firing from the belly of the jet aircraft. 

Large aircraft rocket missiles and large guided torpedoes are suspended 
from the aircraft in a manner similar to that used for ordinary bombs. 

Let us pause briefly on the subject of launchers for anti-aircraft rockets. 

The launching installations for the unguided anti-aircraft rockets differ 
little from the launching installations for firing against ground targets. 
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Fia. ii cel on their launchers.{(Official U.S. Navy photograph) 





If the firing is done by relatively small rounds, the tubular and rail guides 
are assembled into cells and the firing is conducted by salvos. As opposed 
to field installations, the anti-aircraft installations have a greater rotational 
freedom to facilitate aiming at the target. 

The large unguided rockets are launched singly from installations simi- 
lar to the one shown in Fig. 10.3 (the structure of the rocket on the 
launcher is shown in Fig. 2.24). Prior to launching, the operator adjusts 
the azimuth and elevation angle to correspond to the motion of the target. 

It is not necessary for the launching installations of guided rockets to 
be trainable. The launching of a guided rocket is done, in the majority of 
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cases, in a vertical or almost vertical direction with further direction to the 
target by radio beam. 
In the majority of cases the simplest guide rails are used as launchers. 
The launching of an anti-aircraft barrage rocket is somewhat unique. 
The problem consists of lifting a cable to an altitude of § to 1 mile and its 
suspension at the highest point from a parachute. The cable is packed in 





Fic. 10.3. Anti-aircraft launching installation. 


a box which is placed on the ground, and the rocket, during its climb, takes 
up the cable and pulls it. To prevent entanglement of the cable, the launch- 
ing installation is made extremely simple. In one case it is a rod with the 
rocket mounted on one end. 


B. Hoisting, Transporting, Fueling, and Launching Equipment 
of Long Range Rockets 


The organization for the launching of a long range, especially a long 
range liquid, rocket presents incomparably more complicated problems 
than the launching of short range rockets. 

First of all, long range rockets in their physical dimensions and weights 
are such that they require special transportation and hoisting equip- 
ment. They are equipped with delicate, sensitive instrumentation whose 
more delicate components require a special environment for storage and 
transportation. 

The weight of a long range rocket fully prepared for launching is meas- 
ured in tens of tons. To transport and hoist such a rocket would be ex- 
tremely difficult. Therefore, it is clear that a long range rocket must be 
fueled immediately prior to launching, in take-off attitude. But even if the 
weight of the fuel is not considered, there still remain the several tons of 
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the rocket structure itself and of its warhead. In this case, understandably, 
it is not possible to avoid special hoisting and transporting equipment, 
since this is required by the dimensions of the rocket, and especially by its 
handling as a structure quite sensitive to accidental damage. 

For this reason, in some cases, the rocket is delivered to the launching 
site in parts. There the rocket is assembled in the vertical position. 

After the body of the rocket is elevated, the nose section will be installed 
with the aid of a crane, and the stabilizer will be attached to the lower 
section. Then, the fueling of the rocket will begin. In the majority of cases, 
however, the rocket is delivered to the launching area in an assembled 
form as is done, for example, with the V-2 rocket. 

The rocket is transported by means of a surface carriage (Figs. 10.4- 
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lig. 10.4. A rocket on its carriage (V-2 long range ballistic rocket): 7, brace; 2, jack; 
3, folding platform; 4, erecting rail; and 4, ladder. 


10.6). At the launching site, the launching platform is brought to the 
carriage on a dolly. This platform is mated to the carriage, after which the 
rocket is elevated. 

The carriage rail, with the aid of hydraulic jacks, may be elevated to a 
vertical position together with the rocket (Fig. 10.5). In this manner the 
rocket is placed on the launching platform on its four fins. Stabilizing fins 
and the nozzle exit plane are specially reinforced to take the total weight 
of the rocket and to transmit these forces to the main lower bulkhead of 
the airframe, which is butt-jointed to the tail section. During rocket flight, 
the thrust force is transmitted to the same bulkhead. 

A hydraulic telescoping jack, 2, serves to elevate the carriage rail, 4 
(Figs. 10.4 and 10.5). Since, during elevation, the center of gravity of 
the rocket is displaced relative to the points of support (wheels) backward, 
two pressure pads, 1, are provided in the design of the carriage. They may 
be rotated about the vertical axis, thereby moving the support points 
backward. The support pad, /, is shown in Fig. 10.5 in contact with the 
ground. 

After completion of rocket installation on the platform, the rocket is 
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Fig. 10.6. Rocket ready for transportation. (Official U.S. Air Force photograph: 
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Fra, 10.7. U.S. Air Force Atlas on launching pad. (Official U.S. Air Force photograph) 


freed of the straps tying it down to the rail. Later the, rail is moved slightly 
to the side and, during further preparation for take-off, it serves as a ladder 
along which the servicing personnel may reach the instrumentation and 
warhead sections. Servicing platforms are extended for convenience. One 
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of the platforms is at the level of the motor installation and the other at 
the level of the instrumentation section. 

Figure 10.7 shows a much larger rocket (Convair Atlas) on its per- 
manent launching pad. A series of servicing platforms is clearly visible on 
the mobile Gantry crane. 

Figure 10.8 shows a rocket on the launching platform. A supply cable 
from a mobile electric station for ground testing the guidance system is 
supported by an arm attached to the pad. During 
launching the cable is released and hangs from the arm, 
while the guidance system is switched to the internal 
power. 

Figure 10.9 shows a drawing of the launching 
platform. 

The upper part of the platform consists of two 
rings, 1 and 2. The lower ring rests on four supports 
and has a groove on top filled with ball bearings which 
support the upper ring. Therefore, both rings form 
a structure similar to a thrust ball bearing. The upper 
ring, with an installed rocket, may easily be turned 
to the required angle and the rocket plane I-III (see 
Fig. 9.17) may be made to coincide exactly with the 
plane of the trajectory. The rotation of the ring is 
done by means of a ratchet mechanism, 4, connected 
to the chain, 10. 

Each of the fins is located between supports, 9. 

A flame deflector, 3, in the shape of a conical plate 
is located below the launching platform. The rocket 
jet, during launching of. the rocket, impinges against 
the plate and is dispersed in all directions below the 
platform. 

Tubes 5 and 6, attached to the upper ring, are used 
for installing the ignition system. The cable arm is 
attached to the bracket, 8. 

After the rocket is installed on the pad, its fueling 

Fig. 10.8. Rocket begins. The alcohol is delivered in tank trucks. The 
on a launching plat- liquid oxygen is brought up in large Dewar flasks 
form. . 

mounted on truck trailers. 

The oxygen is supplied through a flexible hose which is covered with a 
layer of white frost. The oxygen vaporizing in the tank escapes through a 
vent tube under the platform. 

After fueling with alcohol and oxygen, the motor installation is charged 
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Fig. 10.9. Launching platform (dimensions in millimeters): 7, upper ring; 2, lower 
ring with four legs; 3, flame deflector; 4, ratchet mechanism for rotation of the upper 
ring; 5, 6, fittings; 7, friction lock; 8, mast support; 9, fin rests, and 10, chain. 
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with hydrogen peroxide, also brought up in a special cistern. Before fueling, 
the hydrogen peroxide is preheated to a temperature of 100° to 105°F. This 
is necessary for the proper generation of steam and assurance of motor 
operation in a given regime. 
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Fie. 10.10. Over-all view of the launching platform with the rocket in firing position. 
(Official U.S. Air Foree photograph) 





Potassium permanganate is poured from portable cans, also in a pre- 
heated state. 

The over-all view of the launching platform, with the Bomare rocket in 
firing position, is shown in Fig. 10.10. A typical schematic of the location 
of servicing machinery during fueling is shown in Fig. 10.11. 
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Fre. 10.11. Schematic of the servicing equipment location in the launching area 
during fueling: 1, cable truck; 2, tank truck with oxygen; 3, oxygen pump; 4, launching 
platform; 5, pumping station; 6, tank trucks with alcohol; 7, pneumatic panel; 8, per- 
oxide heater; 9, vehicle with additional equipment and spare parts; 10, electric supply 
vehicle; 11, motor-generator trailer; 12, carriage; and 13, armored truck. 
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As is seen from the presented schematic, during servicing of a long range 
ballistic rocket many auxiliary pieces of equipment are required. Besides 
the above-mentioned surface carriage, launching platform, and cisterns 
with alcohol, oxygen, and hydrogen peroxide, the launching equipment 
includes a pumping station for oxygen, a coupling with a pump for pumping 
alcohol, a heater for hydrogen peroxide, cans for potassium permanganate, 
a heater for potassium permanganate, a compressor capable of delivering 
up to 3300 psi, apparatus for warming up the motor with hot air, an armored 
vehicle for directing launching, a vehicle with a console for horizontal 
testing, a vehicle with a console for verticle testing, a vehicle with spare 
parts, a mobile electrogenerator, a cable coupler, a transformer coupler, an 
accumulator station, a caterpillar tractor for the carriage, a tractor for the 
launching platform and crew, a tractor for the oxygen cistern, a bus for 
the crew, a vehicle for the crew making measurements, liaison vehicles, and 
a supply truck. Therefore, from the above listing for a V-2 rocket, it is seen 
that the aggregate of ground equipment is quite complicated. 

The purpose of the various units comprising the system is understand- 
able without explanations. It is necessary only to point out, in addition, 
that the complicated system of steering and atuomation aboard the rocket 
is subjected, prior to launching, to a thorough checkout. The first stage of 
the checkout is the so-called horizontal tests, 1.e., tests conducted while 
the rocket is still on the dolly. The vertical, or general, tests of the rocket 
already poised on the platform are more thorough. The operation of valves, 
fuel supply system, steering mechanisms and, in general, the entire auto- 
matic system, are checked out at this point. 

During launching, the electricity and pneumatics are supplied by a 
transportable electric generator and compressor. 

Immediately prior to launching, all auxiliary apparatus, with the excep- 
tion of an armored car which is located at a distance of 500 ft from the 
platform for a V-2 launching, are removed. ‘lhe launching command is 
given from the armored car. 

The launching of a long range rocket from a surface platform is not the 
only possible method. During the Second World War, attempts were made 
to launch a rocket directly from a railroad flat car. 

The railroad flat car is equipped with elevatable service platforms and 
a transportable launching platform. During elevation, the platform is 
placed on an erected rail and thus prepared for launching. 

The described method of rocket launching from a platform may also be 
used for large anti-aircraft rockets. 

Figure 10.12 shows a launching area of an anti-aircraft rocket with a 
servicing and guidance system. 
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Fic. 10.12. Control system for a battery of anti-aircraft rockets; 1, launching platform; 2, range finder; 3, control panel; 


4, computer; 5, search radar; 6, launching platform; and 7, guidance radar. 
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The target is located by radar, which transmits a signal to the range 
finder. When the range to the target is determined the rocket is launched. 
During the initial phase, in the first 6 to 8 sec, the flight is vertical with 
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the internal stabilizing system. Then begins the steering of the rocket 
toward the target, whose range is constantly being determined by the 
range finder. Signals are transmitted from the range finder to the computer, 
which determines the necessary direction of the rocket flight for target 
intercept. The results obtained by the computer are transmitted to the 
steering mechanism and, depending on the location of the rocket and the 
computer data, the necessary signal is transmitted to the antenna. In this 
manner the rocket is given the required direction of flight. 

In conclusion let us pause briefly on the launching equipment of pilot- 
less aircraft. 

One of the methods of launching pilotless aircraft is the use of special 
catapults (Fig. 10.13). The pilotless aircraft is accelerated by means of a 
moving piston located underneath the rails of the catapult and actuated 
by the products of hydrogen peroxide decomposition. 
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Fre. 10.13. Catapult for launching pilotless aircraft. 


This method of solving the problem of launching pilotless aircraft can- 
not be regarded as successful. The catapult becomes complex and un- 
wieldly and, most important, not maneuverable, which, under war condi- 
tions, 1s its most serious defect. Much more suitable is a light launching 
installation with guide rails and acceleration by means of powder rockets. 

As was mentioned in Chapter II, the pilotless aircraft may be launched 
from mother planes. 


C. Testing of Rocket Motors and Rockets 


Each rocket, especially a long range rocket, during its development 
passes through an extensive program of tests of individual components, 
as well as in its entirely. A long range ballistic rocket is too complex and 
ex pensive a piece of apparatus to be built without exhaustive experimental 
development of components and assemblies. 

As in aircraft design, the development of a new prototype is accom- 
panied, first of all, by the testing of a rocket model in the wind tunnel in 
order to determine all of the aerodynamic coefficients. 

In fabricating an experimental prototype, the more critical components 
are fabricated in larger numbers than required and part of them are struc- 


306 GROUND EQUIPMENT AND LAUNCHING DEVICES 


turally tested, resulting in an experimental determination of their limiting 
loadings. 

The steering system is subjected to a detailed check-out. The steering 
mechanisms and gyroscopes are tested on special stands. The operation of 
components under conditions of high vibration levels, which develop in the 
body of the rocket during flight, are especially tested. These vibrations are, 
usually, duplicated on special shake tables. 

A separate, and more responsible, problem is the testing and develop- 
ment of motors. 

As with all other details, the components of the motor, first of all, are 
subjected to thorough control during fabrication and assembly. This con- 
trol consists of dimensional inspection and then testing of connections, 
piping, chamber components, and fuel tanks for pressure tightness. 

An assembled motor, with the delivery system, is subjected to pneu- 
matic and hydrostatic tests. During pneumatic tests, the system is filled 
with an excess pressure of 14-28 psi of air. All connections and welded 
joints are covered by brush with a solution of soap. Air leaks through the 
seams are immediately spotted by the formation of soap bubbles. During 
hydrostatic testing, the system is subjected to the full operating, and some- 
times exceeding the operating, pressures. The appearance of leaks is moni- 
tored by the pressure drop within the system and by leakage at the joints 
and seams. 

The motor, as a unit, is subjected (just as all the other rocket assemblies 
are) to functional tests. The proper functioning of the separate sections of 
the motor, the degree of proper functioning of the automatic system, the 
performance of the valves, and the performance and reliability of the liquid 
components supply are determined during this test. 

After a full, thorough check-out, static firing tests are begun. The motor 
is installed on a special stand and is started. The performance of the operat- 
ing motor is measured and the results compared with the calculations. 
Design changes are made if required. 

The firing stands are not used solely for the design development of the 
motors. Investigation of processes within the liquid fueled rockets requires 
extensive experimentation with models (small motors), and with full size 
motors as well. The same stands, with corresponding changes in the pro- 
gram of measurements, are used for this purpose. 

Motors of already developed production rockets are also subjected to 
static tests for the purpose of determining (or refining) the basic rated 
characteristics of the motors. 

The experience of many years shows that, contrary to a widespread 
belief, the testing of rocket motors does not present serious hazards to the 
servicing personnel. However, in conducting static tests, a whole series of 
safety measures must be observed. 
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First of all, caution is necessary in handling fuel components, many of 
which are caustic and are harmful to the human organism. Among these 
are nitric acid and hydrogen peroxide. The handling of liquid oxygen also 
presents a certain danger. Therefore, the preparation for the test itself 
requires observance of definite safety rules such as, for instance, the use of 
special garb for personnel taking part in the charging of fuel and auxiliary 
tanks. 

During improper operation of the motor there is a possibility of combus- 
tion chamber explosion. Such explosions are very powerful, are accom- 
panied by extensive destruction, and quite often cause fire. Therefore, the 
test bay, i.e., the location of the motor during tests, is usually made semi- 
enclosed, surrounded on three sides by concrete walls in order to reduce the 
effect of fragmentation and to localize the center of the fire resulting from 
the explosion. The servicing personnel are located in shelters some distance 
from the motor. The starting and stopping of the motor are remotely 
controlled. 

The administration of static tests also requires the introduction of safety 
measures. The servicing personnel are thoroughly indoctrinated in rules of 
conduct and in the sequence of operations. Immediately prior to start, a 
general signal is given (usually a siren), warning of the beginning of the 
static firing. 

Test stands for large rocket motors, as a rule, are remotely located 
with respect to populated areas. This is necessary if only because the jet 
exhaust of a large rocket motor is accompanied by loud noise. A full descrip- 
tion of the noise of an operating rocket motor is difficult to convey in words. 
Suffice it to say that the noise of an operating 25-ton liquid rocket motor 
many times exceeds the noise of the propellers of a powerful airplane and 
at a short distance is tolerated with difficulty by man. 

The following basic parameters of an operating motor are determined 
and measured on the stand: 

1. Thrust. 

2. Consumption of: (a) combustible; (b) oxidizer; (c) auxiliary com- 
ponents. | 

3. Pressures: (a) in the combustion chamber; (b) ahead of injectors; 
(c) at the various points of the intake and delivery lines of the fuel system. 

4. Temperatures: (a) of the gas in the combustion chamber and in the 
nozzle; (b) of chamber walls; (c) of the cooling liquid. 

5. Velocity of exhaust gases. 

6. Chemical composition of the combustion products. 

7. Various parameters of the tubropump assembly or of the displace- 
ment delivery system. 

The thrust, component consumption, and pressures are determined 
with a satisfactory degree of accuracy. The thrust is usually measured by a 
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hydraulic dynamometer; the consumption, with the aid of flowmeters. 
The pressure is measured by manometers. 

The gas and wall temperatures are determined with a considerably 
lesser degree of accuracy. The chamber wall temperature is usually meas- 
ured by means of thermocouples imbedded in the wall in advance. The 
accuracy of these measurements is not great, inasmuch as a considerable 
temperature gradient is established across the wall and the dimensions of 
the thermocouple are close to the thickness of the wall. Even greater diffi- 
culties are present in measuring the temperature in the chamber and the 
temperature of the high velocity jet in the motor nozzle. In these cases, 
approximate temperature measurements obtained with the aid of an optical 
pyrometer must be satisfactory. 

The exhaust velocity of gases varies across the nozzle exit plane and its 
measurement at various points with the existing instrumentation is ex- 
tremely difficult. Usually, one must be satisfied with the determination of 
the average effective exhaust velocity obtained from thrust and fuel con- 
sumption data. 

The static tests yield the most extensive and fruitful data for the final 
development of the motor. 

After the motor has been developed it is possible to begin the testing of 
the rocket as a unit. These tests consist, first of all, in static tests of a rocket 
with an operating motor, for the development of the automatic system in 
the environment approximating operating conditions. 

The next step is the flight tests of the rocket. This stage is more serious. 

The first launchings of military rockets are conducted, naturally, with- 
out a warhead. Instead, recording and radio transmitting instrumentation 
are installed in the prototypes, and a system of end instruments for the 
measuring of the most important flight parameters is located at various 
points on the rocket. 

The measurements made in flight are: 

1. Operating parameters of the motor: (a) chamber pressure; (b) com- 
ponent consumption; (c) parameters of the delivery system (number of 
turbine revolutions, tank pressure above the component surface). 

2. Component level in tanks. 

. Temperature at various points on the rocket. 

. Axial and lateral accelerations. 

. Aerodynamic pressures at various points on the fuselage. 

. Loads and stresses at the more critical structural components. 
. Amplitude and frequency of developed vibrations. 

8. Guidance system parameters: (a) angle of rocket rotation about 
gyroscopic axes; (b) vane rotation angles; (c) the magnitude of steering 
signals at certain intermediate sections of the automatic stabilizer; etc. 
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For the measurement of the indicated quantities a multitude of end 
instruments are used; specifically, pressure transducers with different 
ranges; level indicators; temperature transducers whose operation is based 
on change of electrical resistance with temperature. Slide wire potentiom- 
eters which generate voltage proportional to the angles of rotation are 
used as end instruments for the angle of vane rotation and of the rocket 
as a whole. Accelerometers (of the type consisting of a mass supported by a 
flexible element) are used for recording of accelerations and vibrations. 

All mechanical end instruments (pressure, acceleration, level) are 
equipped with slide wire or potentiometer devices, which transform end 
instrument indications into electrical signals which are transmitted to the 
ground receiving station through radio channels. 

In this manner, as a result of flight test, it becomes possible to deter- 
mine the variation of the rocket basic flight parameters with time. 

A measurement system of this type is called a telemetering system, and 
represents one of the most important branches of contemporary rocket 
technology. At the present time it is impossible to conceive of a develop- 
ment of a large rocket without the use of a telemetering system. On the 
basis of telemetered data the experimenters obtain refined information on 
flight conditions, and in the case of an unsuccessful launching, by analyzing 
the received data, they can determine the reasons for this failure and elim- 
inate them in the future. 

During flight tests the observation of flight may also be done by means 
of a cine-theodolite. The rocket is observed up to the limit of visibility and 
its coordinates with respect to time are registered simultaneously with 
photographing the rocket on a movie film. The data supplied by cine- 
theodolite photography, supplementing the telemetering data, are quite 
important for the design development of a rocket. In the absence of direct 
visibility the observation of the rocket along the trajectory is performed 
with the aid of electronic tracking devices. 

The launching of large rockets is a very important undertaking. We 
already know that the rocket, because of its high velocity, even though 
not equipped with a warhead, can cause a great deal of damage on impact. 
During experimental launchings, failures are always possible. The malfunc- 
tioning of the motor or guidance system will result in rocket impact either 
in the vicinity of the launching area or with considerable deviation from 
the target. Finally, it is possible that the rocket will leave the predeter- 
mined course completely. In such cases populated areas may suffer. Not 
only experimental but operational launchings of the V-2 rocket in Germany 
serve as an example of this. It is known, for instance, that during the 
bombardment of London by the already developed at that time V-2 rockets, 
from the suburbs of Hague, out of a thousand launchings 80 were unsuccess- 
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ful, while five of the rocket fell in Hague, causing great destruction in the 
city. 


From everything that has been said above it must be concluded that the 


testing of long range and sounding rockets must be done in deserted areas, 
and the basic path of the rocket must lay over unpopulated regions. Finally, 
rocket launchings themselves may be undertaken only after a long prepara- 
tion and a thorough and most careful check-out of the equipment. 
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Axial compressor motors, 67 


Benzine, 35, 36, 70 
Bomarc, 331 
Boosters, launching, 80 
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Catapults, 335 
Center of gravity, 236-238, 275 
Center of pressure, 236-238, 275 
Centrifugal compressor motors, 67 
Cinetheodolite photography, 339 
Circuits 
amplifier, 308-309 
converter, 306-307 
differentiating, 304-306 
Combustion 
basic equation of, 153-154 
effects of temperature and pressure on, 
157 
pulsating, 139-140 
solid fuels, 140-146 
Combustion chamber processes, 127-160 
in liquid-fuel motors, 127—140 
products of combustion, 147—154 
in solid-fuel motors, 140-146 
Combustion chambers, design, 133-138 
Combustion products, parameters of, 
147-148 
Congreve, Col., 4 
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Cooling, design of in motors, 197-199 
Coridis acceleration, 239-240, 285-286 
Critical velocity in gas flow, 174-175 
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Delta wings, 282 
Demodulator, see Circuits, amplifier 
Design index, 17 
Diffusers, multishock, 70 
Discharge, equation of, 166 
Dispersion 
in guided rockets, 320-321 
measure of, 315-316 
in unguided rockets, 317-320 


Efficiency of rocket motors, 23-24 
Energy, equation of, 167 

Eno-Peltri, R., 9 

Equilibrium constant, 155 
Equipment, auxiliary, 333 
Euler-D’Alembert paradox, 218 
Euler's equation, 282 

Exhaust velocity, 172 
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Flame holders, 69 
Flame-out, 69 
Flight, forces affecting, 200-206 
aerodynamic forces, 202-203, 212- 
241 
ballistic coordinates, 200-201 
equation of motion, 205-206 
frontal resistance (or drag), 202, 
217-221 
inertial forces, 203-205 
stabilizing and damping moments, 
235-241 
steering force, 203, 241-245 
thrust and weight, 245-248 
Fluoride combustion products, 101 
Fluorine, 100 
Fuel, see also Specific impulse 
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Fuel displacement by powder charge, 
76 
Fuel, method of feeding, 73-77 
displacement, 73-76 
pump systems, 76-77 
Fueling, 329-332 
Fuel systems, automation of, 90 
Fuels, 26, 28, 35, 40, 42, 51, 59, 62, 
73, 79 
basic requirements, 96—99 
calorific value, 95—96 
chemical energy of, 91-95 
classification tables, 105, 106 
combustion of, 91-95 
experimental, 117-118 
gas constant, 98 
heat of formation, 93, 102 
nuclear energy, 91, 119-126 
solid, 104, 114-116 
specific heat, 98-99 
specific volume, 98 
table of properties, 100 
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Gas flow velocity, 173 
Gauss law of random distribution, 316 
Glushko, V. P., 8 
Goddard, R., 9 
Graphite vanes, 244 
Gravity, influence of, 21-22 
Ground equipment, 26, 322-323, 325- 
335 
Guidance systems, 26, 32-33, 43, 49, 50, 
58, 281 
canard type, 50 
homing, 49, 58 
Gunpowder, use in rockets, 3—4 
Gyroscope, 282-287 
application, 287—294 
correction of, 288-289 
horizontal free, 291-293 
integrator, 312-315 
properties, 282—287 
vertical free, 293-296 
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Heat absorption, 103 

Heat capacity of gases, 150-151 

Heat content of gases, 151-153 

Heat exchange in rocket motors, 193-197 
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by convection, 193-195 
by radiation, 195-196 
Heating, aerodynamic, 233-235 
Hoisting equipment, 325-326 
Hydrazine hydrate, 86 
Hydrogen, 17 
Hydrogen peroxide, 38, 42, 58, 73, 83, 
86, 112-114 
Horizon, artificial, 288-289 


Ignition, 77-78, 138-139 

electric spark, 78 

pyrotechnic, 78 

self-igniting fuels, 84 

solid fuels, 145-146 

time delay in, 103-104 
Injection heads, design, 136-138 
Integrator, see Gyroscope 
Internal energy of gases, 148-151 
Interplanetary flight, 269-270 
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Jamming, of guidance systems, 281 
Jet-propelled, pilotless aircraft, 25-61 
advantages, 40 
disadvantages, 40 
Jet stream 
effect on aerodynamic forces, 232-233 
Jet vanes, 241-245 
Jets, centrifugal, 128-130 
injection, 127-128 
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Kerosene, 35, 73, 79, 110, 111 
Kondratiuk, U. V., 8 
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Launching devices, 26-27, 59, 322-335 
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Mach angle, 227 

Mach cone, 223 

Mach numbers, 171, 214, 215, 216, 229- 
232, 237, 239-240 

Margin of stability, 237 

Mass rate of discharge, 11 

Matador, 39 

Mescherski Equation, 10 

Michailov, Onisim, 5 
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Missiles, bombardment 
intermediate range, 35-36 
long-range, 33—36 
short-range, 27-33 
Modulator, see Circuits, converter 
Moon rockets, 270-271 
Motion, equations of, 250, 259, 296 
Motors, auxiliary, 271 
reaction, 62—90 
air-breathing, 1 
ducted, 68-70 
energy conversion in, 62-64 
non-air-breathing, 2 
pulse jet, 70-72 
table of, liquid-propelled, 72 
types of, 64-65 
rocket 
altitude characteristics of, 190-191 
throttling characteristics of, 191-193 
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Navigation, celestial, 46 
Newton, law of momentum, 10 
Nitric acid, 108 
Nitric oxides, 109 
Nitrogen oxide, 35, 73 
Nozzle expansion ratio, 181 
Nozzle, Laval, 176, 229-230 
Nozzle of rocket motor 
cooling of, 193-197 
gas flow through, 161-179, 184-185 
performance of, 179-187 
properties of, 187-193 
Nozzles, supersonic, 175-179, 187-190 
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Oberth, H., 9 
Oscillations, 275—277 
phugoid, 276 
Oxidizers, 73, 100 
classification table, 105 
Oxygen, 79 
combustion products, 101 
liquid, 17, 35, 73, 80, 111-112 
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Perelman, J. I., 7 
Peter I, 5 

Pitch angle, 201 
Pitching, 281 
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Powder grains, shapes of, 142-144 
Precession, law of, 284-286 
Projectiles, rocket-sustained, 35 
Proximity fuse, 50, 58 

Pulsejet, 36 


Ramjet, 35, 36, 38, 45 
Range, formulas for, 266—269 
Range control, 309-315 
Reactive motion 
definition, 1 
theory of, 10-24 
References, 340 
Regulus missile, 55 
Residual time, in combustion, 132 
Reynold’s number, 216 
Rocket technology, history, 3—9 
Rockets, aircraft, 54-59 
anti-aircraft, 46-54, 79 
anti-tank, 28-29 
ballistic, 40-61 
advantages, 44 
disadvantages, 44 
intermediate range, 40-44, 80-84 
long range, 45-46 
interplanetary, 79 
motors, 79 
meteorological (or sounding), 59 
underwater, 59 
Roll, 277, 281, 293 
Roll angle, 202 
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Satellite, earth, 269-270 
Shaped charges, 28-29 
Shock waves, 223-229 
Skip bombs, 59 
Sound, velocity of, in gases, 168—172 
Specific impulse, 15 
Stabilization, 30-81, 273-289 
by automatic pilot, 278-282 
feedback mechanisms, 280-281 
spin, 31, 277, 287, 319 
Stabilizers, intermediate devices, 301-309 
Stages, separation of, 20 
Static stability, 275 
Steering, see Stabilization 
Steering mechanisms, 203, 251, 281, 289- 
291, 294-296 
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Streamlined shapes, 219, 228 
Submarine, missile launching from, 38 
Supersonic flow theory, 221-229 
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Talos missile, 52 
Tangential acceleration, 204, 265 
Telemetering system, 339 
Terrier missile, 54 
Testing, motors and rockets, 335-340 
safety measures, 336—337 
Tetranitromethane, 109 
Thermal dissociation, 154—160 
Thermodynamic process equation, 162 
Thrust, 11-16, 67, 69, 79, see also Flight 
forces affecting 
effective gas velocity, 15 
equation of, 13 
regulation of, 180 
relation to altitude, 14 
Tihonravov, M. K., 8 
Trajectory, in rocket flight, 249-272 
coasting phase, 249-250, 258 
of earth satellite, 260-262 
flight beyond the atmosphere, 259-262 
powered phase, 251-259 
Transportation equipment, 325-326 
Tsander, F. A., 8, 17 
Tsiolkovski, C. E., 6, 17, 96, 272 
Tsiolkovski equation, 16-17 
Turbojets, 36, 40, 65-68 
Turboprop engine, 3 
Turbopump, 77-78 
Turbulent flow, in combustion chamber, 
130-131 
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V-1 rocket, 36-87, 71 
V-2 rocket, 9, 17, 22, 40-43, 239 
auxiliary equipment, 333 
launching, 43 
launching history, 339-340 
motor chamber, 135-137 
power plant, 80-84 
stabilizing system, 301 
steering, 289-291 
thrust, 80 
trajectory, 254 
transportation, 326 
velocity, 43 
Valve, pressure reducing, 74 
Vanes, see Jet vanes, Graphite vanes 
Velocity, escape, 264-265 
ideal, 16-21 
of single-staage rocket, 16-18 
of multi-stage rocket, 18-21 
Velocities, subsonic, 217 
supersonic, 214, 217 
Viscosity, of gases, 164 
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Warhead, 26, 29, 36, 42, 51, 58 
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World War II, Soviet rockets in, 9 
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Yaw angle, 201, 277 
Yawing, 279, 280, 281 
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